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Philip S. Carter 


1897-1961 


On April 20, 1961, the field of antennas lost one 
of its senior practitioners, and the Professional 
Group on Antennas and Propagation one of its 
elder statesmen. On that day Philip S. Carter died 
at his home in Port Jefferson, New York, after an 
extended illness. Phil Carter has been sorely missed 
at the PGAP meetings, where he long has con- 
tributed to the solution of major problems and to 
the fun. 

He was born in Glastonbury, Connecticut, on 
July 22, 1897, and graduated in Mechanical Engi- 
neering from Stanford University, California, in 
1918. He served in the Signal Corps in World War 
I, was employed by the General Electric Company, 
and in 1920 transferred to RCA, where he remained. 
He worked during these years with Alexanderson, 


Beverage and Hansell. 

He contributed to the development of the Wave 
Antenna, participated in the invention of many 
directive antennas using long wires, and himself 
invented the very popular folded-dipole antenna. 
In the course of this work he obtained 100 patents. 

In his many published papers in the IRE Pro- 
CEEDINGS, he established an international reputa- 
tion as an antenna expert and a brilliant analyst. 
He served on many IRE Committees, including 
the Antenna and Waveguides Technical Commit- 
tee and the Administrative Committee of PGAP. 

Those of us who worked so closely with Phil 
Carter on PGAP problems share our loss with his 
many other associates in the IRE and at RCA. 

Li CoN 
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contributions 


Dipole Antennas Coupled Electromagnetically 
to a Two-Wire Transmission Line* 


K.-M. CHEN? anp R. W. P. KING{, FELLOW, IRE 


‘Summary—The circuit properties of a dipole antenna coupled 
electromagnetically to a two-wire transmission line are studied 
theoretically and experimentally. Some arrays constructed of anum- 


- ber of dipole antennas coupled electromagnetically to a two-wire line 


are considered and the radiation patterns are found. 


INTRODUCTION 


N antenna array consisting of dipole antennas 
' coupled electromagnetically to a single two- 


wire transmission line was first proposed by 
Sletten [1]. Since then, experimental studies of this 


array have been made by several authors [2], but their 
_ results are neither comprehensive nor entirely consist- 
ent. The research reported in this paper is intended to 


provide a reasonably complete theoretical and experi- 
mental picture of the properties of such arrays. 
The array consists of several dipole antennas situated 


- in a plane parallel to that of a two-wire transmission 
line. The centers of these dipole antennas are located 


opposite positions of voltage maximum along the trans- 
mission line. When the axis of a typical dipole antenna 


* Received by the PGAP, November 3, 1960; revised manuscript 
received, February 27, 1961. This research was supported by the 
U. S. Air Force under Contract No. AF19(604)4118. , 

+ Rad. Lab., University of Michigan, Ann Arbor, Mich. For- 
merly at Gordon McKay Lab., Harvard University, Cambridge, 


Mass. : 
t Gordon McKay Lab., Harvard University, Cambridge, Mass. 


is parallel to the direction of the line, no currents are 
excited, since the antenna is in the neutral plane of the 
transmission line. However, when the antenna is ro- 
tated from that neutral position, the antenna is excited 
by the charges on and currents in the line. The magni- 
tude of the excitation which induces currents in each 
antenna varies with 1) the angle of rotation of the an- 
tenna from its neutral position, 2) the length of the 
antenna, 3) the distance between the antenna and the 
line, and 4) the spacing of the transmission line. It is the 
dependence of the antenna current on these many 
parameters that makes the problem both difficult and 
interesting. 

Some features of an antenna system of this type are: 


1) In constructing a high-gain antenna array for 
high-power VHF-UHF communication purposes 
in a conventional manner, numerous transmission 
lines must be used to feed the antenna elements 
with voltages of proper amplitude and phase. In 
such an arrangement the transmission lines are 
easily unbalanced owing to interactions between 
them. This causes difficulty in adjusting the feed- 
ing system and makes the realization of a prede- 
termined field pattern difficult. The proposed an- 
tenna system sidesteps this difficulty by having 
only a single two-wire transmission line. 
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2) An antenna array of the proposed type can be de- 
signed to radiate more power than a slot array fed 
by a waveguide. The simplicity of the construc- 
tion is another convenience. 


In the analysis the antenna current induced by the 
near-zone field of the two-wire line is obtained first, and 
then the equivalent input impedance of the antenna as 
seen from the line. In the last part of the paper, several 
antenna arrays with different properties are suggested 
and the corresponding radiation patterns are obtained. 
The theoretical results are compared with experiment 
and found to be in good agreement. 


ANTENNA CURRENT! 


The first topic to be investigated is the current in- 
duced in a dipole antenna when it is coupled to a two- 
wire line. Before the dipole is coupled to the line, which 
is short-circuited at one end, the distributions of charge 
per unit length and current along the line havea stand- 
ing-wave pattern. After the dipole is coupled to the line, 
the wave pattern on the line tends to change from a 
standing wave to a traveling wave since the line trans- 
fers power to the antenna. 

In general, the distribution of the charge per unit 
length on the line can be written as 


g(z) = [qi cos Boz — jq2 sin Boze? (1) 


From the equation of continuity of charge, the corre- 
sponding distribution of current in the line is expressed 
as 


T(z) = [—jvogi sin Boz + v0g2 cos Bozle%*. (2) 


For the TEM mode, these charges and currents main- 
tain a transverse electric field in the transverse plane 
perpendicular to the axis of the line. If the geometry of 
the assembly of the dipole antenna and the line is as 
shown in Fig. 1, the tangential electric field at the sur- 
face of the antenna can be found [3] to be 


= [b? + d? — s* sin? 6] 
E,(s) = : : 
reo [(ssin9 + d)? + 8?][(s sind — d)? + 3? 
X [gi cos (Bos cos 8) sin 6 — jqz sin (Bos cos 8) sin 6], (3) 


where s is the distance measured from the center along 
the dipole antenna, @ is the angle which the antenna 
made with the axis of the line, b is the distance between 
the dipole antenna and the line, and 2d is the spacing of 
the two-wire line. Also shown in Fig. 1 are the charge 
and the current distributions when g.=0 and the con- 
figuration of the transverse electric field. 

E.(s) may be resolved into symmetrical and _ anti- 
symmetrical components with respect to s =0 as follows: 


‘The details of this part are available in a previous work of 
Chen [3]. 


ae me 
a 


Transmission |) 
Line 
Dipole 


(b) 


Fig. 1—Geometry of antenna assembly (transmission line 


has a standing wave pattern). 


fae —dq [b2 + d? — s? sin? 6| 
re, [CS Stn Boe Ne |e 
-cos (Bos cos @) sin 6; (4) 
id b? + d? — s* sin? 6 
Eg(s) = —— ! 


wey [(ssin@ + d)? + 6?|[(s sing — d)? + 82] 
-sin (Gos cos 8) sin @. (5) 


E,{(s) induces a symmetrical current, E,*(s) an anti- 
symmetrical current in the antenna. In the following | 
analysis, the symmetrical and antisymmetrical parts of © 
the antenna current are found separately. 

Since the tangential component of the electric field at — 
the surface of the perfectly-conducting antenna van- 
ishes, a differential equation governing the symmetrical — 
vector potential at the surface and in the direction of 
the antenna as maintained by the symmetrical current 
in the antenna [4] is 


0° A,? 
Os? 


5 oid Sen (6) 
(63) 


A,* is defined as 


A a*(s) 5: ns 
Tv 


T.*(s")ds’, 


ip e-B0./(5 — 5/)2 + a? 
- = V(s—s') +@ 


where /: is the half length and a is the radius of the 
dipole antenna. The general solution for A,*(s) is found 
with symmetry condition A,*(s) =A,*(—s). It is 


A,*(s) = oe ge cos Bos — 6.(s) |, (2). 3 


_ ud sin °f [b? + d? — w? sin? 6] 
Teo o [(wsind+d)?+82][(w sin 6 — d)? + B2| 
Cos (Bow cos 8) sin Bo(s — w)dw. (8) 


a eondition J,*(s= + hy=0. 


In the antenna of length 2h, with Boh near 1/2, it is 


4 reasonable to assume the zeroth-order distribution of 


"e 


BE 


a mrithe cg = 


- 
a 
‘ 


_ the symmetrical current to be given by 


[7.*(s) lo = Ji eed eu Bos — cos Boh). 


(9) 
; 0Ys 
; © The substitution of (9) in (7) with s=h leads to 
A,*(h) = as [C, cos Bok — 0,(h)| 
Vo 
e780 h— ss)? 2 
=f C,(cos Bos — cos Boh) wera 
Sait V(h a8 s)F Fig a? 
After rearrangement, C, is found to be 
s00(h 
ae W04(1) Pao 


Cos Boh. ~~ Ca(h, h) ais cos BohE.(h, h) 


_ The expansion parameter y, is taken as the ratio of the 
_ maximum vector potential to the maximum current at 
_ s=0on the surface of the antenna. It is found to be 


aH, 0) — hE,(h, 0 
a ee) (11) 
1 — cos Boh 


_ The constants which appear above are defined as 


eee Go) at 
Cally w) = fo oe 


-cos Bosds, 
h g-iBox/(s — w)? + a 
Eth-%) ail CAS ele ins ds, 
ee AL 6) oo 


1207 ohms. 
In the case By =2/2, which is of particular interest, 
the symmetrical antenna current is 


jAr 
foCa(h, h) 


6,(2) cos Bos. (12) 


[75*(s) ]o = 


; Note that the first-order current is available in a pre- 
vious. work of Chen [3]. 


Following the method used in the analysis of the 
symmetrical current, the antisymmetrical current is 


found as follows: 

The antisymmetrical vector potential at the surface 
and in the direction of the antenna as maintained by the 
ee ee current in one antenna is 
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4Q7 
2A, Bo? 
+ fo A= =f —— ents): (13) 
Os? 
With the symmetry condition A,*(s) = —A,*(—s), the 
general solution of A,*(s) is found to be 
A (ss) = = [C. sin Bos — 6.(s)], (14) 
i) 
where 
6,(s) = -f E(w) sin Bo(s — w)dw 
0 
—jged sin 6 [b? + d? — w® sin? 6] 
Teo 0 [(wsiné + d)? + b?|[(w sing — d)? + 52] 


‘sin (Gow cos @) sin Bo(s — w)dw. 


(15) 


From the boundary condition J,*(s = +h) =0, the anti- 
symmetrical current in an antenna with B)h=7/2 can 
be found to be 


[7.2(s) ]o = 


* [6a(h) sin Bos — 64(s)], (16) 


i 


where 


1 
Yo = 


(aati) sin Bo 
h\ Ja 
0.7070.(h) — 0, (=) 


— 04(s)| K (> os as 


The resultant current in the dipole antenna with 
Boh =1/2 is obtained by adding up the symmetrical and 
the antisymmetrical currents as follows: 


[7;(s) Jo = [7.*(s) Jo te [7.2(s) ]o 
_ ey 0 
0 Calh, h) 


Gone = a(9) (17) 
Wy a\le 0 a . 


a 


cos Bos 


The ratio of the amplitudes of [Z,*(s) ]o and [J,*(s) Jo 
is 


8,(h) 
iii ee HY On ee 
[Z.2(s)]omax _0.70784() 5 
Va 


It can be shown that [J,°(s) ]o and [J,*(s) ]o are approxi- 
mately 90° out of phase. The magnitude is 


iat | 183) A |r 2(5) |? 
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A dipole antenna with 6)4=7/2 carries both sym- 
metrical and antisymmetrical currents, but the latter is 
usually much smaller than the former. Therefore, in 
practice, only the symmetrical current needs to be con- 
sidered. 

Fig. 2 shows the theoretical and the experimental dis- 
tributions of current along the antenna. The latter were 
measured with a small loop protruding from a slot in the 
antenna. The distribution is very close to a cosine curve. 
Figs. 3 and 4 show the theoretical and the experimental 
variations of the antenna current at the center of the 
antenna with the rotation of the antenna. In Fig. 3, the 
distance between the antenna and the line is 6.5 cm and 
the distribution pattern on the line is essentially a 
standing wave. In Fig. 4, the distance from the antenna 
to the line is 0.75 cm and the distribution of current 
along the line is maintained as a traveling wave in the 
course of the antenna rotation by loading the line with 
another antenna coupled to the line a wavelength away. 


Relative amplitude of antenna current 
to) 


2 3 4 
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The optimum orientation of the antenna at which the 
maximum antenna current is obtained depends greatly 
on the distance between the antenna and the line. Figs. 
5 and 6 show the theoretical and the experimental varia- 
tions of the antenna current with the length of the an- 
tenna. The resonant length at which the antenna current | 
reaches its maximum is found to be slightly shorter than 
a half wavelength [5]. Figs. 7 and 8 show the theoretical 
and experimental variations of the antenna current with 


the distance between the antenna and the line. It is in- 


teresting to observe that when the orientation of the 
antenna is great, the antenna current decreases quite 
slowly when the separation is increased. Fig. 9 shows the 
theoretical and the experimental variations of the an- 
tenna current with the spacing of the line. The antenna 
current is increased as the spacing of the line is in- 
creased. It is to be noted that the power loss due to 
radiation from the line is also increased by an increase 
in the line spacing. 


o'* ORDER THEORY } 
EXPERIMENT 
o'” ORDER THEORY } 


6 = 23.6° 


EXPERIMENT AES 


o' ORDER THEORY 


EXPERIMENT US 


Ke 
226m 


‘s 
+e 


5 6 G 8 3 10 


Distance from center of antenna s in cm 


Fig. 2—Current distribution along the dipole antenna. 
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Fig. 3—Variation of antenna current I(s=0) with 
the rotation of the antenna. 
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Fig. 4—Variation of antenna current I(s=0) with 
the rotation of the antenna. 
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—o— Theory 


~------ Experiment 
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Fig. 5—Variation of antenna current J(s=0) with 


the length of the antenna. 
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Fig. 7—Variation of antenna current with the separation distance 


when the voltage on the line is kept constant. 
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Fig. 6—Variation of antenna current J(s=0) with 
the length of the antenna. 
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Fig. 8—Variation of antenna current with the separation distance 


when the voltage on the line is kept constant. 
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Fig. 9—Variation of antenna current with the spacing 
of the transmission line. 
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An important problem relating to the antenna cur- 
rent is the significance and effect of the unbalanced cur- 
rent in the line. Since the discussion is rather involved, 
it is omitted here. It may be found in a previous work of 
Chen [3]. 


EQUIVALENT INPUT CONDUCTANCE OF A DIPOLE AS A 
LOAD FOR THE LINE 


When a dipole antenna is coupled to a two-wire line 
with its center opposite a voltage maximum along the 
line, it is equivalent to a shunt admittance across the 
line at that voltage maximum. If the line is assumed to 
be lossless, the power dissipated in that equivalent 
admittance is 


P= Re ($VI*) = 3G, (18) 
where V is the line voltage across the equivalent ad- 
mittance and G is the equivalent conductance. This 
power dissipated in the equivalent admittance must be 


equal to the power radiated from the dipole antenna. 
That is, 


P= 31 0"Ro, (19) 
where J) is the input current and R, is the input re- 
sistance [6] of a dipole antenna. 

I) is found in (9) as a function of g;. V can be related 
to gi as ' 


2d 


In — 


q1 
Cc TTEQ 


C is the capacitance per unit length of a two-wire line 
and R is the radius of the line wire. 

When (18) and (19) are set equal to each other, an 
expression for G is obtained as follows: 


(20) 


where 
g 4(1 — cos Boh) 
[v. cos Boh — Cu(h, h) + cos BohEa(h, h)| 


h 2 aes 2 a] 
x dsing f [b> +d w? sin? 6| 
o [(wsin 6 + d)?+ b2][(w sin @ — d)? + 6? | 


Cos (Bow cos ) sin Bo(h — w)dw. 


The theoretical results are compared to the experi- 
mental results in [1]. The agreement is quite satisfactory. 
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RADIATION PATTERN? 


The fundamental unit of the antenna array described 


in this paper is a pair of tilted dipoles as shown in Fig. | 


10. In a dipole antenna with B)4=7/2, the current can 
be assumed to be given approximately by 


I(s) = Io cos Bos + jla sin 266s, (21) 
where 
jar 
On 6,(h), 
FoCalh, h) 
—2r 
1, = 6k) < — Lo 
CoWa 


@=90° PLANE 


Fig. 10—Coordinate system for determination 
of radiation pattern. 


With this antenna current the radiation pattern of a 
pair of tilted dipoles may be evaluated. In the interest 
of simplicity and practical design, only the radiation 
patterns in the £ plane (or 6=0° plane) and in the H 
plane (or 6=90° plane) are presented here. 

Let the pair of dipoles be located in the plane ¢ =0° 
with the tilt angle 6,. The radiation function in the E 
plane maintained by the symmetrical antenna current 
is found to be 


Tv 
cos E cos (0 — 64) 


F,(6, 6.) = e 7 (7/2) cos 8 
sin (06 — @,) 
us 
cos E cos (0 + 6.) 
ef r/2) cos ae (22) 
sin (6 + @,) 
and 
—Jeolo |e ioRo 
Eo"(6, 0.) = F,(8, 02). (23) 


0 


* The details of this part are available in a previous work of 
Chen [7]. 
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The antisymmetrical antenna current’ maintains a 


radiation field as follows: 


Tv 
cos =. cos (8 — @,) 


F,(8, 84) = sin (6 — 6,)e7i(#/2) cos 8 
[4 — cos? (6 — 6,)| ( ye 


T 
cos |= cos (8 + 6.) 


sin (@ + 6, ef (7/2) cos y 
[4 — cos? (8+ 6) | ( ) 
and 


Cols e iBoRo 


us 0 


Eo,"(6, 8) = F.(6, 8). 

However, since InIo, Eo,(0, 0.) may be ignored so 
that £,"(6, 6.) alone, as given in (23), is adequate in 
practice. 

From (22), it is observed that the E plane radiation 
pattern has four major lobes; two along the line 6=90° 
(the broadside effect), the other two along the line 
6=0° (the end-fire effect). Therefore, it is clear that this 
arrangement of antennas functions as a combined 
broadside and end-fire array. The following formula 
expresses the ratio of these two effects: 


T 
cos (= sin 0.) 
F (90°, 62) D 
———— = tan@, : (24) 


= ta 
7(0%0,) T 
cos (= cos ) 


It is seen from (24) that when 6, is small, this ratio is 
large and the broadside effect is predominant; while 
when @, is large, the ratio becomes small and the end- 
fire effect becomes significant. When 0,=45°, the ratio 
is unity and the two effects are equal. 

The H plane radiation pattern is found to be: 


Tv . 
cos E sin 6, Cos 6 


G(¢, 94) = COS 64, (25) 
(, 60) 1 — sin? 0, cos? 


and the corresponding electric field is 


—f&olo e iBoko 


G(d, 6a). (26) 


Es"(¢, 9.) =e 
0 
From (25), it is found that the pattern is nearly 
independent of ¢ when the tilt angle 0. of the antenna 
is small. In practical application, @. is usually set to be 
small (say 20°) so that the radiation pattern in the 
H plane is nearly omnidirectional. 


RADIATION PATTERNS OF ARRAYS 


Certain antenna arrays which are composed of a 
number of dipole antennas coupled electromagnetically 
to a two-wire line are shown in Figs. 11 and 12. By an 


SHORT - SHORT - 
CIRCUIT CIRCUIT 
5 af 
A 
aN 
45 


DIPOLE 
ANTENNA 


(a) (b) (c) 
Fig. 11—Broadside arrays. (a) Broadside and end-fire array. 


DIPOLE 
ANTENNA 


TWO- WIRE 
LINE 


(a) (b) (c) 
Fig. 12—End-fire arrays. 
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appropriate arrangement of the dipole antennas coupled 
to the line, an antenna array can be made to function as 
1) a combined broadside and end-fire array, 2) a broad- 
side array, or 3) an end-fire array. The radiation pat- 
tern for each array is found as follows: 

1) An array with broadside and end-fire properties 
as shown in Fig. 11(a). The radiation pattern in the E 


plane is 
N 
sin & mw cos 0 
2; 


F(@, 6.) = F (0, 9a ) 
oe bee sin (a cos @) 


(27) 


where F,(6, 0.) is expressed in (22), and N is the total 
number of the dipole antennas. 

2) A broadside array as shown in Fig. 11(b). The 
radiation pattern in the £ plane is 


Tv 
sin (w — cos a) 
2) 
Tt 
sin (= cos s) 
Dae 


+ f2(0, 0a)2 sin (= cos a) 


F(6, 02) = f2(0, 9a) 


; (28) 
cos (m cos 0) 
where 
Tv 
cos (= cos 6a cos ) 
2(0, 0g) = | ———--————_ | cos 9, sin 0 29 
F210, 90) 1 — cos? 0, cos? 6 (29) 
Tv . 
cos (= sin 6, sin ) 
f2(, 02) = | ————————————_ | sin 0, cos 0. (30) 


1 — sin? 6, sin? 6 


3) A broadside array as shown in Fig. 11(c). The E 
plane radiation pattern is 


3m 
| sin (w — cos a) 
4 


F(6, 6a) ad F, (8, Cm) 
ee Ea 
sin | — cos @ 
2 


(31) 


IRE TRANSACTIONS ON ANTENNAS AND PROPAGATION 


4) An end-fire array as shown in Fig. 12(a). The E 
plane radiation pattern is 


Tv 
cos (= sin a) 
2 


cos 0 


sin (w= (fies 0) | 
Cae 


F(6, 0a) = 


(32) 


5) An end-fire array as shown in Fig. 12(b). The © 


FE plane radiation pattern is 


J sin (v= (1 — cos 0) 


7 
cos (= cos s) 
v 


6) An end-fire array as shown in Fig. 12(c). The & 
plane radiation pattern is 


Tv 
—2 cos (= cos ) 
y) 


sin E - (1 — 2 cos 0 | 


F(8, 8a) = 2fz(8, 6a) (33) 


F(6, 9a) =o 


= f2(9, 92) 


cos (7 cos 4) 


sin E = (1 — cos 0 | 


Tv 
cos (= cos a) 
2 


The H plane radiation patterns for these arrays are 
assumed to be nearly omnidirectional. 


+ fe(8, 8a) (34) 


REFERENCES 


[1] C. J. Sletten, “A New Antenna Radiator for VHF-UHF Com- 
Het | AFCRC, Cambridge, Mass., Tech. Rept. No. 57- 
114; 1957. 

[2] S. R. Seshadri and K. Iizuka, “A Dipole Antenna Coupled Elec- 

tromagnetically to a Two-Wire Line,” Cruft Lab., Harvard Uni- 

versity, Cambridge, Mass., SR No. 21; 1958. 

[3] K.-M. Chen, “A Dipole Antenna Coupled Electromagnetically 

to a Two-Wire Transmission Line,” Cruft Lab., Harvard Univer- 

sity, Cambridge, Mass., SR No. 6, Ser. 2; 1960. 

[4] R. W. P. King, “Theory of Linear Antennas,” Harvard Univer- 
sity Press, Cambridge, Mass., p. 72; 1956. 

5] Ibid., p. 168. 

6] Ibid., pp. 168-179 or pp. 232-238. 

[7] K.-M. Chen, “Radiation patterns of antenna arrays consisting 
of dipoles coupled to a two-wire line,” Cruft Lab., Harvard 
University, Cambridge, Mass., SR No. 8, Ser. 2; 1960. 


September — 


ao & nt. Stas 
a TA Se ie x 


bs 


\ 


Saath i th 


pgm 


Je 


Ns 


le ale ah i ie ll A peas hh a bh NTE Ne eee Peete A eels oe ad eh Ae 


A ailing 


as Ae 


IRE TRANSACTIONS ON ANTENNAS AND PROPAGATION 


433 


Mutual-Coupling Effects in Scanning 
Dipole Arrays* 


L. A. KURTZ{, memper, mre, R. S. ELLIOTT, Fettow, ire, S. WEHN}, stupENT MEMBER, IRE, 
AND W. L. FLOCK||, Memper, IRE 


Summary—tThe perturbations in mutual coupling due to finite- 
ness of a dipole array and its ground plane have been studied. 
Experimental procedures were evolved to measure these perturba- 
tions and their influence on antenna performance as a function of 
beam position. Data for one- and two-dimensional dipole arrays are 
presented and interpreted, particularly with reference to edge effects. 


INTRODUCTION 
|e ees scanned antenna arrays 


have been receiving increasingly favorable at- 

tention in recent years because of several claimed 
advantages. These include capabilities to achieve rapid 
scanning and to create a pattern with a very low side- 
lobe level. Some workers have questioned the validity 
of this latter claim, pointing to the feeding difficulties 
caused by the changing influence of mutual coupling on 
aperture distribution as the beam is scanned. This 
paper reports the results of an investigation of the effect 
of mutual coupling in various scanning dipole arrays. 
This investigation was conducted in conjunction with a 
theoretical study on the same topic carried out at the 


~ Massachusetts Institute of Technology Lincoln Lab- 


oratories. Another theoretical study on mutual cou- 
pling effects was undertaken at about the same time at 


~ the Stanford Research Institute.’ 


The term mutual coupling as used in the present 
context needs some elaboration. Imagine a uniformly- 
spaced array of identical radiators, each fitted with a 
pair of terminals. Let the impedance as seen at one 
pair of terminals be measured with all other terminals 
open-circuited. We shall call this impedance the self- 
impedance of the element. Next, let the impedance be 
remeasured with currents flowing at all other terminals. 
The result, called the active impedance, will differ from 


* Received by the PGAP, November 8, 1960; revised manuscript 
received, May 3, 1961. The work reported in this paper was per- 
formed by Rantec Corp, under subcontract to Lincoln Laboratories, 
Mass. Inst. Tech., Lexington, Mass., with the support of the U. S. 
Air Force.— : 

+ Rantec Corp., Calabasas, Calif. 

t University of California at Los Angeles. Consultant to Rantec 
Corp., Calabasas, Calif. 

b tha versity of Alaska, College, Alaska. Formerly consultant_to 
Rantec Corp., Calabasas, Calif. . 

1S, Edelberg and A. A. Oliner, “Mutual coupling effects in large 
antenna arrays: Part I—Slot arrays,” IRE TRANS. ON ANTENNAS 
AND PROPAGATION, vol. AP-8, pp. 286-297; May, 1960. 

2S, Edelberg and A. A. Oliner, “Mutual coupling effects in large 
antenna arrays: Part II—compensation effects,” IRE TRANS. ON 
ANTENNAS AND PropaGATION, vol. AP-8, pp. 360-367 ; July, 1960. 

3P. S. Carter, Jr., “Mutual impedance effects in large beam 
scanning arrays,” IRE TRANS. ON ANTENNAS AND PROPAGATION, 


vol. AP-8, pp. 276-285; May, 1960. 


the self-impedance because the currents flowing in the 
other radiators create electromagnetic fields which in- 
duce voltages in the radiator being measured. It is 
obvious that the phases and amplitudes of these in- 
duced voltages will depend on the phases and ampli- 
tudes of the currents flowing in the other radiators. 
Thus, the active impedance is a function of the aperture 
distribution and will vary with scan. Since all spacings 
in the array, considered as fractions of a wavelength, 
are frequency-sensitive, it follows that the active im- 
pedance also varies with frequency. 

If the array and its ground plane were infinite in ex- 
tent and if the elements were excited with a common 
amplitude and a uniform progressive phase, symmetry 
conditions would require that the active impedance of 
all elements be the same. These active impedances 
would vary with scan angle and frequency, but always 
in synchronism. 

Because the array and its ground plane are finite, con- 
ditions are somewhat different. The finiteness of the 
ground plane alters the image of each element, the al- 
teration being greatest for elements near the edge. 
The finiteness of the array alters the total mutual 
coupling to an element, this effect being greatest for 
elements near the periphery of the array, since they 
are closest to the elements missing from the infinite 
array. Further, the finiteness of the array and its ground 
plane influence the coupling between an element and a 
distant antenna’ used to measure the array pattern. 
Here again, the edge elements are most affected. 

These edge effects combine to cause the active im- 
pedances of the array elements to differ from each 
other and to vary unequally with frequency and scan 
angle. If the edge effects are ignored and the harness 
which feeds the array is provided with network com- 
pensation based on the assumed behavior of an infinite 
array, perturbations in the aperture distribution will 
result, and pattern performance will be affected. Ad- 
mittedly, the larger the array, the smaller the per- 
turbations, because the edge elements then comprise a 
smaller fraction of the total elements. But the larger 
the array, the greater the potential for extremely low 
side-lobe levels. A small perturbation has more effect 
on a 60-db side-lobe level design than it does on a 40-db 
design. Thus, the need arises to assess the seriousness of 
the edge effects and the limitations they put on pattern 
design. It is to the attention of this problem that. the 
present study has been directed. 
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ANALYSIS 


Assume an array of N radiating elements, not neces- 
sarily identical, arbitrarily spaced. Let N cross sections 
be taken in the transmission lines which connect to the 
elements, far enough back from each radiating element 
to assure that, at each cross section, a pure, single 
transmission-line mode exists. In addition to the array, 
assume a distant antenna, also fed by a transmission 
line, at some cross section of which a single transmission- 
line mode can be identified. Then a simple extension of 
a proof given by Silvert® shows that an equivalent 
multi-mesh circuit is a valid representation for the 
array plus distant antenna. The terminals of the circuit 
correspond to the cross sections in the transmission 
lines, and the voltages and currents active in the circuit 
correspond to the voltages and currents of the trans- 
mission-line modes as they exist at the cross sections.® 
Such a formulation avoids the gap problem,’ and since 
a knowledge of conditions in the gap is not needed for 
network compensation, this is a desirable simplification. 

The mesh equations for the equivalent circuit can be 
written 


VigesleZon a diZoy a 8 a Zon 
ee AZ Ae ic LN iN: 


Vy ona IoZno + LhZni + a +P IyZyn. (1) 


where VY, and J, are the voltage and current at the termi- 
nals of the distant antenna; V,, and J, (n=1,2,---,WN) 
are the voltage and current at the terminals of the uth 
element of the array; Z,, is the self-impedance of the 
nth element, and Z,» is the mutual impedance between 
the nth and mth elements. By reciprocity, Zam=Zmn- 

For a transmitting array, Vo= —IpZox, in which Zoz 
is the load placed across the terminals of the distant 
antenna. Then, from (1), 


N 
1G pa ahvAve (2) 


n=1 


As the distant antenna is moved over the surface of a 
sphere whose center lies at the array, the current J 
traces out the transmitting pattern of the array through 
the dependence of the impedances Zo, on angular posi- 
tion of the distant antenna relative to the array. The 
currents J, which appear in (2) can be found by solving 
(1). In any practical case, the terms IpZ,) which ap- 
pear in this solution can be neglected, giving 


4S. Silver, “Microwave Antenna Theory and Design,” Mass. 
Inst. Tech. Radiation Lab. Ser., vol. 12, pp. 53-60; McGraw-Hill 
Book Co., Inc., New York, N. Y.; 1949. 

°R.S. Elliott, “A Critical Review of the Reciprocity Theorem,” 
Rantec Tech. Memo. No. P183-1; September, 1957. 

°P. S. Carter, Sr., “Circuit relations in radiating systems and 
applications to antenna problems,” Proc. IRE, vol. 20, pp. 1004— 
1041; June, 1932. 

TR. W. P. King, “Theory of Linear Antennas,” Harvard Univer- 
sity Press, Cambridge, Mass., pp. 844-848; 1956. 
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By appropriately varying the voltages V;, the currents 
I, in (3) can be made to vary in such a way as to make 
the pattern (2) scan. Eqs. (2) and (3) point up the fact 
that the transformations which connect the pattern cur- 
rent J) to the array currents J,, and these currents in 
turn to the driving voltages V;., depend upon the mutual 
impedances between each array element and the dis- _ 
tant antenna, upon the self-impedances of the array 
elements, and upon their mutual impedances with each 
other. 

As an illustration of what this implies, consider the 
case of a linear array of identical elements, uniformly 
spaced and backed by a ground plane. If the array 
were infinite in extent, then by symmetry, 


Lon = F(a, (4) 


in which ¢ and w would be factors dependent on scan 
angle 0, but independent of m. Further by symmetry, 


Zinn = Lmtp.nt+p; (5) 


in which m, n, and p would be any integers. Thus, the 
manner in which the currents J, differed from each 
other would be totally dependent on the voltages V;,. 
So, too, would the pattern be controlled entirely by the 
voltages V;.8 Indeed, this is the assumption which is 
normally made in designing the aperture distribution of 
a uniform linear array. The finiteness of the array in- 
validates (4) and (5) in three distinct ways: . 

1) ZonAfe”Y for € and W independent of n. The de- 
parture is most pronounced for elements near the edge. 
It is due to the finiteness of the ground plane which 
provides only partial images and to the fact that an ele- 
ment is no longer immersed in a sequence of similar 
elements extending infinitely in both directions. 


Amel 
1< n+ pe 


mn m+p n+p 1 < pe an p = a 


This departure is also most pronounced for elements 
near the edge, and for the same reasons. 
3) Zm+pn+p does not exist for 


m=1,---,N n=1,---,N 
Mm 7p = I Or m+ p> 
and/or n+ <1 or Nea pea 


_ * For example, if V,;=e*® (uniform amplitude, uniform progres- 
sive phase), then solving (3) would yield J,2ein@, and solving (2) 
would yield a conical pattern of zero beamwidth at the scan angle @ 
for which y(@) +9=0. 


a 


« 
D061 


5 
i These three types of deviation will be called jointly 
_ the edge effects. Referring to (2), we see that the first 
__ edge effect imposes a need to compensate the array 
_ curents J, in order to achieve the desired pattern. Total 
K- compensation may not be possible, since ¢ and wy may 
~ depend on scan angle and the index m in a manner which 
_ prevents it. 

Turning to (3), we see that all three edge effects im- 
pose a need to compensate the driving voltages V; via 
_ the array impedances Z,,, (through their appearance or 
“nonappearance in the determinants), and via the im- 
 pedances Zon (through their implicit effect on the cur- 
E rents J,). These conclusions are obviously equally valid 

- for two-dimensional transmitting arrays. 

If the array and its harness are reciprocal, everything 
which has been said is equally applicable to a receiving 
_ array. If the harness is nonreciprocal, one can still draw 

these same conclusions about the array and sub-harness 

consisting of that part of the harness extending back 
from the array elements to the positions of the non- 
reciprocal elements. Thus, the need for compensation 
_ of the edge effects appears to be a real one which cannot 
be avoided. What remains to be determined is its im- 
portance. 
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ANTENNA RANGE 
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The problem of providing a suitable antenna range 
for the experimental determination of the importance 
of the edge effects was a perplexing one. The measure- 
ments were to be made at UHF on arrays of dipoles 
with the practical restriction that the array dimensions 
were never larger than 10 by 10. A conventional ap- 
proach for taking patterns would be to mount the array 
vertically and rotate it. To remove the effects of the 
earth and the mounting structure, a high tower of awk- 
- ward shape would be required. The resulting inaccessi- 
bility of the array elements seemed too great a disad- 
vantage. Thus it was decided to mount the arrays 
horizontally and use a “sky hook” antenna as the other 
end of the range. 

An 18-foot square ground plane, elevated six feet off 
the ground and rotatable about a vertical axis, was 
therefore constructed to accommodate a variety of 
arrays. Forty feet to the west of the ground plane, a 

tower 120 feet tall was erected, and 125 feet to the east 

of the ground plane, a pulley was attached to an im- 

bedded concrete stanchion. A traveling rope was in- 
stalled over this pulley and two others were located at 

the base and top of the tower. Thus, the rope formed a 
- triangle with a sagging hypotenuse and vertices defined 

by the three pulleys. The plane of the rope bisected the 
ground plane. The shortest distance from the center of 
the ground plane to the sagging hypotenuse was 60 feet. 
A battery-operated transmitter could be secured to 

- the traveling rope and equipped with either a dipole or 
a loop. This assured a uniform pattern in the plane of the 
rope for either polarization. A motor drive attached to 
the pulley at the tower base, with a synchro takeoff to 
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a pattern recorder, completed the range. 

The rotating feature of the ground plane permitted 
patterns to be taken through any principal plane of the 
array and also permitted conical cuts. Placing the 
ground plane six feet off the ground allowed the im- 
pedance-measuring bridge and operating personnel to be 
under the ground plane, which proved to be essential. 


MuTuAL COUPLING IN ONE-DIMENSIONAL 
DIPoLE ARRAYS 


A variety of dipole arrays, mounted on the 18-foot 
square ground plane, were studied at the single fre- 
quency 592 Mc and at the single interelement spacing 
of one-half wavelength. The dipole element common to 
all these arrays is sketched in Fig. 1. It is constructed of 
one-half-inch diameter tubing and is fed from a conven- 
tional balun composed of the same tubing. The shorting 
plate of the balun is used to mount the dipole a quarter 
wavelength off the ground plane. 


Fig. 1—Basic dipole element. 


In the first experiment, a one-by-ten array was 
arranged along one centerline of the ground plane with 
one end element a half wavelength from the edge. This 
caused the other end element to be practically in the 
middle of the ground plane and permitted separation of 
the various edge effects. With the elements oriented 
parallel, then 45° echelon, then colinear, a bridge and 
slotted line were used to measure Zan, Znm, and Zno for 
the elements of the array.? The complete data for all 
these experiments is too voluminous to be included in 
this paper but is available elsewhere.'° Significant 
excerpts of the data under discussion are shown in 
Tables I and II. 

Interpretation of this data can be simplified by re- 
turning to the discussion of the edge effects in the 
analysis section. For the one-by-ten arrays described 
by the data in Tables I and II, we see that dipole No. 1 
experiences edge effects due to the proximity of the 
ground-plane edge and due to being an end element in 
the array, whereas dipole No. 10 does not experience 
edge effects due to the proximity of the ground-plane 
edge, but only due to being an end element. Yet for 
parallel and echelon orientation, Z,, for the two end 


9 The reader is reminded that the zero subscripts refer to the dis- 
tant antenna; all other subscripts refer to elements in the array. 

10 “Study of Mutual and Self-Impedance Problems in Large 
Scanning Arrays,” Rantec Corp., Final Rept. P183; September, 
1959. 
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elements is about the same and only slightly higher than 
Znn for a central element. Some variance is seen for the 
colinear case, but a positive interpretation is clouded by 
the variability of tip-to-tip spacing between ends of 
adjacent colinear elements. It was found that as little as 
a sixteenth-inch variation in this spacing could cause a 
considerable change in Z;; and Zio. In any event, even 
for the colinear case, the effect of the ground-plane edge 
on Znn seems to be quite small. This was borne out by 
later experiments on two-dimensional arrays. 

The above remarks apply also for Zn,,41. By the time 
one gets to Zn.ny2, even the colinear data have settled 
down to a steady value, supporting the idea that tip-to- 
tip variability is accountable in the colinear case. 

The effect of the finiteness of the array itself shows 
up in the Z,, and Z,,n41 data for all three orientations, 
but appears negligible for Z,.42 and beyond. It ap- 
pears greatest for the colinear case but in general does 
not cause more than a few per cent variation in the 
impedance values. 

Zo is only slightly affected by the edge at broadside, 
but develops more and more unbalance as the beam is 
scanned further off broadside. The variation is great 
enough to require compensation by the time the beam 
is out to 60°. 

A suitably terminated parasitic dipole was placed 
beyond dipole No. 1, and even though the parasite was 
not in contact with the ground plane, its presence 
caused the impedance data for dipole No. 1 to approach 
closely to that of a central dipole, thus enforcing the idea 
that the finiteness of the array itself is the most im- 
portant edge effect. 

The meaning of these edge effects can be appreciated 
from a discussion of the active impedance. For a one- 
dimensional array, this can be written 

Ts Ts 
Ziutive= lit Li t= Zia tt (6) 
I; I; 

We can see from Table I that for a colinear array, 
Z12/Z1221/2, Z13/Z1152—1/16, etc. Therefore, if 
I,;=I12=Ts, etc., (a good assumption in a small segment 
of a large array), it follows that Zactive ranges roughly 
from 2Z,; down to an almost negligible value for an 
element with neighbors on both sides, as the beam is 
scanned from broadside to endfire. However, for an 
element with a neighbor on only one side, the variation 
is only from 1.5Z;, down to 0.5Z1;. This effect on 
Zactive is therefore much greater than the effect caused 
by the variations in Z,, and Z,m noted in Table I. How- 
ever, until we relate effect on pattern to variations in 
mutual coupling, it is not fair to conclude that the 
variations seen in Table I are small enough to be 
ignored. This is a point which shall be discussed further, 
after the data on two-dimensional dipole arrays has 
been presented. 

At this juncture in the experimental program, the 
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Lincoln computer was supplied with smoothed data 


taken from Table I, plus idealized data for Zn0, for the | 
case of the colinear array. Calculations were made, using — 


(1), of the receiving currents in a ten-element linear 
dipole array in which each element was terminated in a 
matched load. Concomitantly, the receiving currents 
were measured vs scan angle for the one-by-ten colinear 
array. The comparison of computed and measured 
values is shown in Table III. The agreement is quite 
good, close to broadside, becoming poorer-as the beam 
is scanned further out. This is attributed to the fact 
that the variations in Z,) data were not included and 
suggests this correction cannot be ignored in wide-angle 
scanning designs. 

Concerning the amplitudes displayed in Table III, 
one can conclude that the variations are chiefly con- 
fined to the two elements closest to each end, even out 
to 50° scan. Even the phase errors are not too bad in 
the middle six elements. 


MUTUAL COUPLING IN Two-DIMENSIONAL 
DIPoLE ARRAYS 


Five-by-five and five-by-ten dipole arrays of various 
element orientations were tested. Only two of these 
arrays will be discussed, since the conclusions drawn 
from the others are in complete agreement. 

Five linear arrays of ten parallel dipoles each were 
mounted on a small ground plane, such that all fifty 
dipoles were a half wavelength on centers and all edge 
dipoles were a half wavelength from the corresponding 
ground-plane edge. Measurements yielded data of which 
the sample in Table IV is typical. Interestingly, the 
corresponding Z, data, when normalized to the respec- 
tive self-impedances, compare quite closely. This was 
found to be true for all fifteen sets of data and indicates 
that the edge effects alter the mutual impedances in 
much the same ratio that they alter the self-impedances. 

The bridge and “sky hook” range were used to meas- 
ure Zno and J; for all fifty dipoles as a function of scan 
angle. Once again, the case of all elements terminated in 
a matched load was used in determining J;. Samples of 
this data are shown in Tables V and VI. 

Comparing the data in Tables I and IV, we see that 
mutual impedance between parallel dipoles in a°com- 
mon row is essentially unaffected by whether there is a 
single row or several rows. In the 45° direction, the com- 
parison must be inferred by interpolation because the 
spacing is 1/2 greater in the two-dimensional array than 
it was in the corresponding one-dimensional array. 
When the interpolation is made, it is found that the 
mutual-impedance values for the one- and two-dimen- 
sional arrays agree quite closely. 

It is a little unfair to compare rows of colinear di- 
poles in the two cases of Tables I and IV, because such 
a row is only five elements long in the two-dimensional 
case. However, the agreement is fairly good. This same 
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TABLE IV TABLE V 
Znm DATA FOR Five BY TEN DipoLE ARRAY Zno AND Ij, DATA FOR FIVE-BY-TEN DIPOLE ARRAY 

0.019[683° | 0.039/526 | 0.043/584 | 0.039/526 | 0.019/683 1.05[5° | 1.05/10 1.03[6 1.03/15 1.07| 
1.00 [0 0.54 |9 0.074{133 | 0.014/283 | 0.004533 — 1.18/31 | 0.97{41 0.95{26 0.90/44 13:2213T 
0.023[524 | 0.059/443 | 0.071/406 | 0.059/443 | 0.023|524 1.04/5 | 1.04[10 1.03 (7.> 1.02/12 1.04/8 
0.380/95 0.179/115 | 0.031/205 | 0.004/318 | 0.005/606 1.29/16 | 1.12/30 1.07/18 1-12(3t 1.29/19 
0.020/379 | 0.074[283 | 0.155[245 | 0.074/283 | 0.020/379 1505165. “l= de02 tae {0156 1.01/9 1.03|7 
0.133/246 | 0.083/269 | 0.023{380 | 0.008/559 | 0.008/694 £,15{15. | -1.02(27- 1.01/10 1.00/30 1.21/18 
0.035/197 | 0.156/101 | 0.380/95 | 0.156/101 | 0.035|197 — 1401125 ee Oa i 0.99|/7 — 1. O1"S +01 See 
0.068/416 | 0.046/439 | 0.025/530 | 0.017/667 | 0.010/787 1.19/10 | 1.05/28 0.99/12 1.02/28 1.23{15 
0.090/142 | 0.549|17__| 1.00 {0 0.549|17 | 0.090/142 1702 Five O25 ne Vee [ss 0.89[7 1.01/4— 
0.038[587 | 0.030/611 | 0.022/683 | 0.017/793 | 0.010]911 £.19/4> | 1.02/23 - P0565. 1.00/23 1.26/12 
0.035/197 | 0.156/101 | 0.380)95 | 0.156/101 | 0.035197 1.04/0 | 1.04/6 1.00/0 © 1.01/6 | 1.04[4- 
0.024/770 | 0.022/787 | 0.017/848 | 0.014/931 | 0.010{1033 (181s Oe 1.00|0 — 1.02/25 1.23/10 
0.020/379 | 0.074{283 | 0.155/245 | 0.074/283 | 0.020/379 1.00/0 | 1.04/6— 1.04[2 1.00/4 — 1.02/1 
0.016/953 | 0.014/965 | 0.013/1014 | 0.011]1085 | 0.008 {1179 123i = Uc 10s its 1.05150 1.02/22 1.22[6 
0.023[524 | 0.059/443 | 0.071/406 | 0.059/443 | 0.023[524 — B02] ad OO A 1.02/2— 10312 = 1.05|2 
0.011/1135 | 0.010/1146 | 0.010/1180 | 0.008/1243 | 0.007/1334 = ay 

1.16]3 1.05/16 0.98|3 _ 1.02/16, sesh se 
0.019/683 | 0.039/526 | 0.043/584 | 0.039/526 | 0.019/683 — : 

1.03|4_ | 1.07{2 F038 tod 1.03/0 — 1.06|3_ 
0.007|1319 | 0.007|1320 | 0.008|/1347 | 0.007|1412 | 0.006/1491 Eye. ¥ 

1725160 Lie kO [08'a 11270" 1.28|3 
0.016/850 | 0.024[778 | 0.030{756 | 0.024778 | 0.016[850 a orien t ‘ ; 

1.04/6 | 1.051 1.03|2_ 1.03) 15 1.047 
0.005/1499 | 0.006/1493 | 0.006/1531 | 0.005/1586 | 0.004/1654° "¥4 ey 

1.19|2 0.98/19 0.93|4— 0.98/21 LATO 

Upper values: WRT central element. 
Zsnir=2.20| 14° g=0°. 
Uppe lues: Zyo. 
Lower values: WRT corner element Tee wee. fe 
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TABLE VI 


~ Zno AND I; DATA FOR FIvE-ByY-TEN DiIpPoLE ARRAY 


0.93|15° | 1.00/16. 0.98/19 1.0017. 0.93] 15 
2.26|31. )91.87|23. 1.78|34. 1.91|26, 2.26|34, 
0.86|15 | 0.96/15, 0.93|16 0.96|17. 0.88|16 
POOL US h! 1.52 10. 1.39] 13° 1.46/1— 1.94|20 
0.83|12 | 0.95/14 0.96|12. 0.9411. 0.88|14 
Pars eh 1238 9-)7 | > 4227|10 1.37415 Leh eh8 
0.82/12 | 0.96|10 0.99|8 0.96|7__ 0.88|10 
1.64/3. | 1.30/17 121717 1.26/26 1.67|33 
0.88|11 | 1.00[6 1703| 0 Ol}. 0.94|4 

Pe ets d0| 0.) 1.08] 1.25/29 1.59/10 
0.89/10 | 1.01|6 1.02|0_ OE Oe 0.91/38 
1.60/6 | 1.26/28 1.02/0 125/35° 1.65/20 
0.92|3_ | 0.97/4 1.00/9— 1.02/10 0.93/12 
1.58/16 | 1.23/33 0.96{7 | 1.18[39 1.63/29 
0.88/6 | 0.98{13 0.98/15) 1.01/15 0.92/17 
t.55/21 | 2.23(38° 0.96{10° 1.15/44 1.62(33 
0.93/8 1.00/14 1.03/18 1.03/23 0.96/24 
ROS Sa pa ee 1.00]10° foi4i 51 1.67/38 
0.92/20 | 1.03/26 1.12/29 1.03/30 0.96/29 
1.50/25 | 1.04[37 1.22/14 1.03/46 1.56[37. 


6=60°. 
Upper values: Zno. 
Lower values: J;, 


five-by-ten array was also tested with all dipoles 
rotated 90°, and the colinear data agrees even more 
closely with the one-dimensional data. Thus, we can 
conclude that one-dimensional measurements of the 
mutual impedance between dipoles will give a good in- 
dication of the mutual-impedance relationships in a 
two-dimensional array. 

Comparing Table II with Tables V and VI reveals 
little difference between the one- and two-dimensional 
arrays. For a broadside beam, Z, is essentially the same 
for the two cases. With the beam scanned off 60°, the 
average agreement in Z,» values is not quite so good, 
but this effect is insignificant, because neither set of 
Zno values has departed drastically from unity. How- 


ever, the diagreement in J; values 7s significant. We 


cannot compare Table III with Tables V and VI be- 
cause the latter deal with scan in the H plane. How- 
ever, tables similar to V and VI can be constructed from 
I, data taken on the five-by-ten array with all dipoles 
rotated 90°. When this is done,!! the one-dimensional 
I, values do agree with the two-dimensional values out 
to 30° scan. However, beyond 30° scan, the disagree- 
ment between the two cases becomes marked. This is 
evident when Table III is compared with Table VII 
wherein the J, data for 45° scan is listed. Table VII 
shows a depression in J, magnitude extending in all 
directions from the central elements, whereas Table III 
shows a tilted distribution in the direction of scan. The 
difference can be attributed to the additional mutual- 
impedance terms which are present in the two-dimen- 
sional case. 

The slotted line was used to measure the active im- 
pedance vs scan of a central element in the five-by-ten 
dipole array. This was done by recording a sequence of 
partial-reflection coefficients and summing the data 
with proper regard for phase. The procedure is some- 
what equivalent to measuring the terms in (6) one by 
one. The interested reader can find the details of this 
procedure elsewhere.!® The results are shown in Table 
VIII. Scanning in the E plane causes an increase in the 
real component and a decrease in. the imaginary com- 
ponent, whereas scanning in the H plane causes a de- 
crease in the real component along with a decrease in 
the imaginary component of the active impedance. As a 
consequence of this, compensation causes diagonal 
scanning to affect the active impedance relatively little. 
Thus, a 45° echelon array, scanned in either principal 
plane, is attractive from the point of view of active im- 
pedance. Two crossed 45° echelon arrays, used for 
polarization diversity, with one-dimensional scanning, 
are much to be preferred over the conventional com- 
bination of horizontal and vertical dipoles. Experiments 
have confirmed these assertions. 


4 “Study of Mutual and Self-Impedance Problems in Large 
Scanning Arrays,” Rantec Corp., Progress Rept. No. 17; July, 1959. 
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TABLE VII 
I, DATA FOR TEN-BY-FIVE DIPoLEe ARRAY 
0.53/12° 0.50|3 0.54[1 0.54/0 0.56]10 
0.67/17 | 0.60/5 0.65(9 — 0.60|7 0.68|15 
0.75/21. 0.64|6 0.68] 10 OrSotan 0.76/18 
0.81/21 | 0.69[0 0.822 0.66/1 | 0.83/19 
0.87/26 | 0.78|6 1.00/0 One 2m 0.86{24 
0.83/30 | 0.69/16 0.956 0.69/13 0.84/28 
0.72/27 | 0.47/9 OWESICES 0.47/11 0.74[28° 
0.72[17 0.61/14 0.90] 25° 0.60) 13 0.70/16 
0302/23)" e086 A 2015, 0.83/6 — 0.9618 
0.86{57 | 0.75[46 0.79{33° 0.74/44 0.88/67. 
E-plane scan 6=45°. 
TABLE VIII 
NORMALIZED ACTIVE IMPEDANCE VS INTERELEMENT ngs Eh FOR CENTRAL ELEMENT OF A FIVE-BY-TEN DIPOLE ARRAY 
Ga onms 
E-Plane Scan— 
0° CO ihc 60° 90° 120° 150° 180° 
0° 0.80+ 70.77 0.78-+70.64 0.80+ 70.42 1.16+ 70.30 2.95-+70.12 15 —740 Negative 
SOE 0.83+ 70.78 0.81+70.66 0.81+70.47 1.15-+-70.37 2.8 +70.40 28 =+710 Negative 
2) A ; n : i 
‘° 60 0.77+ 70.67 0.80+70.68 0.70+ 70.45 0.96+ 70.36 2.15-+-70..33 10 — 71.00 Negative 
= 90° 0.70+ 70.52 0.67+ 70.50 0.60+ 70.36 0.73+ 0.25 1.49+70.13 4.40—71.30 3—710 
120 0.68-+70.42 0.60+ 70.45 0.46+ 70.30 0.50+ 70.13 1.0 —j70.20 2.30—72.10 Negative 
150 0.43+ 0.38 0.33+ 70.36 0.22-+-70.17 0.25—j0.05 0.66 —j0.43 1.53—72.35 Negative 
180 0.26+ 0.37 0.21+ 70.24 0.12+ 70.03 0.21—70.15 0.79—70.57 2.0 —73.8 Negative 


ear rae ae 


= 


~~ 


1961 


It is interesting to note from Table VIII that for far 
out E-plane scanning, the real part of the active im- 
pedance becomes negative. This agrees with the inde- 
pendent calculations of Blasi.” 


ASSESSEMENT OF EpGE EFFECTS 
IN DIPOLE ARRAYS 


From a study of Tables III, V, VI and VII and the 
parallel discussions, one can reach the conclusion that, 
except for wide-angle E-plane scanning, the errors in 
receiving currents are principally confined to the first 
few dipoles near the edges of the array. Furthermore, 
these errors in receiving currents do not, in the range 
quoted, exceed one hundred per cent, and reach that 
value only in the end elements. Accordingly, one can 
make an estimate of the effect these errors might have 
on the antenna pattern. The procedure adopted was to 
assume a 100 per cent error in the receiving current in 
each end element, a 75 per cent error in the penultimate 
elements, a 50 per cent error in the next elements in, a 
25 per cent error in the next elements in, and no errors 
in all remaining internal elements. One-dimensional 
arrays of various numbers of elements were considered, 
always with half-wavelength spacing. A perfect net- 
work attached to the array, which would convert a 
uniform amplitude distribution in the array to a 
Dolph-Tchebycheff pattern of desired sidelobe level, 
was also assumed. It was further assumed that the eight 
error currents, in the four elements at each end of the 
array, had such phases as to cause the greatest sidelobe 
level distortion after being operated on by the network. 


_ The results of the ensuing calculations are shown in 


Fig. 2. A twelve-element array designed for 20-db side- 
lobes gives only 7-db sidelobes, whereas a 144-element 
array designed for 20-db sidelobes gives 15-db side- 
lobes. However, a 144-element array designed for 40-db 
sidelobes gives only 32-db sidelobes. 

It should be emphasized that this is the most pessi- 
mistic case and would not be expected to occur over 
more than a small segment of the total scan range. It is 
evident that the larger the array, the smaller the effect. 
This is consonant with earlier results on random errors® 
and leads to the important conclusion that all sources 


2. A. Blasi and R. S. Elliott, “Scanning antenna arrays of 
discrete elements,” IRE TRANS. ON ANTENNAS AND PROPAGATION, 
vol. AP-7, pp. 435-436: Octoher, 1959. 

13 R. S, Elliott, “Mechanical and electrical tolerances for two- 
dimensional scanning antenna arrays,” IRE TRANS. ON ANTENNAS 
AND PROPAGATION, vol. AP-8, pp. 114-120; January, 1958. 
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Fig. 2—Worst possible degradation in sidelobe level for a dipole 
array with tapered current errors in the four elements on each 
end. 


of error in antenna arrays become less significant as the 
array size is increased. 

In applications where the edge effect must be taken 
into account, we conclude that the effect of missing 
mutual-impedance terms is the most important, that 
Z,9 Variations are unimportant except at wide scan 
angles, and that unbalances between the values of 
Za nek aNd Zn nk are insignificant. 

Other element types have been studied for edge ef- 
fects in arrays, notably helices and log-periodic ele- 
ments. Though the data amassed was not so extensive 
as for dipoles, the conclusions drawn were the same, 
and there is thus strong evidence to believe that these 
conclusions about edge effects have general applicability. 
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Optimum Feeds for All Three Modes of a 
Monopulse Antenna I: Theory* 


PETER W. HANNAN}, SENIOR MEMBER, IRE 


Summary—In a monopulse antenna for use in a tracking radar, 
the requirements for the sum and difference modes are not the 
same. For the amplitude-comparison type of monopulse having an 
antenna whose main aperture is illuminated by a feed, these inde- 
pendent requirements can be met by a feed designed to fulfill two 
conditions. One is that excitation of the feed aperture in the dif- 
ference modes be effectively about twice as wide as in the sum 
mode. The other is that the shapes of all the feed excitations be free 
of any avoidable irregularities. 

The amount of improvement available when these conditions 
are met depends on the optimization point desired for each mode, as 
well as the design of the reference antenna. In a typical case, the 
gains and slopes in the difference modes should increase by several 
db. Furthermore, the near sidelobes and the spillover radiation in 
the difference modes should decrease by an order of magnitude. A 
small increase of gain in the sum mode is also available, and finally, 
the criticalness of positioning the feed may be appreciably reduced. 


I. INTRODUCTION 


N THE FIELD of tracking radars, the need for 
|t See of increased angular accuracy has led 

from sequential lobing and conical scan to simul- 
taneous lobing or monopulse. As the requirements be- 
come even more demanding, it becomes necessary to in- 
corporate advances in the design of monopulse antennas 
to more fully realize their potentialities. One such ad- 
vance is possible in the amplitude-monopulse type of 
antenna consisting of a reflector or lens illuminated by 
a feed. 

For an ordinary single-mode antenna, the optimum 
design of a simple feed for maximum antenna gain is well 
known [1], [2]. In the case of a monopulse antenna, 
there are usually a sum and two difference modes, and it 
has not been customary to design a feed for simul- 
taneous optimum performance in all three of these 
modes. The particular compromise that is made may 
depend on the system requirements, and the relative 
importance of the various modes. One choice often 
made involves a “crossover level” which maximizes a 
certain product of sum gain and difference slope [3], 
[4]; in this case, neither of these quantities achieves its 
optimum value. In addition, some of the other antenna 
characteristics, such as sidelobe levels, spillover radia- 
tions, and criticalness to misalignment, are consider- 
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ably poorer than they need be with an optimum design 
in every mode. te 

It is the purpose of this paper to discuss in detail the 
limitations involved in the usual compromise design, 
and to describe the degree of improvement available if 
each mode can be optimized independently. In a com- 
panion paper, means to achieve this independent opti- 
mization in all three modes will be presented. 


II. DESCRIPTION OF MONOPULSE ANTENNA 


It is assumed that the reader has some familiarity 
with monopulse [5], [6]; however, it is desirable to 
describe the basic monopulse antenna system which is 
considered in this paper. This is shown in. Fig. 1, to- 
gether with the terminology to be used. The essential 
components are a comparator, a feed, and a focusing 
element. 

There are three regions of the antenna system which 
are of special interest. When signals are provided to the 
feed by the comparator, there exists a field across the 
aperture of the feed which is herein termed the feed 
excitation. This, in turn, creates a field across the aper- 
ture of the focusing element which is called the main 
aperture illumination. The resulting field at a large dis- 
tance from the antenna may be described by its dis- 
tribution as a function of angle; this is the antenna 
pattern. 


A. Modes of Operation 


In describing the operation of a monopulse antenna, 
there are two different concepts which are commonly 
employed. One concept considers the signals which may 
exist within the comparator that would excite different 
portions of the feed aperture and create antenna pat- 
terns consisting of simple lobes displaced from the an- 
tenna axis. This may be called the Jobe concept; it stems 
from sequential-lobing and conical-scanning antennas. 
The lobe concept is helpful in understanding the prop- 
erties of a monopulse antenna, and for the ordinary 
monopulse antenna it is often employed as the principal 
means of analysis. However, it does not provide a 
convenient basis for the optimization discussed in this 
paper. In addition, the lobe signals represent an inter- 
mediate stage in the system, and do not have as direct 
an operational significance as may be desired. 

The concept which is employed in this paper con- 
siders the signals existing in the three channels at the 
input of the comparator. Since these three channels cor- 
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Fig. 1—Monopulse antenna system. 


respond to three distinct modes of operation, this view- 
point may be referred to as the mode concept. These 
modes are called the swm mode, the azimuth difference 
mode, and the elevation difference mode, in accordance 
with the type of action occurring within the comparator 
during reception. 

_ When coupled to the transmitter of a radar system, 
the sum mode provides illumination of the distant tar- 
get; when coupled to a receiver it provides range in- 
formation and a reference signal. The azimuth and ele- 
vation difference modes are coupled to receivers whose 
signals, when combined with the reference sum signal, 
provide azimuth and elevation angle information, re- 
spectively. Although only the sum mode actually exists 
in transmission, it is common practice to consider all 
three modes in transmission for ease of analysis; by 
reciprocity the antenna patterns are the same whether 
obtained in transmission or reception. Shown in Fig. 1 
are a sample set of field distributions during transmis- 
sion for the sum, azimuth, and elevation modes, in each 
of the three regions of interest. 


B. Assumptions and Limitations 


A number of assumptions are made in this paper 
which simplify the calculations and presentation. As 
shown in Fig. 1, the feed and main apertures are rec- 


-tangular, and each field distribution is assumed to be 


separable into the product of a horizontal and a vertical 
distribution. In the sum mode, both the horizontal and 
the vertical distributions are of the even type. In each 
of the difference modes, the distribution in one direction 


is of the odd type while in the other direction it is of the 
even type. The antenna is, of course, assumed to be free 
of any dissipation losses. It is also assumed that both 
the feed aperture and the main aperture are large com- 
pared with a wavelength, and that the feed is much 
smaller than the main aperture. As a result of all these 
assumptions, both the main illumination and the an- 
tenna patterns are relatively easy to calculate, and can 
be presented in a single dimensionless form. Many 
practical antennas do not correspond exactly to the 
assumptions made, and in these cases the numerical 
results obtained herein are applicable only as approxi- 
mations; however, the basic concepts and trends are 
generally valid. 

The analysis presented in this paper is specifically de- 
veloped for the amplitude-comparison type of mono- 
pulse, as typified by the antenna shown in Fig. 1. A 
similar approach could be applied to the phase-com- 
parison type [7]. However, in a phase-comparison an- 
tenna of the kind involving a feed and a focusing ele- 
ment, improvements analogous to those suggested in 
this paper are difficult to implement, and are beyond 
the scope of this presentation. 


III]. NEED FOR DIFFERENT WIDTHS OF FEED EXcI- 
TATION IN SUM AND DIFFERENCE MOpDES 


One of the requirements for optimizing a monopulse 
feed in all three modes is that the width of feed excita- 
tion be different in the sum and difference modes. This 
can be shown in several ways, each furnishing a dif- 
ferent insight into the basic processes involved in ob- 


446 


taining good antenna performance. In each case, the 
sum and one difference mode will be investigated; the 
other difference mode would have the same requirement 
as the one considered. 


A. Focal Field 


The first viewpoint is shown in Fig. 2, which illus- 
trates the region around a feed during reception. When 
a plane wave is incident on the focusing element, it is 
focused to a compact area [8] which allows a relatively 
small feed to capture most of the available power. It can 
be seen intuitively (as well as proven rigorously) that 
for maximum signal, the size of the feed aperture should 
be about equal to the size of the main lobe of the focused 
wave, when the power in the sidelobes of the focused 
wave is intentionally wasted in the interests of a prac- 
tical design. For the sum mode, indicated in Fig. 2(a), 
this situation corresponds to the standard design of a 
simple single-mode antenna for maximum gain. 

In the difference mode, the antenna pattern has two 
main peaks of opposite polarity which are displaced 
equal amounts off the antenna axis. When a wave is 
incident on the antenna from the direction of one peak, 
it is focused to an area off the center of the feed, as in- 
dicated in Fig. 2(b); a similar area on the opposite side 
would correspond to a wave incident from the direction 
of the other peak. To obtain the proper signal from 
each wave, the feed aperture may be divided into two 
halves, connected to the symmetrical arms of a hybrid 
. junction. It can be seen that for maximum signal at 
each peak, each half aperture should be about equal to 
the size of the main lobe of one focused wave, just as 
it is in the case of the sum mode. Thus the over-all feed- 
aperture dimension in the odd plane of the difference 
mode should be about twice as great as that for maxi- 
mum sum gain, and in the even plane it should be about 
equal to that for maximum sum gain. It is also interest- 
ing to note that, upon analysis of the output from the 
complete circuit including the hybrid junction, the 
maximum difference gain should be about 3 db below 
the maximum sum gain. 


B. Main Aperture Illumination 


Another viewpoint considers the illumination of the 
main aperture by the feed during transmission; this is 
illustrated in Fig. 3 for a hypothetical feed having 
“constant” excitation in both planes. In order to obtain 
maximum gain in the sum mode, the feed size should be 
such that the illumination is tapered down at the edge 
of the main aperture by about 10 db, as shown in Fig. 
3(a). This is the well-known [1], [2] optimum taper 
which balances excessive spillover loss against excessive 
waste of the main aperture, to achieve maximum gain 
in a simple, single-mode antenna. In some designs, the 
illumination is tapered further down so that the dis- 
continuity at the edge of the main aperture is reduced: 
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Fig. 2—Optimum feed sizes; focal-field viewpoint. 
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Fig. 3—Optimum feed sizes; illumination viewpoint. 


this lowers the level of the sidelobes in the antenna pat- 
tern, at the expense of a small reduction of gain. 

For the case of the difference mode, the considera- 
tions of maximum gain and low sidelobes lead to a sim- 
ilar conclusion: namely, the illumination of the main 
aperture should be appreciably tapered down at the 
edge. In addition, some of the special problems of the 
difference mode, such as criticalness to feed tilt and 
edge asymmetries, place a premium on low illumination 
at the edge of the main aperture. For simplicity, it may 
be assumed that the difference illumination should be 
tapered down by about the same amount as the sum 
illumination, say 10 db, as shown in Fig. 3(b). How- 
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ever, in the odd plane of the difference mode, the feed 
is only about half as directive as in the sum mode; 
therefore the feed aperture dimension in this odd plane 
should be about twice as great as it is for the sum mode. 
In the even plane, the feed aperture dimensions for the 
two modes should be about equal. This conclusion is, of 
course, the same as that reached previously. 

It is instructive to consider the situation in which the 
feed size is optimized for the sum mode, but cannot be 
increased for the difference mode. Under these circum- 
stances the difference illumination reaches a maximum 
close to the edge of the main aperture, as shown in Fig. 
4(a), and the large discontinuity at this location creates 
high sidelobes in the corresponding antenna pattern. At 
least half of the power in the difference mode goes into 
spillover, so that there is about 3-db loss in the dif- 
ference signal compared with the optimum condition; 
thus, the difference peak gain would be about 6 db below 
the sum gain. The high illumination of the edge of the 
main aperture makes the difference mode sensitive to 
antenna misalignment and edge asymmetries, and the 
large amount of spillover is likely to contribute to 
additional sidelobes. These “spillover sidelobes” may be 
especially offensive in permitting spurious signals to 
enter the difference channel; examples of this are co- 
herent signals such as those caused by reflection from 
the ground and reflection from asymmetrical antenna 
surfaces, or incoherent signals such as thermal radiation 
from the ground. 

If the feed size were optimized for a difference mode 
without decreasing it for the sum mode, the sum il- 
lumination would be too narrow. As indicated in Fig. 
4(b), optimization in one difference mode causes about 
half the main aperture to be wasted in the sum mode; 
thus a reduction of about 3 db in sum gain would result. 
Attempting to optimize the feed size in both difference 
modes would create additional losses. 

While it is true that a feed size might be utilized 
which strikes a compromise between optimum sum- 
mode and optimum difference-mode performance, the 
defects mentioned above would still be present to a 
large degree. Thus, it is evident that the ordinary feed 
design having the same effective size in the sum and 
difference modes imposes a limitation which degrades 
the antenna performance in a number of ways, many of 
which are likely to have a significant effect on system 


performance. 


C. Antenna Pattern Characteristics 


Up to this point, the viewpoints chosen have merely 
illustrated the need for different feed-aperture sizes in 
the sum and difference modes, and have not yielded 
quantitative estimates of antenna performance. The 
final viewpoint to be presented is one whereby numerical 
values are obtained for a number of significant charac- 
teristics of the antenna pattern. These values are plotted 
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(a) Difference mode when feed is optimized for sum mode. 
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(b) Sum mode when feed is optimized for one difference mode. 


Fig. 4—Non-optimum feed sizes; illumination viewpoint. 


in Figs. 5-8 (pp. 448-449) as a function of normalized 
feed aperture size, for the case of “constant” shapes of 
excitation of the feed aperture. As will be discussed in 
Section IV, these are not the best shapes for the feed 
excitation. However, they are relatively simple to calcu- 
late, and provide a typical case for study of the effect of 
feed aperture size on the antenna patterns; in addition, 
the useful concept of effective aperture is ordinarily re- 
lated to a “constant” shape. 

For the exact definitions of the various antenna pat- 
tern characteristics, the reader is referred to Section 
VII. The curves themselves are drawn as accurately as 
possible, because they are likely to be helpful during the 
design of an antenna. Each curve originates from an 
explicit formula; these formulas appear in the report 
[9] from which this paper stems. A few of the formulas 
also appear in [10] and [11]. 

Sum Mode: In Fig. 5, a set of curves are shown for 
several characteristics of the antenna pattern in the 
suin mode. These characteristics include gain ratio 
(efficiency), sidelobe ratios, and beamwidth ratio. A 
curve for spillover power ratio is also presented; since 
the spillover cannot be included in a general calculation 
of sidelobes, some indication of the “spillover sidelobe” 
effect is desirable. Also shown for reference is a dashed 
curve of edge taper voltage ratio for the illumination of 
the main aperture. 

It should be understood that all of these quantities 
are determined on the basis of a one-dimensional geom- 
etry. However as: mentioned in Section II-B, the field 
strength of the antenna radiation is separable into the 
product of a horizontal and a vertical distribution, in 
the usual manner for a large rectangular aperture. As 
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Fig. 5—Sum-pattern characteristics vs feed size. 


a result, the two-dimensional sum-mode value for gain 
ratio may be obtained by multiplying together the ap- 
propriate one-dimensional sum-mode values for each of 
the two planes involved. In the case of spillover, this 
multiplication rule applies instead to the fraction of 
power not spilled over, 7.e., one minus the spillover 
ratio. The other pattern characteristics given in Fig. 5 
are also valid for the two-dimensional case, when re- 
‘stricted to the appropriate cardinal plane. 

It may be seen that the curve of gain ratio reaches a 
maximum at a normalized feed size of 0.68, while the 
curves for the relative levels of both the first and second 
sidelobes reach a minimum at a normalized feed size of 
1,00. If one were to define an optimum feed size for the 
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Fig. 6—Difference-pattern characteristics vs feed size. 


sum mode in terms of the best combination of these two 
properties, then an intermediate size, say about 0.80, 
would be chosen. Assuming that the antenna require- 
ments are the same for both planes of the sum mode, 
the same result would apply to both dimensions of the 
feed aperture. 

Difference Mode: In Fig. 6, all but one of the curves 
show properties for the difference mode which are sim- 
ilar to those properties given in Fig. 5 for the sum mode. 
As before, these are given for a one-dimensional geom- 
etry; in this case, only the odd plane of the difference 
mode is considered. The two-dimensional difference-' 
mode values for gain ratio (difference peak efficiency) 
and one-minus-spillover-ratio are obtainable by mul- 
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Fig. 7—Effects of misalignment on patterns vs feed size. 


tiplying the one-dimensional difference-mode value in 
the odd plane (Fig. 6) by the appropriate one-dimen- 
sional sum-mode value in the even plane (Fig. 5). The 
beamwidth ratio, ‘sidelobe ratios, and edge taper are 
valid for the two-dimensional case when restricted to 
the odd plane of the difference mode. 

The curve of gain ratio reaches a maximum at a 
normalized feed aperture size in the difference plane of 
about 1.37, which is approximately twice the size for 
maximum sum gain. The first and second sidelobe ratios 
are minimum at a normalized feed size of 2.00, so an 
optimum size in the difference plane, considering both 
gain and sidelobe level, might be at about 1.60; these 
feed sizes also are twice those obtained for the sum 
mode. In the other plane, the feed aperture size would 
be optimized on the basis of the sum-mode curves, and 
would be about equal to the size obtained previously for 
the sum mode. 


As mentioned in Section III-B, spillover may be an - 


important contribution to sidelobes in the difference 
mode. It may be seen from the curve of spillover ratio 
that at a normalized feed size of 1.60 the spillover power 
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Fig. 8—Round-trip pattern characteristics vs feed size. 


is less than 1/3 of its value at a feed size of 0.80, and 
has been reduced to an approximate, practical minimum. 
Actually, as will be discussed in Section IV, a more ideal 
shape for the feed excitation would permit a further re- 
duction in spillover. 

There is one additional characteristic of the antenna 
pattern in the difference mode which is shown in Fig. 6; 
this is the difference slope, or the slope of the difference 
pattern on the axis of the antenna. This characteristic 
is of fundamental importance in a monopulse antenna 
because it determines the angle sensitivity, or the ability 
to accurately track the angle of a distance source under 
conditions in which errors may be caused by spurious 
signals of either an incoherent nature (receiver noise, 
antenna noise, etc.) or a coherent nature (plumbing 
asymmetries, antenna asymmetries, ground reflections, 
etc.). The difference slope is presented in Fig. 6 as a 
voltage ratio, normalized with respect to the maximum 
possible difference slope (see Section VII). As before, 
it is given for a one-dimensional geometry; the two- 
dimensional case may be determined for the odd plane 
of the difference mode by multiplying the difference 
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slope ratio by the square root of the sum-mode gain 
ratio in the other plane. It may be observed from the 
curve that the difference slope reaches a maximum for a 
normalized feed aperture size of about 1.23 in the odd 
plane. This is somewhat smaller than the size for maxi- 
mum gain at the difference peak, but is still consider- 
ably larger than that for maximum sum gain. 

Misalignment Factors: In Fig. 7 some quantities are 
shown which are concerned with the degradation of the 
monopulse characteristics caused by certain misalign- 
ments or misadjustments within the antenna, While 
these do not have the basic significance of the charac- 
teristics discussed so far, they are, nevertheless, of 
great interest to the antenna designer. One of these 
quantities is called the feed-tilt ratio, and indicates the 
amount by which the central null of the difference 
pattern is filled in when the feed is tilted away from the 
direction of the center of the main aperture. The dif- 
ference-signal component thus created is in quadrature 
with the normal signal and does not shift the angular 
position of the central minimum of the difference pat- 
terns; however, as has been described elsewhere [12], it 
can combine with a phase difference between the chan- 
nels of the plumbing or IF circuits to cause a boresight 
error. The feed-tilt ratio is presented here as the rate of 
change of spurious difference voltage with feed-tilt angle 
for a small amount of tilt, normalized with respect to 
- the maximum possible value (see Section VII). Under 
the assumptions of Section II-B, this ratio is applicable 
to both a one-dimensional and a two-dimensional geom- 
etry without any reference to the other plane. In the 
two-dimensional case, the tilt angle to be considered is 
the component (see Fig. 1) in the odd plane of the par- 
ticular difference mode of interest. 

It is apparent from the curve of feed-tilt ratio that if 
the feed aperture dimension in the difference plane is 
small, the feed-tilt ratio is large and the feed may have 
to be carefully aligned to obtain good performance. 
However, if the normalized feed size increases to a value 
of 2.00, the feed-tilt ratio is zero and the alignment 
becomes noncritical. This is a consequence of the zero 
in amplitude of illumination at the edge of the main 
aperture; under this condition a small sideways shift of 
the illumination causes no significant change. (The 
second-order shape of the curve at this point corre- 
sponds to the second-order shape of the illumination 
curve, and the latter is an unusual property of the 
“constant” shape of difference excitation assumed at 
the feed aperture. A more customary shape would be 
predominantly first-order and would go through a re- 
versal of sign at this point, but the noncritical property 
would still exist.) In many ordinary monopulse anten- 
nas, it is necessary to provide a feed-tilt adjustment 
because tolerances cannot be held tight enough to rely 
on proper alignment of the feed. If the feed aperture 
were enlarged toward the noncritical point, it might be 
possible to eliminate the feed-tilt adjustment. On the 
other hand, one should remember that certain errors in 
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the comparator may also be compensated by a feed-tilt 
adjustment; if this is eliminated, another type of com- 
pensation may be required. 


The other quantities shown in Fig. 7 describe an ef- 4} 


fect of defocusing the antenna, or of moving the feed 
toward or away from the main aperture. The particular 


effect of interest is the change of phase between the on- — 


axis sum and difference signals. In a pure amplitude- 


comparison monopulse antenna, an incoming wave — 
creates sum and difference signals which are in-phase | 


or out-of-phase with each other at the output of an 
idealized comparator. If the antenna is then defocused, 


this phase relationship is modified, and the monopulse 


system becomes a mixture of amplitude-comparison 


and phase-comparison. As described elsewhere [13]; ‘| 
is desirable to keep this phase change, called the hybrid- — 


monopulse angle, within limits. In Fig. 7, two curves 
are shown which are labeled difference defocus ratio and 
sum defocus ratio, respectively. These ratios yield the 
rate of a phase change with axial feed displacement for 
small departures from the focal point, normalized with 
respect to the maximum possible rate of change in the 
difference mode (see Section VII). The difference be- 
tween the two phase changes is the hybrid-monopulse 
angle. Under the assumptions of Section II-B, this is 
applicable to both a one-dimensional and a two-dimen- 
sional geometry, without any reference to the other 
plane. It may be seen that by choosing certain feed aper- 
ture dimensions in the sum and difference modes, the 
values for the difference and sum defocus ratios can be 
made equal, and therefore the hybrid-monopulse angle 
can be made noncritical to small amounts of defocusing. 
The ratio of feed sizes is, however, appreciably greater 
than two; for instance, a sum size of 0.68 would require 
a difference size of about 2.09. 

The properties shown in Fig. 7 may be considered to- 
gether with all of those of Fig. 6 in determining an 
optimum size for the feed aperture in the difference 
mode. The size chosen will, of course, depend on the 
relative importance of the various properties in any 
particular application. In a typical case, a consideration 
of all the properties might give the same result as that 
obtained with only the first six properties in Fig. 6. 
Thus, the normalized feed size in the odd plane would 
have an optimum value of about 1.60 for the difference 
mode, roughly twice as great as the optimum value for 
the sum mode. 

Round-Trip Performance: Earlier in this paper, 
mention was made of a compromise feed size which 
would provide the best performance for the ordinary 
situation in which the sizes cannot be different in the 
sum and difference modes. It is instructive to examine 
such a design, in order to compare it with the more ideal 
one in which the size for each mode may be optimized 
independently. For the sake of simplicity and correla- 
tion with other investigations of this subject, the 


optimization will be confined to two particular quanti-- 


ties; one is the on-axis round-trip sum voltage (sum 
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voltage squared) and the other is the on-axis round- 
trip slope (sum voltage times difference slope). Fig. 8 
shows curves of these two quantities in normalized 
form, as a function of normalized feed aperture size, for 
a two-dimensional case in which the sizes are the same 
for the sum and both difference modes, and the hypo- 
thetical “constant” excitation is assumed for both 
planes. 

A feed size midway between the peaks of these two 
curves may be chosen as representing the best com- 
promise, and it appears that only a small degradation is 
obtained [3]. However, in the case of the round-trip 
slope curve, one factor is the one-way one-dimensional 
difference slope, and the latter is far from its true 
optimum value. 

The values which the two round-trip quantities would 
attain when the feed size is optimized independently in 
the sum and difference modes are indicated in Fig. 8, by 
means of the black dot and the dashed horizontal line. 
For the round-trip slope a significant increase is thereby 
obtained; the actual amount is approximately 2.5 db. 
There would also be, of course, a marked improvement 
in other properties discussed previously, such as side- 
lobes and spillover in the difference modes. Further- 
more, these latter benefits could be considerably greater 
if the feed size were optimized with all of these quanti- 
ties in mind. 


1V. NEED FOR EFFICIENT SHAPES OF FEED EXCITATION 
IN SUM AND DIFFERENCE MODES 


The preceding section has shown the desirability of 
allowing the feed aperture size to be different for the 


~ sum and difference modes. There is another aspect of 


the feed excitation which has an important effect on 
antenna performance: namely, the shape of the excita- 
tion. This will be illustrated by the same three view- 
points employed previously. 


A. Focal Field 


In Fig. 9, the field of an incoming wave in the focal 
region near the feed is shown. The shape of the major 
lobe of this field is smooth and tapers in a simple man- 
ner to zero on each side. It is natural to expect that for 
high gain, the shape as well as the size of the feed 
excitation should be similar to that of the focused field. 
(As a matter of interest, the “matched filter” criterion 
[14] could probably be employed to prove that the 
greatest signal is received when the shapes are iden- 
tical.) By way of example, two simple hypothetical 
shapes of feed excitation may be considered. In Fig. 
9(a), cosine and sine shapes are shown for the sum and 
difference modes, respectively. Since these closely ap- 
proximate the shape of the major lobe of the focal field, 
the gain should be close to the maximum obtainable by 
complete utilization of this major lobe. In Fig. 9(b), 
double-cosine shapes are shown; here the gain should 
be significantly reduced because of the dissimilarity in 
the region of the “holes” between the double cosines. It 
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(a) An efficient set of shapes. 


(b) An inefficient set of shapes. 


Fig. 9—Comparison of feed-excitation shapes; focal-field viewpoint. 


might also be added that this situation is characterized 
by an inherent reflection of the incoming wave by the 
feed, and the reflected power is capable of causing addi- 
tional deficiencies in the performance of the antenna 


[15]. 


B. Main Aperture Illumination 


Fig. 10 presents the viewpoint which evaluates the 
illumination of the main aperture, and the same ex- 
amples of feed-excitation shapes are considered. With 
the cosine and sine shapes illustrated in Fig. 10(a), the 
radiation pattern of the feed has low sidelobes, and the 
power wasted in spillover beyond the edges of the main 
aperture is minimized. With the double-cosine shapes 
illustrated in Fig. 10(b), the sidelobes of the feed pat- 
tern are high, and it becomes necessary to waste an 
appreciable amount of power in spillover. Thus it is evi- 
dent that the former shapes are more efficient than the 
latter; in addition, the excess spillover existing in the 
latter case may aggravate certain problems, as men- 
tioned in Section III-B. 


C. Antenna Pattern Characteristics 


The final viewpoint, given in Fig. 11, presents quan- 
titative values for a few particular properties of the an- 
tenna pattern. These properties are gain ratio in the 
sum mode, slope ratio in the difference mode, and the 
spillover ratios for both modes, for those particular feed 
sizes which yield maximum gains or slopes, as indicated. 
The values are given for a one-dimensional case; as 
before, two-dimensional results are obtainable by multi- 
plying together the proper one-dimensional values. 
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(b) An inefficient set of shapes. 


Fig. 10.—Comparison of feed-excitation shapes; 
illumination viewpoint. 


Three examples are chosen to illustrate the effect of 
shape of feed excitation; these are the single and double 
cosine and sine types mentioned just previously, plus 
the constant types. It is apparent that the single sine 
and cosine shapes (a) are the best and the double-cosine 
shapes (c) are the worst, while the constant shapes (b) 
are intermediate between these two. The differences 
among the various results are large enough to have a 
significant effect on antenna performance. For example, 
in going from the double shapes (c) to the single shapes 
(a), in only the one dimension, there is an increase of 
1.1 db in sum gain (see also [2 ]) and 1.3 db in difference 
slope, and a reduction by a factor of 0.16 in sum spill- 
over and 0.28 in difference spillover. These spillover 
reductions would be even more striking if a larger feed 
size had been chosen. 


V. AN IDEAL FEED EXCITATION 


In Fig. 12, an “ideal” feed excitation is shown in two 
dimensions, for all three modes. This is not to be re- 
garded as a unique ideal, but rather as a typical one 
representing a kind of composite of the previously-dis- 
cussed concepts. For example, the difference excitations 
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Fig. 11—Antenna-pattern characteristics for various 
shapes of feed excitation. 
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Fig. 12—An ideal feed excitation. 


are about twice as wide in the odd planes as are the ef- 
fective sum excitations, and the shapes of the excita- 
tions are free of any avoidable irregularities. If a feed 
were designed to have all these features, the resulting 
antenna performance should be very close to optimum 
in all three modes of the monopulse antenna. 


VI. Discussion 


It is instructive to estimate the improved perform- 
ance obtainable with a feed such as that shown in Fig. 
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12. Suppose that an ordinary feed had an effective nor- 
malized size of 0.80 in both planes and all three modes, 


_and that with the new feed this size becomes 1.60 in the 
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odd plane of the difference modes. Suppose further that 
the ordinary feed had shapes of excitation which were 


double cosine by constant for the even types and sine by 


-odd constant for the odd. The curves and tables in this 


_ paper may now be employed to estimate the improve- 
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ment, in terms of the complete two-dimensional quan- 


tities. 


In the sum mode the gain should increase by about 1.4 


_db and the spillover ratio should (it is estimated) be re- 
duced from about 0.3 to less than 0.03. These changes 


are a result of the improved excitation shape. In the dif- 
ference modes the peak gains should be increased by an 


- average of about 3.5 db and the slopes by an average of 
about 1.9 db, while the spillover ratio should (it is esti- 
mated) be reduced on the average from about 0.5 to 
‘less than 0.05. These changes are mainly a result of the 
doubled feed size in the odd planes. 


The first sidelobe ratios in the difference modes should 
go from about 10.5 db to about 23.7 db, because of the 
doubled feed sizes. This 13-db improvement may be 


complemented by an improved suppression of the far 


sidelobes of about 10 db in all the modes, as a conse- 


quence of the spillover reduction. Finally, the critical- 


ness of positioning the feed should be reduced to about 
1 in tilt and about # in focusing, because of the doubled 
feed sizes. 

The analysis in this paper has been confined to a 


_ rectangular or square main aperture. Although no calcu- 


lations have been made for a circular aperture, it is be- 
lieved that somewhat more improvement in difference 
gain, slope, and spillover, and somewhat less improve- 
ment in the near difference sidelobes, would be obtained. 
It is concluded, therefore, that for any amplitude- 
monopulse antenna involving a large feed and a focus- 


ing element, a feed excitation such as that shown in 


Fig. 12 should provide significantly better performance 
than can be obtained with an ordinary monopulse feed. 


VII. DEFINITIONS 


Gain Ratio: 

For the sum mode, gain ratio is the on-axis antenna gain relative 
to the maximum possible gain. The latter, reference gain, is that gain 
which would be obtained if the main aperture were uniformly illumi- 
nated with no spillover. (An alternate term for gain ratio which is 
often used for microwave antennas is “efficiency”.) Normally, the 
on-axis sum gain is also the peak gain; for the range of feed sizes 
where this is not true, the peak gain is also shown by means of a 
dashed addition to the curve in Fig. 5. A 

For the difference mode, gain ratio is the antenna gain at the 
peak of the difference pattern relative to the maximum possible 
sum gain. (The above definition of difference gain ratio employs the 
maximum possible sum gain as the reference quantity because this 
allows one to compare difference gain with sum gain directly on the 
graphs. It would also be possible to define difference gain ratio em- 
ploying the maximum possible difference gain [16] as the reference 
quantity; in this case the difference gain ratio would be greater than 
that given on the curve in Fig. 6 by the factor 1/0.609, or 2.15 db.) 


Spillover Ratio: 

For each mode, spillover ratio is the ratio of power radiated by 
the feed into the region outside the edges of the main aperture, to 
total power radiated by the feed. (Spillover is likely to contribute 
to sidelobes in the antenna radiation pattern,) 
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Edge Taper: 

For each mode, edge taper is the ratio of voltage at the edge of 
the main aperture to the maximum voltage on the main aperture, 
when the main aperture is illuminated by the feed. (In the case of the 
difference mode, a small feed may radiate a pattern whose peak is 
outside the edge of the main aperture; in this case the value of the 
edge voltage relative to the peak radiated voltage is indicated by the 
dotted portion of the curve in Fig. 6.) 


Sidelobe Ratio: 


For each mode, the sidelobe ratio is the ratio of voltage at the 
peak of the designated sidelobe of the antenna radiation pattern to 
the peak voltage of the main lobe of the same antenna pattern. The 
sidelobe ratios are calculated directly from the main aperture il- 
lumination without including any effect of spillover. 


Beamwidth Ratio: 

For the sum mode, the beamwidth ratio is the angle between 
half-power points of the antenna pattern with a specified feed size, 
relative to the angle with a very small feed (7.e., uniform illumination 
of the main aperture). In the latter, reference case, the angle is 
0.887 \/A radians. 

For the difference mode, the beamwidth ratio is the angle be- 
tween the outer half-power points of the antenna pattern with a 
specified feed size, relative to the angle with a very small feed 
(i.e., odd, linear illumination of the main aperture). In the latter, 
reference case, the angle is 2.03 \/A radians. 


Difference Slope: 

Slope of the center of the antenna pattern in the difference mode, 
assuming no asymmetries. In other words, rate of change of the 
difference voltage with antenna angle, where the angle is small and 
the difference voltage contains no quadrature components. 


Relative Difference Slope: 
Difference slope relative to maximum possible sum voltage. 


Difference Slope Ratio: 

Ratio of the difference slope to the maximum possible difference 
slope. The latter could be obtained if the main aperture were illu- 
minated by an odd, linear voltage pattern [17], [18] with no spillover; 
in this reference case the relative difference slope (see previous 
definition) would be +A/+/3 \ voltage ratio per radian. 


Feed-Tilt Ratio: 

The rate of change of the on-axis difference signal with feed-tilt 
angle, where the amount of feed tilt is small, relative to the rate of 
change which would be obtained with any very small feed (z.e., odd, 
linear illumination of the main aperture). In the latter, reference 
case, the rate of change of difference signal with feed-tilt angle, if 
taken relative to the difference slope, would be 6\F/7A?. 


Difference Defocus Ratio: 

The rate of change of phase of the near-axis difference signal 
with axial displacement of the feed in wavelengths, where the phase 
change is caused only by the phase curvature across the main aper- 
ture resulting from a small amount of displacement of the feed from 
the focal point, relative to the rate of change which would be ob- 
tained with a very small feed (i.e., odd, linear illumination of the 
main aperture). In the latter, reference case, the rate of change of 
difference phase with axial feed displacement in wavelengths is 
37rA2/20F* radians. 


Sum Defocus Ratio: 

The rate of change of phase of the on-axis sum signal with axial 
displacement of the feed in wavelengths, where the phase change is 
caused only by the phase curvature across the main aperture resulting 
from a small amount of displacement of the feed from the focal point, 
relative to the rate of change which would be obtained with a very 
small feed operating in the difference mode (7.e., odd, linear illumina- 
tion of the main aperture). (The above definition employs the small- 
feed rate of change in the difference mode as the reference quantity 
because this simplifies the calculation of the hybrid monopulse 
angle. It would also be possible to define the sum defocus ratio em- 
ploying the small-feed rate of change in the sum mode as the refer- 
ence quantity; in this case the sum defocus ratio would be greater 
than that given on the graph by the factor 9/5, as indicated. 


Round-Trip Sum Voltage Ratio: 

The round-trip sum voltage ratio is the square of the on-axis 
sum-pattern voltage, relative to the maximum possible value. The 
latter, reference value, would be obtained if the main aperture were 
uniformly illuminated with no spillover. For the two-dimensional 
square case given in Fig. 8, the round-trip sum voltage ratio is equal 
to the square of the one-dimensional sum gain ratio given in Fig. 5. 


Round-Trip Slope Ratio: 
The round-trip slope ratio is the product of on-axis sum-pattern 
yoltage and on-axis difference slope, relative to the maximum possible 
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value. The latter, reference value, could be obtained if the main 
aperture had uniform illumination with no spillover in the sum mode 
and in the even plane of the difference mode, and had odd, linear 
illumination with no spillover in the odd plane of the difference 
mode. For the two-dimensional square case given in Fig. 8, the 
round-trip slope ratio is equal to the 3/2 power of the one-dimen- 
sional sum gain ratio given in Fig. 5 times the difference slope ratio 
given in Fig. 8. 
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Optimum Feeds for All Three Modes of a 
Monopulse Antenna II: Practice* 


PETER W. HANNANY, SENIOR MEMBER, IRE 


Summary—In an amplitude-monopulse antenna whose main 
aperture is illuminated by a feed, it is desired to control the feed 
excitation independently in the sum, azimuth difference, and ele- 
vation difference modes. It is possible to describe a feed system, 
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comprising an infinite array of radiators and hybrid junctions, which 


demonstrates the principle of complete independent control of the 
three modes. 


There are several practical feeds which approach a hypothetical 


ideal one in varying degrees. One representative type is the ‘‘twelve- 
horn feed;” although it has some disadvantages, it is versatile and 
quite useful. However the most attractive type for many applications 
appears to be one having a combination of multiple-horn excitation 


in one plane and multimode-waveguide excitation in the other. With 
a simple form of the ‘‘multihorn-multimode” feed, the ideal is sub- 
stantially realized, and a major improvement in monopulse per- 


formance can be obtained, 
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I. INTRODUCTION 


N a previous paper [1] a study was made of the im- 
proved performance available in an amplitude- 
monopulse radar antenna comprising a feed and a 


focusing element. This improvement arises from the 


possibility of independently optimizing the feed in each 


of the three modes of operation, rather than compro- 


mising between them. There was presented a hypotheti- 
cal “ideal” feed excitation which would essentially 
optimize all three modes, but no information was given 


regarding the attainment of such a feed. It is the pur- 


pose of this paper to discuss the means by which this 
theoretical feed may be reduced to practice. 

Some early monopulse antennas have already ob- 
tained improved performance by reducing the amount 
of compromise between the modes. This has often oc- 


curred as a fortunate by-product of a feed design in- 
_tended mainly to achieve some other benefit, such as 
shielding from nearby objects, a simple comparator 
structure, or a short focal length. As a result, the 


optimizing process has not usually been complete. An- 
tennas known to the writer which belong in this cate- 
gory have been developed by Bell Telephone Labora- 


tories [2], North American Aviation Company, Hughes 


Aircraft Corporation [3], and Radio Corporation of 
America [4]. 

About seven years ago, Wheeler Laboratories made 
an investigation of this problem, in connection with the 
design of a monopulse tracking antenna [5], [6], for 
Bell Telephone Laboratories. A principle was con- 
ceived and a number of design techniques were worked 
out for feeds which would give optimum performance 
in all three modes. The principle was subsequently 
applied in a study by Wheeler Laboratories, in which 
the basic merits of the amplitude and phase types of 
comparison were analyzed and compared [7]. 

An investigation of feeds for optimum monopulse 
performance has occupied the attention of several 
other groups. In the first published material known to 
the writer [8], the principle was recognized and ex- 
periments made with a particular technique which ap- 
proaches optimum performance in more than one mode. 
A number of interesting techniques have been studied 
at Naval Research Laboratory and General Electric 
Company. A feed [9] for use-in an operational antenna 
has been developed by Lincoln Laboratory which ap- 
plies the principle to all three modes of operation. 

Several additional monopulse tracking antennas have 
been designed by Wheeler Laboratories for Bell Tele- 
phone Laboratories. Some of these antennas [10] incor- 
porate feeds which have achieved substantially optimum 
performance in all three modes of operation. While it is 
not proposed to describe the details of their design and 
measured performance in this paper, it is possible to pre- 
sent the basic methods which are involved. Various 
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techniques will be discussed, and a comparison of their 
performance will be calculated. 


II. PRINCIPLE OF INDEPENDENT CONTROL OF EXxCcITA- 
TION IN THE SUM AND DIFFERENCE MODES 


It is interesting to inquire whether an excitation of 
the ideal type is theoretically realizable. Ordinary 
methods of exciting a feed yield distributions which 
are limited to having the same width in all modes, and 
having shapes which are good for one mode but not for 
another. Upon consideration of the problem, it becomes 
evident that independent control of the excitation in the 
sum and difference modes is required to fulfill the ob- 
jectives. 

Fig. 1 presents a scheme which, in principle, is ca- 
pable of achieving any degree of independent control 
which may be desired. As shown, the feed aperture is 
divided into many small sections, and each section is 
excited independently in the sum, azimuth difference, 
and elevation difference modes. The amplitude of each 
excitation is determined by the amount of coupling ex- 
isting in the appropriate coupling element of the appro- 
priate distribution network, while the necessary even 
or odd character of the excitation is provided by the 
series of hybrid junctions comprising the comparison 
network. 

It should be mentioned that any undesired irregu- 
larities or ripples in the excitation caused by a nonuni- 
form field of an individual radiator can be effectively 
smoothed by making the spacings small compared to a 
wavelength. Of course the problems of impedance of 
the radiators, as well as coupling between radiators, 
would be very considerable. However, in principle, these 
problems could be solved. 

Since the above scheme enables an ideal feed excita- 
tion to be approximated as closely as is desired, it may 
be stated that such an excitation is theoretically realiza- 
ble. 

Actually, the degree of independent control provided 
by the scheme of Fig. 1 is more extensive than that re- 
quired to achieve the “ideal” excitation previously [1] 
indicated. In that case, the excitation was shown as be- 
ing the product of the horizontal and vertical distribu- 
tions, and the sum distribution was assumed to be 
utilized for the even plane of the difference mode. The 
scheme of Fig. 1, while capable of providing both these 
features, is not limited to doing so. 


III. PracticAL METHODS FOR ACHIEVING 
INDEPENDENT CONTROL OF EXCITATION 


The scheme shown in Fig. 1 is theoretically able to 
provide independent control of excitation in the three 


1 This scheme demonstrates independent control not only for a 
feed in a focusing type of antenna, but also for an array type of 
antenna [7]; in the latter case independent control of the monopulse 
modes is also desirable. An alternate scheme for independent control 
of a practical array has been suggested by Mattingly [11]. 
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Fig. 1—A scheme for demonstrating complete independent 
control of the sum and both difference modes. 


modes of operation. However, the large number of 
closely-spaced apertures and the multitude of circuits 
imply a low probability for achieving a satisfactory 
working design with a reasonable amount of effort. It is 
natural, therefore, to investigate some simpler schemes 
which may be more practical to develop, even though 
their feed excitations are not completely ideal. 

There are two techniques which are considered here: 
one is termed the multihorn method and the other is 
called the multimode method. As an introduction, these 
techniques are first presented for the case in which inde- 
pendent control is the objective in the sum mode and 
only one difference mode, 1.e., control in only one plane. 
For simplicity, it will be assumed that the azimuth and 
elevation planes correspond to the H and E planes, re- 
spectively. 


A. Single-Plane Multthorn Feeds 


In Fig. 2, the multihorn method is illustrated, for 
either H-plane or H-plane control of 8 horns. Also 
shown is a network of hybrid junctions which would 
properly excite these horns. It is evident that in the 
plane of control, the desired two-to-one width of ex- 
citation is available. The shapes, however, are not 
ideal, particularly in the H plane. (This defect is also 
present in an ordinary four-horn cluster having no in- 
dependent control.) Naturally, the two-to-one width 
cannot be obtained in the plane that lacks independent 
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(a) Eight-horn feed for control in £ plane. 


E-PLANE 


(b) Eight-horn feed for control in # plane. 
Fig. 2—Single-plane multihorn feeds. 


control, so that this feed would be far from optimum in 
a system requiring all three monopulse modes. 


B. Single-Plane Multimode Feeds 


Fig. 3 illustrates the multimode method, for either 
E-plane or H-plane control of two horns by means of 
three waveguide modes. (The terms “multimode 
method” and “waveguide modes” refer to the natural 
modes of propagation in a waveguide or horn; this 
should not be confused with the three so-called modes 
of operation, which are the sum, azimuth difference, 
and elevation difference modes.) Also shown in the fig- 
ure are the two waveguide-mode combinations which in 
this case were chosen for the even excitation. As before, 
the desired two-to-one width of excitation is available 
in the plane of control. In addition, the shape of the 
sum mode in the plane of control is rather good. Some 
of the other shapes, however, are not ideal, especially 
those existing in the case of E-plane control. Again, of 
course, the two-to-one width cannot be obtained in the 
plane that lacks independent control, so that this feed 
would also be far from optimum in a system requiring 
all three modes of monopulse operation. 


4" 


ae 


Be ele 


0 2 ove 
I>-p 
: 
E 
aN 
uP 


(a) Triple-mode feed for control in £ plane. 


E-PLANE 


(b) Triple-mode feed for control in H plane. 
Fig. 3—Single-plane multimode feeds. 


C. Dual-Plane Multihorn Feed 


Now consider the application of these methods to- 
wards obtaining independent control in the sum mode 
and both difference modes. Fig. 4(a) shows a multihorn 
method for the case of twelve horns, together with an 
appropriate network of hybrid junctions. In this case 
the desired two-to-one widths of excitation are available 
in both planes, but again, the shapes of excitation fall 
short of the ideal. It is interesting to note that this 
multihorn configuration is basically similar to the 
scheme of Fig. 1; however, the limited number of aper- 
tures and lack of amplitude control prevent complete 
optimization from being realized. On the other hand, a 
very substantial improvement in performance would be 
gained in comparison to the ordinary four-horn cluster; 
this will be illustrated quantitatively in a later section 
of this paper. 

It should be realized that there are a number of other 
multihorn configurations available. For instance, the 
partitions in the outer eight horns of the twelve-horn 
feed may be removed to yield an eight-horn feed. This 
would give somewhat better performance in the eleva- 
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E-PLANE 


(b) Four-horn triple-mode feed. 


Fig. 4—Dual-plane feeds. 


tion difference mode and would eliminate four hybrid 
junctions. However, these advantages would have to be 
weighed against the more complicated excitation shapes 
that would exist in such a feed. Another configuration 
would involve the addition of partitions or more horns 
to the twelve-horn feed. This would permit a closer ap- 
proach to the scheme of Fig. 1, thereby allowing more 
complete optimization; however, the physical structure 
would become more complex. 


D. Dual-Plane Multimode Feed 


Having combined the single-plane multihorn methods 
of Fig. 2 into the dual-plane multihorn configuration of 
Fig. 4(a), it is natural to ask whether the multimode 
methods of Fig. 3 can be similarly combined. On in- 
vestigation of this question, one concludes that, while 
there is no apparent theoretical limitation on such a 
process, it is quite unattractive because of the practical 
difficulties involved in properly controlling the many 
waveguide modes required. In addition, the elevation 
difference mode would still have a rather poor shape, 
unless even more waveguide modes were employed. 
However, it is of interest to note that a hypothetical 
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scheme employing an unlimited number of waveguide 
modes would provide an approach to complete inde- 
pendent control leading to the same theoretical conclu- 
sions as did the scheme of Fig. 1. 


E. Dual-Plane Multihorn- Multimode Feed 


There are still further combinations available among 
the methods presented in Figs. 2 and 3; these combina- 
tions involve a multikorn method in one plane and a 
multimode method in the other. Fig. 4(b) presents such 
a case, in which the better plane of each method has 
been retained. The result is a feed excitation in which 
not only are the desired two-to-one widths available in 
both planes, but the shapes are effectively ideal in the 
H plane and depart only a moderate amount from the 
ideal in the £ plane. 

The feed as shown comprises four identical horns, 
each utilizing three waveguide modes. As might be ex- 
pected, there are a number of alternate multihorn- 
multimode configurations which might be considered. 
For instance, the outer two horns might be reduced in 
width and excited with only the 1-mode; this would re- 
duce the physical size of the feed at the expense of a 
somewhat more complicated set of excitation shapes. 
Another example involves the use of more than the four 
horns shown, thereby permitting greater control of the 
excitation in the E plane. This would allow a closer 
approximation to the ideal shapes, at the expense of a 
more complex excitation network. 

Returning to the simple configuration shown in Fig. 
4(b), it is interesting to mention the factors relating to 
its usefulness in a real antenna system. As regards the 
practicability of this scheme, it may be noticed that the 
configuration of the hybrid network is less complex than 
that of Fig. 4(a), and is confined essentially to a single 
plane. The generation of the proper waveguide modes is 
certainly a problem of design; however, these modes 
lend themselves to techniques which are quite feasible, 
and can yield results very close to the ideal. Thus the 
“multihorn-multimode” method appears to be most 
attractive as a practical scheme for realizing almost 
completely the benefits obtainable with independent 
control of the sum and both difference modes. . 


IV. COMPARISON OF PERFORMANCE FOR SEVERAL 
PRACTICAL FEEDS 


It is instructive to compare a few practical feeds on a 
quantitative basis. Fig. 5 presents such a comparison, 
based on calculated values for three properties: sum 
gain ratio, difference slope ratio, and spillover ratio. 
For simplicity, the size and the mode ratio of each feed 
is chosen to yield maximum sum gain; as discussed in 
the previous paper [1], this is a somewhat smaller size 
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FEED EXCITATION 
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FOR MAXIMUM SUM GAIN 


(c) TWO-HORN TRIPLE- MODE 


(d) TWELVE-HORN 


(e) FOUR-HORN TRIPLE-MODE 


Fig. 5—Antenna-pattern characteristics for 
various practical feeds. 


than that which might be chosen as an optimum when 
several properties are considered. As a result, the listed 
gains and spillovers are somewhat greater than they 
might be in a more typical case. Furthermore, the 
slopes turn out to be appreciably different from those 
which would be obtained with a larger feed. The as- 
sumptions made regarding a rectangular aperture, a 
large feed in wavelengths, etc.,?2 are the same as were 
made in the previous paper, and all the definitions and 
symbols are the same. The values listed are the com- 
plete two-dimensional quantities, as obtained by multi- 
plying together the proper one-dimensional values. As 
before, it is assumed that the azimuth and elevation 
planes correspond to the H and E planes, respectively.’ 


* Since practical antenna defects (aperture blocking, reflections, 
cross-polarization, “space-attenuation,” dissipation, tolerances, etc.) 
are neglected, the listed gains and slopes are likely to be higher than 
the values encountered in actual practice. On the other hand, the 
assumed large feed in wavelengths may have somewhat more spill- 
over loss, and hence less gain, than one about a wavelength in size, 
particularly in those cases involving poor shapes of excitation. 

> In the case of feeds giving lower spillover in the H plane than 


in the £ plane, the harmful effects of the ground might be minimized 
by reversing this assumption. 
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A. Four-Horn Feed 


The ordinary four-horn cluster is shown in Fig. 5(a). 
This has, of course, the defects associated with an ex- 
cessively narrow difference excitation and the poor 
double-cosine shape for the H-plane even excitation. 
The main virtue of this feed appears to be its conceptual 
simplicity [5], [6], including its susceptibility to the 
historical lobing method of analysis. It also has a rela- 
tively small aperture size. 


B. Two-Horn Dual-Mode Feed 


In Fig. 5(b), another common feed is shown; in this 
type the E-plane partition, present in the four-horn 
cluster, has been removed. This allows two waveguide 
modes to propagate, and the H plane is characterized 
by a single cosine for the even excitation and a sine for 
the odd excitation. The single cosine is not only an 
efficient shape, but its effective width is narrower than 
that of a double cosine of the same over-all dimensions. 
As a consequence, the over-all feed-aperture size in the 
H plane can be increased, thereby allowing the azimuth- 
mode excitation to widen toward the desired value. 
However this width falls short of the ideal two-to-one 
difference-to-sum condition, and no further improve- 
ment can be made with only the two modes which are 
available. In addition, of course, the elevation-mode 
excitation remains too narrow. This type of feed [12] 
has found most use [3] in applications which place a 
premium on a simple mechanical structure and rela- 
tively short focal length. 


C. Two-Horn Triple-Mode Feed 


In Fig. 5(c) the feed shown in Fig. 3(b) is analyzed. 
Here the H-plane dimension of the feed has been set to 
maximize the difference slope, and then the mode ratio 
of the even excitation (Z;/E:) has been set so that the 
sum gain is also maximized. While the performance in 
the sum and azimuth difference modes is excellent, 
that in the elevation difference mode, of course, re- 
mains poor. This type of feed is most helpful in applica- 
tions [4] which require a very short focal length; it is 
believed that RCA was the first to incorporate the triple- 
mode technique in a monopulse antenna. 

It is evident that the feeds of both Part B and Part G 
suffer from lack of independent control in the elevation 
plane. It should be noted that in certain antenna appli- 
cations there may also be a problem caused by the dif- 
ferent antenna patterns which would occur in the azi- 
muth and elevation planes. These defects are virtually 
eliminated by the incorporation of independent control 
in both planes. The two practical feeds which achieve 
this are to be shown in Part D and Part E. 


D. Twelve-Horn Feed 
In Fig. 5(d) the feed shown in Fig. 4(a) is analyzed. 
Here the dimensions of every horn are chosen to be 
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.identical, for simplicity and ease of calculation. This 


happens to yield the two-to-one width condition, so 
with the size set to maximize sum gain, the difference 
slope is fairly close to the maximum value obtainable 
with this feed. As expected, the performance is good in 
all three modes, but is still somewhat less than optimum 
because of the inefficient shape of the feed excitation in 
the H plane. In many applications, this defect, together 
with the complexity of the circuitry, would tend to limit 
the attractiveness of the twelve-horn feed. On the 
other hand, it may be the only possible solution in some 
cases. One such case might be a diversity-polarization 
antenna using dual-polarized horns; perhaps another 
case is a low-frequency antenna in which a dipole feed 
array must replace the horn array. There are other, 
special uses for independent control in monopulse an- 
tennas where the twelve-horn method may be the most 
practical solution.‘ It is believed that the twelve-horn 
method was first employed by Lincoln Laboratory in a 
feed [9] developed for a versatile monopulse antenna. 


E.. Four-Horn Triple-Mode Feed 


In Fig. 5(e) the feed shown in Fig. 4(b) is analyzed. 
Here the H-plane dimensions and mode ratios are set as 
they are in Part C, and the E-plane dimensions as in 
Part D. It is evident that this feed yields the best per- 
formance in all three modes when compared with any 
of the others. Furthermore, the performance is rather 
close to that obtainable with the “ideal” feed [1]. 

The numerical results given in Fig. 5 illustrate the ad- 
vantages to be gained with a practical form of inde- 
pendent control. In going from the ordinary four-horn 
cluster to the four-horn triple-mode feed, for instance, 
the sum gain has increased by 1.1 db and the difference 
slopes have increased by 3.9 db, while the total sum 
spillover has decreased to about 3 and the total differ- 
ence spillovers to about } of their original values. The 
reduction of difference spillover is particularly signifi- 
cant, since originally it was so great. Although not 
listed in Fig. 5, various other properties of the difference 
modes, such as gain, sidelobe level, and criticalness to 
misalignment are, as discussed in the previous paper 
[1], also greatly improved. 


4 For example, certain antenna systems obtain sequential lobing 
by combining the sum and difference signals of a monopulse antenna. 
Where independent control is available, the lobing performance may 
be greatly improved. With a twelve-horn feed and comparator, this 
feature is relatively straightforward to provide, as follows. Those 
ports of the comparator which would furnish even excitation of the 
outer horns in Fig. 4(a) are not terminated but are properly coupled 
to the difference ports, and these latter ports are combined alter- 
nately in-phase and out-of-phase with the sum ports in the proper 
manner. The result is a sequential lobing system having high cross- 
over level and gain together with low sidelobes and spillover, com- 
pared with a lobing system using a four-horn feed. Lincoln Labora- 
tory has worked out [9] some alternate schemes which instead pro- 
vide independent channels for sequential lobing and communication. 
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V. DISCUSSION 


The analyses in these papers have assumed that the 
feed aperture is large compared with a wavelength. In 
many antennas this simplifying condition does not 
exist, and it is appropriate to ask whether the conclu- 
sions obtained herein would be significantly altered. It 
happens that when the feed aperture becomes small 
compared with a wavelength, one particular theoretical 
approach indicates that the main illumination would 
appear to resemble the optimum, regardless of the 
widths and shapes of the feed excitation. In such a case, 
therefore, one might be led to expect that a “simple” 
feed, not incorporating independent control, would yield 
good monopulse performance. While it is possible that 
some of the benefits of independent control could be ob- 
tained in this manner, the electrically small size of such 
a feed would introduce difficulties and defects which 
would considerably degrade the antenna performance. 
In the ordinary single-beam antenna, these problems 
are well-known; they include narrow bandwidth of 
impedance match, high-power limitations, and the like- 
lihood of excessive cross polarization. For a monopulse 
antenna, these difficulties are usually aggravated, and 
additional feed problems, such as defocusing errors and 
coupling between radiators, are encountered. Thus it is 
believed that an electrically small feed should be 
avoided. 

Another possibility exists for achieving the benefits 
of independent control without the additional circuits 
or waveguide modes described in this paper. Certain 
devices, when placed just in front of an ordinary mono- 
pulse feed, have the property of affecting the radiation 
of the feed differently for the various modes of opera- 
tion. For example, a thin metal fin in the central H 
plane of a feed affects only the difference mode having 
odd symmetry about the fin. Such a fin was employed 
about seven years ago in a monopulse antenna [5], [6] 
designed by Wheeler Laboratories; it narrowed an ex- 
cessively wide radiation from a four-horn feed in one 
difference mode, as well as improving the impedance 
match in this mode. However, it was found that the fin 
could not be designed to provide an effect equivalent to 
the ideal two-to-one feed size without introducing severe 
defects in various other properties of the nonopulse an- 
tenna. More recently, a similar limitation on the bene- 
fits of a fin was found in the case of a two-horn triple- 
mode feed. Although there may be some exceptions, it is 
probably fair to consider these experiences typical. 
Therefore it is believed that “parasitic” feed devices, 
while quite useful for certain limited purposes, are not a 
reliable means for designing an optimum monopulse 
feed. 

Returning to the question of feed size, there remains 
to be considered the more common case, in which the 
feed aperture is neither large nor small compared with a 
wavelength. Under these conditions, the ordinary feed 
design may incorporate to some degree those desirable 
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radiation properties obtainable with independent con- 
trol. Nevertheless, when an independently controlled | 
feed is designed, an improvement approaching that in- 
dicated in previous portions of this paper is likely to be 
obtained. About five years ago, for example, measure- 
ments were performed at Wheeler Laboratories with a 
monopulse antenna which initially utilized an ordinary 
four-horn cluster as the feed. The over-all aperture di- 
mensions of this feed were about 3/2, which was elec- 
trically just large enough to yield acceptable perform- 
ance in a particular application, and the edge taper of 
the main illumination was about 11 db in the sum mode 
for maximum gain with a circular main aperture. Upon 
substitution of a feed having the same effective aperture 
size but a more efficient excitation shape, the antenna 
gain in the sum mode increased by about 1.0 db. When 
a further substitution was made in which one dimension 
of the feed aperture was increased to approach the de- 
sired two-to-one condition, the gain in the correspond- 
ing difference mode increased by more than 3.0 db and 
the sidelobe level changed from 12 db to 18 db below 
the difference peak. Thus the expectation that a wave- 
length-size feed aperture would encounter almost the 
full benefit from independent control has been experi- 
mentally confirmed for a particular case. 

In some more recent antenna designs by Wheeler 
Laboratories, the multihorn-multimode technique has 
been applied to obtain simultaneous optimum perform- 
ance in the sum and both difference modes. One case 
[10] involves a feed aperture very large compared with 
a wavelength, while another has over-all aperture di- 
mensions of about 2 wavelengths. Both cases have 
yielded the benefits expected from independent control, 
and have proven satisfactory from a practical stand- 
point. It is concluded that the incorporation of a prac- 
tical form of independent control in the feed of a 
monopulse antenna provides a worthwhile degree of 
improvement. In addition, the low spillover and non- 
critical nature of an antenna of this type suggests that 
independent control may permit some unusual antenna 
applications in which monopulse would ordinarily not 
be feasible. 
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A New Way of Solving Maxwell’s Equations* 


V. H. RUMSEY, FELLOW, IRE 


Summary—A general solution of Maxwell’s equations for a single 
frequency can be expressed as the combination of two types of 
solution, each of which is characterized by an electric vector which 
is equal to the magnetic vector times the intrinsic impedance of free 
space but a quarter cycle out of phase. These two types are mathe- 
matically orthogonal, and each can be conveniently expressed by 
means of a single scalar. The method leads to circularly polarized 
surface waves propagating along anisotropic sheets which are per- 
fectly conducting in one tangential direction and perfectly transparent 
in the orthogonal tangential direction. 


I. INTRODUCTION 


HE remarkable development of microwave tech- 
‘lk nology in recent times depends toa great extent on 
a certain way of solving Maxwell’s equations, 
namely the decomposition of the field into transverse 
electric (TE) and transverse magnetic (TM) parts. The 
importance of this technique rests on the fact that it is 
practically the only way of getting a complete set of 
solutions. In this paper we will work out an alternative 
way of solving Maxwell’s. equations which may be 
described roughly asa decomposition of the field into 
right-handed and left-handed circularly polarized parts 
(it is a rough description because it is not literally true, 
but for the present let it serve as a fairly meaningful 
name). 
To understand how this compares with the TE, TM 
method, suppose we take the z axis as the reference so 
E,=0 for the TE part and H.=0 for the TM part. Thus 
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the TM part can be considered as due to electric dipoles 
parallel to the z axis and the TE part as due to magnetic 
dipoles parallel to the z axis. From this point of view, 
the circularly polarized parts can be considered as due to 
right- and left-handed circularly polarized dipoles 
parallel to the z axis, a circularly polarized dipole con- 
sisting of a certain combination of an electric and a 
magnetic dipole. Hence, we can see in a rough way that 
the two methods are equivalent in generality. The 
choice between the two methods therefore depends on 
the nature of the problem at hand. It turns out that 
the circularly polarized method has the advantage in 
connection with problems involving propagation over 
sheets made of perfectly conducting filaments so that, 
at any point on the sheet, the conductivity is infinite in 
one direction but zero in the direction at right angles. 
The simplest example is a sheet made of parallel copper 
wires extremely close together, but not touching, so that 
the sheet would be “invisible” to an electric field 
parallel to the sheet but perpendicular to the wires. 


II. SoLUTIONS FOR wHIcH E= +jnH 


Let us begin with a statement of the classical tech- 
nique (usually called the Hertz potential method). It 
says that Maxwell’s equations, in the form 


VX H=jouek and VX E=—jonH, (1) 
(assuming the e*' time convention) are solved by 
putting 
anGweee—- i "La, (2) 
and tH, = VX 3V; (3) 


= hE 
Fe=V KX 2u 


I 
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where 4 is a unit vector parallel to the z axis and U and 
V satisfy the scalar wave equation 


V+ BY =0 (8? = w'ne). (4) 


It can be shown that if U and V are general solutions of 
(4), the Hertz potential method gives a general solu- 
tion of (1).+? 

As an alternative let us now consider solutions of (1) 
for which E=CH, C being some scalar constant. Sub- 
stitution in (1) shows that this is possible if 


C= + jn (5) 
where 
n= Vu/e. (6) 
Thus we get two types of solution denoted by E; and 
Ee: 
Ey = jnH, and E, = — jnHs. (7) 
Note that the Poynting vector is expressed by 
f= (Re He) x (Re He"), (8) 
which gives 
P, = 7H, X Hi; 


and P», = _— nH, os in hep (9) 


H, and H; being the real and imaginary parts of H, or 
Hy», respectively. Thus the instantaneous Poynting vec- 
tor is independent of time. 
For a lossy medium with conductivity o we have 
jou 
que +o 
IE => Re (n) H, x H; and IPo => Re (n) H, x es 


2 


(10) 


(11) 


While (4) immediately gives us H,; if we know Fi, 
it does not tell us how the three components of F, are 
related. Instead of working with the three scalar func- 
tions which constitute these three components, it is ob- 
viously much better to express the whole field by means 
of only one scalar function. This can be done by apply- 
ing the Hertz potential technique to E;. The results are 
as follows: 


E:=VXVX#Ui— BV X 8U; (12) 
au, 

meilen OY AU Ui, (13) 
ra 


the formulas for E» being obtained by reversing the sign 
of 8. The function U; from which E, is derived is any 
solution of (4). Now let us put 


E= E,+ kp. (14) 


1K. T. Whittaker, Proc. London Math. Soc., vol. 1, pp. 367; 1903. 
2 P. M. Morse and H. Feshbach, “Methods of Theoretical Phys- 
ics,” McGraw-Hill Book Co., Inc., New York, N. Y.; 1953. 
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Then substitution for E; and E2 from (12) gives 


E=VXVX 4&Ui + U2) + BV X &U2 — Ui). 


We see that (U:+ U2) and (U2— Uy) are the TM and 
TE Hertz potential functions for EZ. Now it is evident 
that we can always find U; and U; to fit any such pair 
of Hertz potential functions and therefore in view of the 
generality of the Hertz potential method, we can repre- 
sent any solution of (1) as a linear combination of solu- 
tions of types one and two. In short, the representation 
(14) is complete. 

The orthogonality of the E; and E, types can be 
quickly demonstrated by means of the reaction method.? 
We consider the reaction (or coupling) between source 
a which radiates a field of type Z; and source 6 which 
radiates a field of type E». The reaction, represented by 
the symbol (ab), is expressed by 


(ab) = $f (bs xX H, — Ey X Ha): nds, 
S 


where S is any closed surface separating sources a and 
b, n pointing from a to b. On putting 


(15) 


(16) 


E, = jnH. and E, = — jn, (17) 
we see that (ab) =0, which says there is no coupling or 
reaction between the two sources. (The physical inter- 
pretation of (ab) =0 is this. If source a is a unit current 
generator connected to some antenna which radiates a 
field of type £1, (ab) is the open circuit voltage this an- 
tenna receives from source 6. Thus (ab) =0 means that 
the antenna which radiates type EF, is “blind” to an in- 


cident field of type E».) 


III. SourcEs oF THESE FIELDS 


We saw in (12) that we can express the field of type 
F, in the form 


Hy = VX VX 8, — pv x ee (18) 
Now suppose 
; e-iBr 
U, = rae (19) 


which is the elementary source function. Then we find 
that the second term of H represents the field of an 
electric dipole of current moment 


Il = — 88, (20) 
and the first term represents a magnetic dipole of cur- 
rent moment 


M1 = Sjom. (21) 


*V. H. Rumsey, “Reaction concept in electromagnetic theory,” 
Phys. Rev., vol. 94, pp. 1483-1491; June, 1954. : 
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Thus (18) represents the field of the combination of an 


electric and a magnetic dipole at the same point, aligned 
in the same direction, with 


(M1) = — jn(1l). 


By working out (18) it will be found that such a com- 


(22) 


- bination gives a circularly polarized field at every point 


on the sphere at infinity: it may therefore be called a 
circularly polarized dipole. Note, however, that it does 


_ not give circular polarization in the near field. 


Now we can show that any field of type Z; must be 


due to a collection of these circularly polarized dipoles. 


To do this we write Maxwell’s equations for the field 


_ of an electric dipole density represented by J and a 
magnetic dipole density represented by K: 


WO aT tijuek (23) 
=VXE=K + jopl. (24) 

- Substitution from (4) immediately gives 
Bo ns 5 (25) 


which is precisely the condition (22) for a circularly 
polarized dipole. 

The Huygens source for a plane surface turns out 
very simply in terms of the scalar function U in (18). It 


can be derived by using the method described by 


Rumsey.’ The result is this: If the source of some type 
E, field is on the left of the plane z=0, the field on the 
right is the same as would be generated by a density 
circularly polarized dipoles 
spread over =0, a unit circularly polarized dipole hav- 
ing J/=1 [see (22) ]. 


IV. APPLICATION TO ANISOTROPIC SURFACES 


The kind of application where a field of the type 
E=jnH has advantages is obviously one for which the 


- boundary conditions on £ are the same as the boundary 


conditions on H. Such boundary conditions, although 
not common, occur in connection with propagation over 


- anisotropic surfaces. For example, we may take the 


problem of reflection from a surface on which the normal 
components of E and H vanish, a problem which turns 
out to be a key point in connection with Huygens 
sources consisting of normal dipoles.* Propagation over 
a plane perfectly conducting anisotropic sheet is an- 
other example: a sheet of parallel insulated copper wires 
is a simple case, being perfectly conducting in the direc- 


tion of the wires. In this case E parallel to the wires 


must be zero and H parallel to the wires must be at least 
continuous, since the discontinuity in tangential H is 
perpendicular to the current. Since the conditions on 
either side of the sheet are the same, tangential H on 


“Some new forms of Huygens’ principle,” IRE 


AV skd. R , 
samaed PropaGATION, vol. AP-7, pp. S103-S116,; 
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one side must be equal and opposite to the value on 
the other side. Now H parallel to the wires can only 
satisfy this condition and be continuous if it is zero, so 
here also the boundary conditions on £ and H are the 
same. 

To express this explicitly, suppose that we have a 
field which is of type £; for z>0 and of type E» for 
z<0, E,; and E» being expressed in terms of the.two 
scalars U; and U2 according to (13). Let us make 


ou, 


dU2 ee 
= (26) 


Oz (hee 


i = (0); 


at 


Then it will be found from (13) that this makes tan- 
gential E continuous but tangential H discontinuous 
by the amount J where 
22 x Fy 


If = Tig 36 HEL, ; a) 
jn 
Therefore if 


Eie= 0) vation), (28) 
where / is some vector in the xy plane, J must be parallel 
to I. The physical interpretation of these results is that 
there is an electric current sheet J on z=0 which flows 
along perfectly conducting filaments .pointing in the 
direction of J. One application of this was given by 
Rumsey,® and a simpler application is described in the 
next paragraph. 

To consider the simplest kind of waves that can be 
transmitted along an anistotropic sheet made of per- 
fectly conducting filaments parallel to the x axis we take 


Uy = evertyuutyse (0) (29) 
Uy = — evertw-vez (gz < 0) (30) 
with 
Ve bi bie = —8 (31) 
Note that this fits (26). From (13) we find that 
jn, = Ev = Wy. + 2 X YB + 26"), (32) 
where the vector 
Y= bye t Ivy + ae. (33) 


From (32) we find E-E=0 which means that E (and of 
course H) is circularly polarized. The component of E 
parallel to the filaments is zero in this case if 


= B, = Yes = 08s (34) 


5 \V. H. Rumsey, “A solution to the equiangular spiral antenna 
problem” (Abstract), IRE TRANS. ON ANTENNAS AND PROPAGATION, 
vol. AP-7, p. $117; December, 1959. 
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which gives on substitution in (31) and (32) 
x = Bs = 

, and 
>= lB) = 
Thus the solution is fixed by y,, which may therefore be 


considered to represent the mode in which the anisotropic 
sheet is excited. 

Now Re (£) and Im (£) represent the instantaneous 
fields at =0, and at wt=—7/2. From (29) and (35) 
with 


8 
eh ky 


(8? oe Yy?) U1 
+ jE, 


(35) 


Yu = JB (36) 
we find 
Sums, = Ree, 
= (6? + By”) cos (Byy + Bx)et*u? with z > 0, (37) 
+= Re E, = Im E, = (6? + B,”) sin (Byy + Bx)e*u?, (38) 


from which it can be shown that the instantaneous elec- 
tric field is tangential to the curves 


(39) 


Byz + log sin (v8, + Bx) = constant. 


Fig. 1—The instantaneous field of a circularly polarized 
surface wave traveling over a sheet of wires. 
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Fig. 2—The field of a circularly polarized wave which increases 
exponentially with distance from the anisotropic sheet. 
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These are plotted in Figs. 1 and 2 for x=0. A three- 
dimensional picture is shown in Fig. 3. Note that the 
“phase vector” (the vector normal to planes of constant 
phase) lies in the (x, y) plane, having components 
(+8, B,), but the electric vector lies in the (y, z) plane, 
and consequently the waveform moves obliquely over 
the anisotropic sheet with a phase velocity less than c_ 
(the velocity of light) by the factor (1+ 8,?/6?). Note | 
also that the Poynting vector is along the x axis: . 
P= + £(8? + B,?)?e?u7, (40) 
Thus if we always adopt that choice of the signs in (37) 
to (40) which makes the field decrease exponentially — 
with distance from the anisotropic sheet, reversing the — 
sign of B, reverses the phase vector and the Poynting vec- _ 
tor but does not change the sense of polarization with 
reference to the direction of the Poynting vector, as 
illustrated in Fig. 4. j 


Fig. 3—The waveform moves obliquely over the sheet of 
wires in the direction of the phase vector. 
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Fig. 4—The Poynting vector is inclined to the phase vector, and the 


polarization is right-handed about the Poynting vector for z>0. 
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When 7, is real, corresponding to an exponential driv- 
ing voltage distribution, the solution is essentially sim- 
ilar if we interchange the y and zg axes. However, the 
results look quite different as we can see from Figs. 5 
and 6. Note that the waveform shown in Fig. 6 now 
moves away from the sheet, but still the Poynting vec- 


tor is parallel to the wires and the phase velocity along 


the wires is c. These limitations also apply to the more 
general case where , is complex. In all these cases we 
get slow waves, circularly polarized in the plane normal 
to the wires, and traveling along the wires with the ve- 
locity of light. 

Turning now to a different kind of example let us 
consider a conical surface which is perfectly conducting 
along its generators and perfectly nonconducting in the 
perpendicular direction. Here we can use the radial type 
of Hertz potential. Thus if # is the unit vector in the 
radial direction, let us try 


ou; 
B= (= ) + Fx WU + BHU (41) 
Yr 
outside the cone and 
0U, 
Bs = v( F ) = BF x Ws + BHU (42) 
“fi 


inside. These formulas also give 
k= jn, 


and HE, = — 79H», 


Fig. 5—For an exponential distribution of driving voltage, 
the phase vector is inclined to the wire sheet. 


phase 


cross section 
in YZ plane 


Fig. 6—The waveform moves away from the sheet for an 
exponential excitation along they axis. 
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when U,; and U, satisfy 
U U 
vi— + p— =0. (43) 
r r 


Now if we take (using conventional spherical coor- 
dinates) 


6 n 
tani—— 
U; = eiBr eine (44) 
80 
eye) == 
69 \™ 
tan — 
Ue = e—ibr ; eine | (45) 
ein = 


we find that the radial components of E and H vanish 
everywhere, that Ey is continuous and that there is a 
sheet of current flowing in the radial direction on the 
cone 6=6, of density 


sean aa (46) 
nr sino 

Thus we obtain a set of circularly polarized TEM waves 

traveling along the anisotropic conical surface with a 

phase velocity less than c being equal to c in the radial 

direction. 


V. CONCLUSION 


The purpose of this paper is to introduce a new way 
of solving Maxwell's equations which is more than 
merely an alternative to the classical method, as the 
examples show. Although these examples are simple, 
indeed they are the most elementary of all, they give us 
an interesting class of circularly polarized surface waves, 
which are in some respects strikingly different from 
familiar types of electromagnetic waves. This approach 
suggests a number of further applications, the case of 
lossy media being one of the more obvious. Although the 
mathematics for this case are much the same, the phys- 
ical implications are much different. Another applica- 
tion with extensive ramifications is that mentioned at 
the beginning of Section IV. The application to fre- 
quency independent antennas is not discussed in this 
paper, because it requires explanation at great length, 
but it is perhaps the most interesting, as well as being 
the original, motivation for this investigation. 
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Correlation of Wind Shear with Tropospheric 
Scatter Signals* 


LOUIS H. BAUER}, MEMBER, IRE 


Summary—Examination of the hourly values of the 10-90 per 
cent fading range of selected 915 Mc tropospheric scatter records of 
a 400-mile path have shown unusual behavior which correlates well 
with wind shear through the scatter volume. Although the meteoro- 
logical data were of a gross nature and were collected at a point at 
least 70 miles from the common volume, correlation coefficients of 
0.76 to 0.8 were obtained. The data indicate that the fading range 
(10-90 per cent in this case) increases by 1.4 db per meter/sec/km 
increase in average wind shear through the scatter volume. 


INTRODUCTION 


XPERIENCE and experiment have proved that 
iD radio signals at frequencies of 100 Mc and above 

can be transmitted, reliably, to receiving stations 
located far beyond the horizon. This beyond-the-line- 
of-sight communication, known as tropospheric scatter, 
is a field in which a large amount of time has been de- 
voted to experiment and study during the past ten 
years. The phenomenon was observed much earlier but 
could not be explained satisfactorily.'~* However, as 
the quantity of experimental data increased, it was 
learned that the reception of these signals was not due 
to fleeting, isolated events, and that the signals, cover- 
ing a wide frequency range, could be transmitted and 
received, at will, provided that a transmitting and re- 
ceiving system with sufficient power and sensitivity 
was used. 

In 1950, Booker and Gordon‘ published a theory to 
explain the scatter process by postulating that the elec- 
tromagnetic waves were scattered from irregularities of 
the index of refraction in the atmosphere. The theory 
was based on the premise that there are random varia- 
tions (eddies) of refractive index which cause energy to 
be scattered out of the line of sight toward the receiver. 
Since most of the experimental data could be explained 
by this theory it has gained wide acceptance, although 
additional theories have been advanced at various 
times.* 


* Received by the PGAP, December 28, 1960; revised manuscript 
received, March 15, 1961. The work discussed in this report was 
performed under Contract AF 19(604)-1835 for the Propagation 
Lab. of the Electronics Directorate of the AF Cambridge Research 
Center, Bedford, Mass,, by Cornell Aeronautical Lab., Buffalo, 
INESY. 

{ Radiation, Inc., Orlando, Fla. Formerly with the Cornell Aero- 
nautical Lab., Buffalo, N. Y. 

1A. B. Crawford, C. R. England, and W. W. Mumford, Bell 
Sys. Tech. J., vol. 17, pp. 489-519; October, 1938. 

2R. Hull, OST, vol. 21, pp. 16; May, 1937. 

8G. Marconi, “Radio communication by means of very short 
electric waves,” Proc. Roy. Inst. Great Britain, December, 1932. 

4H. G. Booker and W. E. Gordon, “A theory of radio scattering 
in the troposphere,” Proc. IRE, vol. 38, pp. 401-412; April, 1950. 

5 T. J. Carroll and R. M. Ring, “Propagation of short radio waves 
in a normally stratified troposphere,” Proc. IRE, vol. 43, pp. 1384— 
1390; October, 1955. 


After the Booker-Gordon theory was published, 
scientists tried to correlate meteorological parameters 
with median signal level or fading frequency of the re- 
ceived scattered signal. In this effort, measures of 
turbulence (such as Richardson’s number) were studied 
and analyzed but conclusive results were not produced. 
Some favorable results were obtained when changes in 
refractive index® were correlated with the scattered 
signal. Ringwalt, e¢ al.,7 and Ament* have reported 
correlation of signal level and wind shear at the path 
midpoint at 1250 Mc over a 262-nautical-mile link. 

In this paper, an experiment that produced data 
which show correlation of a measure of turbulence 
(wind shear) with the 10-90 per cent fading range® of 
the scattered signal is described. This effort was moti- 
vated by the desire to explain unusual behavior of the 
received signal for several periods of operation. Al- 
though the meteorological data were of a gross nature, 
conclusive results were obtained. These results should 
provide the impetus for further study and experi- 
mentation. 


PATH GEOMETRY AND EQUIPMENT 


The data in this paper were obtained on an experi- 
mental tropospheric scatter link directed by the Propa- 
gation Laboratory of the Electronics Directorate of the 
Air Force Cambridge Research Center, Bedford, Mass. 
The link was an East-West 400-mile path extending 
from Bedford to Wilson, N. Y., with the midpoint at 
Utica, N. Y. Based on the 3-db points of the free-space 
antenna patterns of the transmitting and receiving 
antennas (28-foot dishes), the common volume extended 
from an altitude of about 20,000 feet (6.0 km) to about 
40,000 feet (12.0 km) over Utica. Fig. 1 shows the 
pertinent path parameters. From the geometry of the 
path, the scattering angle @ was calculated to be ap- 
proximately 2°. 

AFCRC operated and maintained the transmitting 
station which was located at the Bedford airport. The 
transmitter, comprising a temperature controlled crys- 
tal oscillator with a stability of 1 part in 108 per day, 


° M. H. Latour, N. E. Rosier, and W. E. Zetroner, “Some Results 
from an Over-Water Tropospheric Propagation Study,” presented at 
URSI Meeting, Washington, D. C.; April 23-26, 1958. 

7D. L. Ringwalt, W. S. Ament, and F. C. MacDonald, “Measure- 
ments of 1250 Mc scatter propagation as function of meteorology,” 
IRE Trans. on ANTENNAS AND PROPAGATION, vol. AP-6 (Com- 
munication), pp. 208-209; April, 1958. 

*W. S. Ament, “Airborne radiometeorological research,” PRoc. 
IRE, vol. 47, pp. 756-761; May, 1959. 

_ * The 10-90 per cent fading range is defined as the difference in 
signal levels at which the received signal is equalled or exceeded 
10 and 90 per cent of the time. 
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Fig. 1—Plan and side view for Bedford, Mass.-Wilson, N. Y., 
troposphere scatter link. 


was operated at a frequency of 915 Mc CW. Output 
power of this equipment, 10 kw +1 db, was fed into a 
28-foot dish antenna, the center of which was about 25 
feet above ground and oriented toward the western end 
of the scatter link. 

The receiving station located at Wilson was operated 
by Cornell Aeronautical Laboratory under Contract 
AF 19(604)-1835. Equipment in this station included a 
conventional heterodyne-type receiver with an over-all 
bandwidth of 100 cps and a measured noise figure of 


9 db. The local oscillator of the receiver was a frequency 


standard of the same model as the one used in the trans- 
mitter. Consequently, its frequency stability was also 1 


‘part in 10° per day. The resultant drift rate between 


the transmitter and receiver was approximately 10 
cycles per day. Input to the receiver was obtained from 
4 28-foot dish antenna, the center of which was 35 feet 
above ground and oriented toward Bedford. Output of 
the receiver was calibrated with the known power out- 


put of a signal generator. 


The scatter link operated forty-eight hours each week. 


During this time, two types of data were recorded at 


the receiving site. The first consisted of Esterline- 
Angus records which provided a qualitative measure of 


_-the character of the signal, and the second consisted of 


SVS 


records from a ten-level totalizer which produced data 


- suited to statistical analysis. 


REDUCTION OF SCATTER DATA 


The totalizer consisted of ten counters which re- 
corded the time that the signal equalled or exceeded 
preset signal levels. Once each hour the counters were 
photographed and this information was transferred to 
IBM cards for processing on a digital computer. The 
computer determined the percentage of time the signal 
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equalled or exceeded each level and then plotted these 
values on a probability scale from which the median, 
10 and 90 per cent values were obtained. The total error 
due to receiver calibration and method of determining 
the median, 10 and 90 per cent values was estimated to 
be no greater than +3 db. 

The Esterline-Angus records of the received signal 
were used to indicate any unusual signal behavior such 
as large amplitude peaks of short duration (generally 
less than five minutes). These peaks were believed to be 
due to aircraft in or near the common volume. Other 
phenomena such as fading frequency and signal en- 
hancement due to lightning strokes'® were obtained 
from the Esterline-Angus charts. 


WIND SHEAR 


The source of wind shear data was that available in 
Northern Hemisphere Data Tabulations obtained from 
the U. S. Department of Commerce. These tabulations 
report wind speed and direction at 1000-meter altitude 
intervals four times each day. The nearest weather 
station to the common volume was at Albany, N. Y., 
about 70 miles to the east. This situation imposed a 
handicap on correlation of the data because weather 
data were not collected in the common volume. 

The measure of wind shear used for this comparison 
was the mean square value of the average of the wind 
shear where the altitude interval was taken at 1000- 
meter increments for the levels from 7 km to 11 km. 
The altitude interval was limited to this range because 
data outside the 7 to 11-km region showed poor correla- 
tion for some periods while good correlation was always 
obtained for the 7 to 11-km interval. For each 1-km 
layer the shear was found by resolving the vector into 
North-South and East-West components and then add- 
ing these components vectorially. The value of shear 
for the four levels (7-8 km, 8-9 km, 9-10 km, 10-11 
km) was then averaged to obtain a single value. For 
comparison, this value was plotted to the same time 
scale as the fading range. No attempt was made to 
take into account the time delay due to difference in 
location of common volume and weather station. 


CORRELATION OF SCATTER SIGNAL 
AND WIND SHEAR 


Three sets of data containing significant changes in 
fading range are presented. These samples do not, how- 
ever, represent the major portion of the data which has 
been collected since the receiving site has been in opera- 
tion. The major portion of the total data collected on 
this link exhibited normal scatter characteristics with 
the 10-90 per cent fading range values occurring at 
about 13 db. Thus, these samples are unusual and are 
not representative of the majority of the data. 

Wind shear and fading range are shown in Figs. 2-4. 


10 L. H. Bauer and W. A. Flood, “UHF forward scatter from 
lightning strokes,” Proc. IRE (Correspondence), vol. 45, p. 1743; 
December, 1957. 
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Fig. 3—Data for March 6-7, 1957. 


On each graph the hours that are free of peak activity 
due to aircraft have been denoted by circling the point. 
This does not mean that all other hours are worthless, 
because the peaks may be contributing only a small 
amount to the fading range (1 or 2 db). On the other 
hand, the very large increase in fading range for some 
hours is due to excessive peak activity and should not 
be considered because of the difficulty in removing the 
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Fig. 4—Data for March 20-22, 1957. 1 


effects of aircraft. The same behavior is found in the plot _ 
of hourly median values. 7 

It can be seen from the graphs that there is good — 
correlation between the two quantities, fading range — 
and wind shear. To obtain a measure of the correlation, | 
the data were analyzed statistically to obtain the corre- | 
lation coefficient for each period. Briefly, a smooth 
curve was drawn through both wind shear and fading © 
range as shown by the line on the plots. Then the corre- 
sponding values of wind shear and fading range were 
picked off at two hour intervals. The correlation coef- 
ficient was then computed for these values using 


ae n( dy xy) — (Do X09») 
V [nd @) — (x)? 2) — (9)? 


where 


r=correlation coefficient, 
x=value of one of the variables, 
y=value of the other variable, 
n=number of xy pairs. 


No effort was made to compensate for time lag be- 
tween the two sets of data. 

Values of correlation coefficient of 0.8, 0.76, 0.78 for 
February 20-22, 1957, March 6-8, 1957, and March 
20-22, 1957, respectively, were computed. It is ap- 
parent that the correlation is quite good. Confidence in 
the correlation coefficient is contingent on the number 
of independent samples and the accuracy of the data. 
The matter of independent samples is somewhat hard 
to define. The rule used here is that the wind speeds 
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encountered for these data would surely flush the com- 
mon volume in two hours; thus, two-hour intervals 
produce independent samples. If the time constant of 
the mechanism is greater than two hours the choice is, 
‘of course, too small. In the absence of adequate in- 
formation the choice seems reasonable. 

To obtain the behavior of the hourly median value 
‘the variation of the median was subjected to the same 
technique as the fading range. That is, a smooth curve 

-was drawn through the data, and the correlation co- 
_ efficient was computed using two-hour intervals. The 
values of correlation coefficient were 0.79, —0.79, and 
0.46 for February 20-22, 1957, March 6-8, 1957, and 
“March 20-22, 1957. When these values are compared to 
those of fading range for these three samples, it is ob- 
vious that the wind shear does not correlate as well with 
“the median as it did with fading range. This result tends 
to contradict the results of Ringwalt? and Ament’; 
Benowever, additional samples with shear data collected 
at the path midpoint may make the median variation 
_ compatible with their data. 
rk From the smoothed data of Figs. 2—4, the values of 
hourly fading range were plotted against the corre- 
. BP oading value of wind shear as shown in Fig. 5. By the 
method of least squares the trend of the data was 
E established and is indicated by the solid line. The plot 
Behows that the fading range increases by 1.4 db per 
“meter/ sec/km increase in wind shear. It is interesting to 
"note that the value of fading range does not go below 
about 12-13 db, which was also true for the majority of 
all data collected on this link. 
4 The recordings of scatter data during the periods of 
large wind shear were insufficient to determine accu- 
“rately the spectra of fading frequencies. However, by 
4 comparison with samples of known fading characteris- 
tics estimates of the range of fading frequencies were 
_ obtained. For the three samples, frequencies from zero 
- to 5-10 cps were present. This is in contrast to a nom- 
“inal range of zero to about 2 cps. Examination of the 
- meteorological data showed that the winds were pri- 
‘ 
4 


marily from the west, ranging from about 25 to 50 
meters/sec. 

An additional plot of fading and wind shear is shown 
in Fig. 6. For this day (July 3, 1957) a thunderstorm 
moved into the common volume between 1400 and 
1500 hours. Evidence of the thunderstorm activity was 
ee Gsiried in the Esterline-Angus record by spikes due 


b to forward scatter from lightning strokes. From the 
"figure it can be seen that the fading range abruptly 
4 increased by about 13 db between 1400 and 1500 hours. 
Likewise, the average wind shear abruptly increased 
during.this same time. Due to rapid change in the data, 
the correlation of wind shear and fading range is graph- 
ically depicted in the plot. The figure also indicates why 
use of meteorological data at shorter time intervals 


would be most desirable. 
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Fig. 5—Variation in fading range with wind shear. 
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Fig. 6—Data for July 3, 1957. 


CONCLUSION 


The data presented in this paper contain two items 
for consideration. One is the correlation of fading range 
and wind shear; the second, the large departure from 
the expected 13.4 db for the 10-90 per cent fading 
range. The explanations for these two results are not 
readily apparent. One would expect the increased wind 
shear to result in increased turbulence, thus producing 
Rayleigh fading characteristics. However, analysis of 
additional data collected over this link"! indicate that 
signal statistics are seldom stationary over a 15-minute 
interval and almost never for hourly periods. Accord- 
ingly, because of the nonstationary nature of the hourly 
values, sequentially adding several Rayleigh distribu- 
tions of different mean value would be expected to in- 


1 W. K. Klemperer, “Tropospheric Scatter Propagation Investi- 
gation,” Aeronautical Lab., Cornell University, Buffalo, N. Y., 
Rept. No. C-1074-P-1, Final Rept. on Contract No. AF 19 (604)- 
1835; February 19, 1960. 
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crease the fading range to a value greater than that ex- 
pected from a Rayleigh distribution. One might also 
argue, on this basis, that the 10-90 per cent fading 
range for all hourly data would be greater than the 
Rayleigh range of 13.4 db. Such was not the case.” 

The physical mechanism which produced the correla- 
tion between fading range and wind shear is not readily 
apparent from these data. One might speculate that dur- 
ing these periods of high wind shear the viscous cutoff 
of the spectrum of mean square refractive index fluc- 
tuations is located at lower wave numbers as discussed 
by Bolgiano.” Under this condition scattering from 
large eddy sizes could be postulated leading to a quasi- 
coherent scattering mechanism. Whether a model of 
this type is valid can only be determined from additional 
data samples and better meteorological data. 

The results show that additional data of short term 
statistics are required. These are summarized in the 
following list: 


1) Obtain meteorological data in the scatter volume 
and at a number of locations removed from the 


2 R. Bolgiano, “A theory of wavelength dependence in ultra-high 
frequency transhorizon propagation based on meteorological con- 
siderations,” J. Res., vol. 64D, pp. 231-237; May-June, 1960. 
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scatter volume. (These data would be useful in 
determining the horizontal extent of the effects of 
wind shear.) 


September 3 


2) Collect meteorological data at least at hourly in- — 


tervals, and more often if possible, to reduce the 


gross features of present scheduled radiosonde | 


data. 


3) Employ refractometer sounding in and near the © 
common volume. This would furnish small-scale _ 


data. 
4) Include fading frequency spectra in the analysis. 


Inclusion of the above items would remove the gross | 
features of the data of this report and would provide | 
data from which the mechanisms of scatter in the tropo- — 


sphere may be specified. 


It is hoped that the results presented here will pro- — 


vide the stimulus for further analysis of troposcatter 
data to specify the dependence of scatter signals on 
meteorological parameters, particularly wind shear. 


ACKNOWLEDGMENT 


The assistance of Dr. Walter A. Flood in interpreting | 
the data presented in this paper is gratefully acknowl- _ 


edged. 


Coherent and Incoherent Scattering of Microwaves 
from the Ocean* 


C. I. BEARD, SENIOR MEMBER, IRE 


Summary—This report, the third in a series,!:2 summarizes 
experimental studies of microwave over-ocean propagation on line- 
of-sight paths. Measurements at 5.3, 3.2, and 0.86 cm wavelength in 
1955, in the Gulf of Mexico, between two oil drilling platforms one 
mile apart are compared with the results of the phenomenological 
model developed earlier. The total field was measured at maxima, 
and minima of the interference patterns using a sequence of re- 
ceiver beam widths. The scattered field alone was measured by 


* Received by the PGAP, February 27, 1961. This work, com- 
pleted under Signal Corps Contract DA 36-039 SC-85402, was begun 
under Contract NOrd 7386, Bureau of Ordnance, Dept. of the Navy, 
while the writer was with the Appl. Phys. Lab., The Johns Hopkins 
University, Silver Spring, Md. f 

} Sylvania Electronic Defense Labs., Mountain View, Calif. 

1C. I. Beard, I. Katz, and L. M. Spetner, “Phenomenological 
vector model of microwave reflection from the ocean,” IRE Trans. 
ON ANTENNAS AND PropaGaTion, vol. AP-4, pp. 162-167; April, 
1956. Published version of Rept. CM-846, Appl. Phys. Lab., The 
Johns Hopkins Univ., Silver Spring, Md.; July, 1955. 

_ ?C. I. Beard and I. Katz, “The dependence of microwave radio 
signal spectra on ocean roughness and wave spectra,” IRE TRANs. 
ON ANTENNAS AND PRopaGaTION, vol. AP-5, pp. 183-191; April, 1957. 


means of a narrow beam. The 1955, Gulf of Mexico data, are consist- 
ent with the previous 1953, Golden Gate data, in regions of overlap of 
certain parameters, and provide extended coverage for other values 
under different conditions. 

Specifically, some of the new findings are: a) The experimental 
values of the coherent reflected field are larger than exp [—2(2rhy 
/1000)?] for values of ‘‘apparent ocean roughness” (hy/) greater 
than 110 milliradians. h=standard deviation of water surface, 
y=grazing angle in milliradians, and \=electromagnetic wave- 
length. b) After rising from zero, the incoherent scattered power 
exhibits a downtrend for hy/) greater than 110. The incoherent 
power divided by the square of the smooth sea reflection coefficient 
shows no polarization dependence. c) The distribution along the 
surface as measured by the narrow-beam antenna shows that the 
coherent and incoherent power are markedly peaked in the specular 
direction. d) The ratios of coherent-to-incoherent power (m2) ob- 
tained from the shapes of the probability distributions agree with the 
power ratios obtained by measurements of the magnitudes of the 
signals. e) Total signal spectra broaden with increasing hy/) as found 
at the Golden Gate; the relative shapes of spectra agree with those 
from the Golden Gate. Horizontally- and vertically-polarized total 
signal spectra are the same. 


1961 - 
I. INTRODUCTION 


| HIS work is a continuation of previous studies of 
microwave propagation on line-of-sight paths over 
| water. As in earlier work!? the procedure of 
Rice,® essentially as developed by Lawson and Uhlen- 
beck? and Goldstein® and others, is used to analyze 
the measured received fields into coherent (“steady 
signal”) and incoherent (“noise”) components. These 
“components, so isolated, are then plotted against hy/n, 
the “apparent ocean roughness.” h is the standard de- 
‘yiation of the water surface, wy is the grazing angle 
(milliradians in this report), and \ is the electromag- 
“netic wavelength. It is found that the data points clus- 
ter around smooth curves. The Gaussian curve obtained 
_for the coherent field by the theoretical treatment of 
_Pekeris, Macfarlane, Davies,® Goldstein and Gold- 
-muntz,” Ament,® and others, fits the measured field 
‘data only up to AY/X=100. More complete theories 
which treat both the coherent and incoherent scattering 
simultaneously? still have to be applied to the Gaussian 
‘distribution of ocean surface height. The spectral 
"widths of the total received field were experimentally? 
and theoretically!® shown to be related to hy/n. 
_ The reduced data for the 1953 Golden Gate experi- 
ments on forward scatter from the ocean indicated a 
relationship between the fields (coherent as well as in- 
coherent) and In)/X. However, more data were desirable 
_to determine how specific the results were to the condi- 
tions of the Golden Gate, and how dependent the fields 
are on other variables such as ocean-wave direction, 
the number of swells in the first Fresnel zone, wave 
spectra, illumination, polarization, etc. Accordingly, ex- 
' periments were planned at the Applied Physics Labora- 
tory of The Johns Hopkins University, and were per- 
formed by the Electrical Engineering Research Lab- 
oratory (EERL) of the University of Texas in 1955, 
in the Gulf of Mexico. The results of these investiga- 


» 35S. O. Rice, “Mathematical analysis of random noise,” Bell Sys. 

Tech. J., vol. 23, pp. 282-332, 1944; and vol. 24, pp. 46-156, 1945. 
See (3.10)-(11) and (4.2)-(3). 

47. L. Lawson and G. E. Uhlenbeck, “Threshold Signals,” 
_ M.L.T. Rad. Lab. Ser., vol. 24, chs. 3 and 6, McGraw-Hill Book Co., 
* Inc., New York, N. Y.; 1950. 

’ See D. E. Kerr, “Propagation of Short Radio Waves,” M.I.T. 
~ Rad. Lab. Ser., vol. 13, pp. 550-562, McGraw-Hill Book os, Ince, 
> New York, N. Y.; 1951. 

6H. Davies, “The reflection of electromagnetic waves from a 
_ rough surface,” Proc. IEE, pt. 4, Monograph No. 90, pp. 1-6; 

January, 1954. : 
7 H. Goldstein and L. A. Goldmuntz, “Mean Square Signal Com- 
-putation,” Nuclear Development Associates, White Plains, N. Y., 
Rept. NDA-18-2; August 18, 1952. : 
8 W. S. Ament, “Toward a theory of reflection by a rough sur- 
face,” Proc. IRE, vol. 41, pp. 142-146; January, 1953. ' 
9V. Twersky, “On scattering and reflection of electromagnetic 
waves by rough surfaces,” IRE TRANS. ON ANTENNAS AND PROPA- 
GATION, ‘vol. AP-5, pp. 81-90; January, 1957. Also, “On scattering 
and reflection of sound by rough surfaces,” J. Acoust. Soc. Am., vol. 
29, pp. 209-225; February, 1957. 
é 10]. M. Spetner, “A statistical model for forward scattering of 
waves off a rough surface,” IRE TRANS. ON ANTENNAS AND PROPAGA- 
TION, vol. AP-6, pp. 88-94; January, 1958. 
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tions have appeared in a series of reports.!'~ 

The wealth of original data obtained in the course of 
the Gulf of Mexico measurements, and the cooperation 
of individuals at EERL in making data records avail- 
able to the writer, motivated the additional detailed 
analysis!’ reported in the present paper. In order to 
sketch the state of the subject, the results of F. E. 
Brooks!® and the present study are summarized as 
follows: 

1) Vertical polarization Gulf data are given!’ which 
indicate that the coherent field follows exp|—2(2rhy 
/1000A)2] up to hW/d of 80 mils. The coherent term was 
extracted from the data by the approximate equation 
(7) (to be given shortly). The present work, using the 
exact equation (6), shows there is a well-founded ex- 
perimental “curve” formed by both the Gulf and Golden 
Gate data for hy/A>100 which is not accounted for by 
the simple Gaussian theoretical coherent expression 
exp| —2(2rh)/1000A)?].. Shadowing appears to be 
negligible for the data points in this “curve” for hy/d 
above 100. Horizontal and vertical polarization data 
agree in certain cases. Discussion and evidence are 
given of the pronounced effects of surf beats and 
changes of index of refraction gradients on the deter- 
mination of the coherent term from Q-band data. 

2) Vertical polarization Gulf data given'® for the in- 
coherent field for hY/A<80 fill a gap in the earlier 
Golden Gate results.! In the present paper this treat- 
ment is extended, and more Gulf data points for both 
vertical and horizontal polarization are given for 
80 <hy/ <275. For the incoherent field the principal 
result is to show that as a function of hy/d the combined 
Golden Gate and Gulf data points cluster in a smooth 
band that rises from zero, reaches a peak at hy/A~100, 
and then decreases gradually to the highest value of 
hy /d of 275. The peak of the incoherent “curve” occurs 


uF, E. Brooks, Jr., G. P.. DuBose, Jr., and C. W. Tolbert, “Pre- 
liminary Survey Report of the 1955 Gulf of Mexico Propagation 
Tests,” Elec. Engrg. Res. Lab., Univ. of Texas, Austin, Tex., Rept. 
EERL 3-17; January, 1956. 

2 F, E, Brooks, Jr., G. P. DuBose, Jr., and C. W. Tolbert, “Pre- 
liminary Summary Report of Radio, Meteorological and Oceano- 
graphic Data Obtained During 1955 Gulf of Mexico Propagation 
Tests,” Elec. Engrg. Res. Lab., Univ. of Texas, Austin, Tex., Rept. 
EERL 3-18; January, 1956. 

13. C, I, Beard and F. E. Brooks, Jr., “Comparison of Predictions 
of Vector Model of Microwave Reflection from the Ocean with Ex- 
perimental Results,” Appl. Phys. Lab., The Johns Hopkins Univ., 
Silver Spring, Md., Rept. CF-2526; June, 1956. 

“ F, E. Brooks, Jr., N. K. Wagner, and J. R. Gerhardt, “Supple- 
mentary Data Report for 1955 Gulf of Mexico Propagation Tests,” 
Elec. Engrg. Res. Lab., Univ. of Texas, Austin, Tex., Rept. EERL 
3-20; March, 1957. 

6 N. K. Wagner and J. R. Gerhardt, “Final Report of Meteorolog- 
ical and Oceanographic Results for 1955 Gulf of Mexico Propagation 
Tests,” Elec. Engrg. Res. Lab., Univ. of Texas, Austin, Tex., Rept. 
EERL 3-22; March, 1957. 

16 F, E, Brooks, Jr., “Final Report of Radio Results for 1955 
Gulf of Mexico Propagation Tests,” Elec. Engrg. Res. Lab., Univ. 
of Texas, Austin, Tex., Rept. EERL 3-21; March, 1957. 

17C. I, Beard, “Coherent and Incoherent Scattering of Micro- 
waves from the Ocean,” Sylvania Electronic Defense Labs., Moun- 
tain View, Calif., Rept. EDL-E53; July, 1960. 
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at the same value of AY/d (100) at which the experi- 
mental coherent “curve” departs from the Gaussian 
theory; thus, a common mechanism is indicated— 
most probably multiple scattering. Horizontal and 
vertical polarization data are shown to be the same 
within the band of uncertainty. The shape of the 
incoherent curve obtained in this paper is different 
from the earlier work! in which the incoherent term was 
drawn as a constant for hy/X>50. Elementary con- 
siderations are used to show that the initial slope of the 
curve should be zero rather than finite. 

3) The effect of decreasing the receiver illumination 
on the coherent and incoherent terms is small until the 
illumination of the first Fresnel zone is appreciably 
decreased.® In this paper, the trends of the “effective” 
coherent terms are compared with the theoretical 
curve. 

4) The coherent and incoherent field magnitudes 
obtained from a narrow 0.3° beam pointing in the specu- 
lar direction!’ (the 0.3° beam subtends the first Fresnel 
zone) generally fall no farther below the average co- 
herent and incoherent curves than do some of the lower 
data points obtained with antennas that fully illuminate 
the water. From this it is concluded that most of the 
coherent and incoherent energies come from the first 
Fresnel zone. 

5) The present paper develops and analyzes the dis- 
tribution of the coherent and incoherent fields along the 
water surface as seen by the narrow (0.3°) beam. For 
low hw/\ the coherent field is sharply peaked in the 
specular direction, but the coherent “image” spreads 
out for greater roughness (larger hW/d). The incoherent 
power relative to that in the specular direction is also 
peaked in the specular direction, but the surface dis- 
tribution does not change over the range of hW/ meas- 
ured (21 to 78). From these data the differential scat- 
tering cross section can be computed. 

6) By pointing the 0.3° beam in azimuth, it is 
shown" that the scattered energy must be concentrated 
in a narrow band (~0.3°) along the path. 

7) Curves of m? (the ratio of coherent to incoherent 
power) vs if/X are computed from the magnitudes of 
the coherent and incoherent data of the combined 
Golden Gate and Gulf results. These new calculated 
curves fit the data values of m? (obtained from the 
shapes of the probability distributions) at an inter- 
ference maximum more closely than the earlier curves. 
Horizontal and vertical polarization values of m? are 
shown to agree in general in the plot at a maximum. 
The results of m? vs hy/) for the reflected signal alone 
show some agreement. Within the experimental error, 
we infer from the plots of m? that the Rice model applies. 

8) EERL has developed a ten-antenna array at Q 
band with a waveguide switch for rapid sampling of 
each antenna.'® The signals from nine of the antennas 
are resolved into a direct wave, a coherent wave, and 
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an incoherent wave. The arrival directions of the inco- 
herent waves are presented as sources in the water 
called “glitter diagrams”; these apparent sources jump 
about in time. However, in the present paper, the ten- 


polarization data; the coherent and incoherent magni- 
tudes and the total field spectra are compared with 
vertical polarization data. 

9) Cross correlation of the ten-antenna signals in-) 
dicates that the orthogonal components of the inco- 
herent component (to be discussed) are uncorrelated.’*’ 

10) Many examples of the power spectra of the Gulf 
received signals are given.® In the present paper an 
analysis is made of the radio signal spectral widths vs 
hy/X following our earlier procedure.? The resulting. 
plots and those of the earlier work? show similar trends; 
the present results extend the curves to values of 
hy/\ <100. Horizontally and vertically polarized total 
signal spectra are shown to be the same. The amount of 
spectral broadening with hW/X is much less than at the 
Golden Gate, and this may be the result of the equip-— 
ment’s attenuating the higher frequencies. The relative - 
shapes of the Golden Gate and Gulf of Mexico total 
signal spectra are shown to agree, however. 


II. VEctoR MopEL RELATIONS 


The nomenclature and vector diagram have been 
given previously in Fig. 1 of Beard, et al.,1 The en-_ 
velope of the instantaneous total field (of magnitude 7) 
at the receiver is the sum of the assumed constant!® di- 
rect field D, a time-average coherent reflected field (of 
magnitude C and phase ¢), and a random incoherent 
scattered field (of magnitude J). It is assumed that the 
incoherent field J arises from a large number of inde- 
pendent random scatterers and that it is Rayleigh dis- 
tributed with an equiprobable phase. For convenience, 
a steady signal vector S is introduced as S= | D+ C| ; 

The coherent field C and the incoherent power J? are 
extracted from the recordings vs time of the envelope 
of the total field T as follows: 

For this probability distribution, Rice® gives: 


T/o = (m/2)[(1 + m2) Io(m?/2) 
+ m?I,(m?/2)| exp (—m?/2), (1) 
where 


Io{x) = Jo(ix), I\(%) = iJ;(—ix), 


‘8 Tf necessary, the effects of atmospheric inhomogeneities can be 
taken into account by the meteorological results of these Gulf of 
Mexico experiments and other studies. See for example, R. B. 
Muchmore and A. D. Wheelon, “Line-of-sight propagation phe- 
nomena—I. Ray treatment,” Proc. IRE, vol. 43, pp. 1437-1449: 
October, 1955; J. W. Herbstreit and M. C. Thompson, “Measure- 
ments of the phase of radio waves received over transmission 
path with electrical lengths varying as a result of atmospheric tur- 
bulence,” Proc. IRE, vol. 43, pp. 1391-1401; October, 1955. 


and 
m? = S?/20? 
= sine wave power/average incoherent power 
552/74: (2) 
Substituting 
T? = 20°(1 + m?) (3) 
into 
jee ae (4) 
yields 
Talo = [2(1 + m*) — (T/o)?]*”. (5) 


From (1) and (5), Tsa/T can now be computed as a 
function of m*. (See Appendix.) 

Since the cumulative distribution analyzer of EERL 
obtains the median of Z rather than the mean, one 
uses the results of Norton, ef al.,!° which tabulate 
Tmea/S in terms of 1/m? (using our notation). Then by 
use of (2) and (5), one obtains values of T3a/Tmea vs 
m?. (See Appendix.) 

By means of curves drawn from the values in the 
Appendix one can find the magnitude of the vector 
signal Smin=(D—C) at a minimum of the interference 
pattern simply from the measured ratio of T'sa/Tmea or 
T,a/T. The same procedure can also be used to find 
Smax=(D+C) at a maximum. Then C/D is obtained 
from 

C/D ae (Sinax Be isin) / (Omax os ae) (6) 
Note that (6) is not the same as the commonly used 
approximation 


Se Ta — 5 EC Me, ae RE (7) 
obtained directly from the data. 

The incoherent power J?=2c? is obtained from the 
Appendix for the measured ratio Tsa/ Timea (at a maxi- 
mum for best accuracy). 

Although (1)—(5) have been expressed in terms of 
total field 7, they are equally applicable to R, the re- 
flected field alone. In this case, the steady signal S is 
only C, the parameter m? equals C?/20?, and R replaces 
T. With these substitutions, the Appendix can be used 
to solve for Cand o and thus 7°. 

Cumulative distributions for the Rice model,’ ob- 
tained by numerically integrating’® the probability den- 
sity function, are plotted in Fig. 1 asa function of m?. 


19 K, A. Norton, L. E. Vogler, W. V. Mansfield, and P. Jaohort, 
“The probability distribution of the amplitude of a constant vector 
plus a Rayleigh-distributed vector,” Proc. IRE, vol. 43, pp. 1354- 
1361; October, 1955. 
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For the total field data, m? may be expressed in terms 
of measurements at a maximum or a minimum by 


m? = (D + C)?/20% = [(1/r)  (C/Dr)|?/2(¢/Dr)2, (8) 


where the plus sign is for a maximum, the minus sign is 
for a minimum, and 7 is the smooth sea reflection co- 
efficient. Empirical curves of C/Dr and o/Dr vs hW/d 
(to be given later) together with (8) enable one to plot 
m? vs hy/® with ry as a parameter. Factors setting 7 are 
the grazing angle, polarization, wavelength, and con- 
ductivity as influenced by temperature and salinity. 
Thus the plots of m? at maxima and minima (to be 
given later) will apply to a wide variety of conditions 
through the parameter 7. 

If the “reflected” field R is measured solely, then 
instead of (8) one has 


m? = C*/2e? = (C/Dr)?/2(o/ Dr)*. (9) 


Again using the empirical curves of C/Dr and o/Dr vs 
hy/X, one can plot m? as a function of hy /X. 


1 
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Fig. 1—Normalized cumulative amplitude distributions of a con- 
stant vector (coherent signal) plus a Rayleigh-distributed vector 
(noise). m?=ratio of coherent power to incoherent power. 


III]. DESCRIPTION OF GULF EXPERIMENTS 


The site selected was the farthest offshore “pair” of 
oil drilling platforms along the United States coastline. 
The platforms, 5425 feet apart, are located 20 miles off 
the Louisiana coast in the Gulf of Mexico. The large 
offshore distance minimizes the disruptive effect of land 
masses on the uniformity of oceanographic and me- 
teorological conditions over the path area. The path 
length is chosen short to minimize meteorological ef- 
fects, yet long enough to include at least several swells 
so that the statistical assumptions of the model (Sec- 
tion II) will apply. The water depth of about 35 feet 
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is “deep water”2° for wave periods of less than 4 seconds. 
For map and photographs, see Beard.'’ Measurements 
were made from April 19 through May 2, 1955. 

Three wavelengths, C band (5.3 cm), X band (3.2 
cm), and Q band (0.86 cm), were used to give a wide 
range of apparent surface roughness (hy/)). The trans- 
mitting antennas were square 20-db horns for C and 
X bands (with beam widths of about 20° in the vertical 
plane) and a circular 27-db horn for Q band. The re- 
ceiving horns for C and Q bands were the same as their 
transmitting horns. At X band, receiving antennas of 
five different beam widths were used. 

For total field interference pattern measurements, the 
transmitter height was first adjusted for minimum 
average received signal and then a 10-minute recording 
was made; the transmitter was then raised or lowered 
to a position giving maximum average received signal 
and another 10-minute recording made. 

Complete details of the radio measurements are avail- 
able.#!:2.14,16.17 Extensive oceanographic and metero- 


logical measurements” were coordinated with the. 


radio measurements. 


IV. EXPERIMENTAL RESULTS 
A. Coherent Field vs Apparent Ocean Roughness 


In Fig. 2 are plotted the values of the normalized 
coherent field C/Dr obtained in the Gulf of Mexico ex- 
periment (except those in Section IV-C) along with the 
three values from the Golden Gate experiment. (In the 
earlier work! only five points were available; two of 
these were airplane runs of considerable uncertainty 
and, hence, are not repeated here.) For the Gulf data 
essentially only 4/X of hY/d is varying since y changes 
only 1 or 2 mils!” about a value of 13 mils.2! The values 
of C/D were extracted from the data by using curves 
drawn from the table in the Appendix to obtain the 
values for (6). The measured field, from which the co- 
herent field was extracted, is listed in the caption of 
Hagee2. 

Difficulty was experienced with the quantity h/ since 
the step guage was damaged by a floating log at the be- 
ginning of the experiment. A detailed analysis of the 
wave data has been made together with correction for 
the errors introduced." It is estimated that the wave 
heights are now within 10 per cent of their correct 
values, and thus, the values of AW/X are within 10 per 
cent of their correct values. 

The discussion of Fig. 2 is divided naturally into two 


*® Deep water is defined as the depth greater than half the water 
wavelength. \»=(g/27)T? tanh (27D/.) where dx is the wavelength, 
T is the period, g is the acceleration of gravity, and D is the depth. 

*! For ¥=13 mils, the values of the smooth sea, vertical polariza- 
tion reflection coefficient 7, used in calculating C/Dr, are 0.81, 0.82, 
and 0.88 for wavelengths of 5.3, 3.2, and 0.86 cm, respectively. See 
J. A. Saxton and J. A. Lane, “Electrical properties of sea water,” 


Wireless Engrg., vol. 29, pp. 269-275; October, 1952. For the three 
Golden Gate points, y is 27 mils, 
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gnness parame- 

ter (hW/)) is less or greater than 100. Re oe 
For hy/\<100, the Gaussian theoretical curve’* 
given in Fig. 2 generally forms an upper envelope of — 
the Gulf experimental points. This theoretical curve is : 


given by the expression 
C/Dr = exp[—2(2rhp/1000))?], 


parts depending upon whether the roughn 


with y in milliradians. The most likely source of error is 
to be off the true minimum. Whichever way the trans- 
mitter is off the minimum, the measured value of C will 
be less than the true value. Hence, the upper envelope 
of the experimental points (in the absence of other 
errors) would be closest to the true values. 

The result obtained by Rice” for slightly rough sur- 
faces is (in the notation used here) : 


C/Dr = 1 — 2(2rhp/1000N)?. 
This expression fits the data of Fig. 2 up to hW/d of 70 


as well as the Gaussian, since it is the first two terms 
of the expansion of the Gaussian. 
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Fig. 2—Coherent field vs apparent ocean roughness. 


As an example of the errors encountered near a mini- 
mum, the surf beat guage records show that the com- 
bination of surf beat™ plus tide caused water level”! 
variations during the 10-minute recordings at minima 
of up to 1 inch (an average value was 0.6 inch). For 
Q band, a one-inch change in water level over the 
900-foot length of the first Fresnel zone changes the 
relative phase between C and D by 30° for this path. 
The mean total signal for the full 10-minute recording, 
which is the quantity used to compute C/D, is conse- 
quently higher than it would otherwise be, and the 
computed C/D is smaller. 

The index of refraction gradient measured between 


*S. O. Rice, “Reflection of electromagnetic waves from slightly 
rough surfaces,” Commun, Pure Appl. Math., vol. 4, pp. 351-378; 
August, 1951. 

** W. H. Munk, “Surf beats,” Trans. Amer. Geophys. Union, vol. 
30, pp. 849-854; December, 1949. =~? 

* For “shallow water” waves (depth&),/2), \w~Is/gD. For 
the 35-foot depth and T=60 seconds (for surfbeats), \»=2000 feet, 
a sizable fraction of the path. 
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the two heights (of 9.5 feet and 35.4 feet) contained 
swings of 0.4 to 0.8 N unit even when averaged over 
periods of 2 to 3 minutes.” If these large-scale refrac- 
tive-index fluctuations are assumed to travel with the 
mean wind speed, their horizontal dimensions range 
from one to two miles. A change of 0.5 N unit corre- 
sponds to 5.6°, 9.3°, and 34.7° at C, X, and Q bands, 
respectively, if one assumes that this change occurs 
between the entire specular and direct “ray” paths. 

The over-all result of these two effects is that the 
measurements were probably made within + 10° error of 
a minimum at C and X bands. At hW/A=50, the 10° 
error at a minimum corresponds to the measured value 
of C/D being only 4 per cent low. This is a satisfactory 
tolerance. However, the Q-band measurements could 
have varied within +35° of a minimum. At hy/A\ =90, 
a 30° error corresponds to the measured value of C/D 
being 16 per cent low; at hW/AX=250, the 30° error 
causes only an 8 per cent error in C/D. 

That the interference pattern is in motion at times is 
demonstrated by taking selected one-minute averages 
from the two Q-band total field recordings made at a 
one-foot height interval (180)° on April 29; the highest 
and lowest values of C/D are computed from these 
averages by (7) in Table I. 


TABLE I 
Tr C/D 
1-Min. Averages From 1-Min. C/D 
Run —————— Averages ren 
: verage 
Highest | Lowest Highest | Lowest 
Q-I-29a 0.240 OL217 0.21 0.07 0.09 
Q-I-29b 0.187 0.156 


The 3-to-1 variation in C/D shows how the antennas 
are “sliding up and down” the interference pattern. In 
addition, the comparison shows how the 5-minute aver- 
age loses most of the interference pattern. 

An attempt was made to see if there was any de- 
pendence of C/Dr on swell direction. The plot of C/Dr 
vs hy/d for swell direction in a 45° sector centered from 
the south was found to agree with another plot for swell 
direction in a 45° sector centered from the southeast. 


’ Thus with the limited quantity of data available, no 


clear-cut effect exists. 

Fig. 2 also contains the values of C/Dr determined 
from the reflected signal alone (measured by the 0.3°- 
beam width X-band antenna). Five of these seven points 
are near the total signal points; the sixth is not, al- 
though it is next to a total signal point; the seventh is 
in complete disagreement. Since this 0.3° beam sub- 
tends just the first Fresnel zone, it appears that the first 
Fresnel zone contributes the major part of the coherent 
term in five out of seven cases. (The only extremely low 


25 N=(n—1)10%, where 1 is the index of refraction. 
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point at C/Dr=0.18 is suspected of misorientation as 
discussed later in Section IV-D.) 

All of the C/Dr points discussed so far are for vertical 
polarization. One X-band horizontal polarization point 
was obtained when the 0.3°-beam width antenna was 
turned on its side. (The vertical plane beam width is now 
approximately 2°.) With the antenna pointed at the 
transmitter in azimuth, the measured C/Dr=0.91 
(r=1). With the antenna pointed one beamwidth (0.3°) 
off in azimuth, C/Dr =0.90, a good check. These values, 
plotted as one point (at AY/A = 38) in Fig. 2, fall on the 
theoretical curve. 

In contrast to the Gaussian theoretical curve follow- 
ing the upper envelope of the experimental points for 
hy /» <100, it falls markedly below the experimental for 
hy/\>100. The experimental “curve” in Fig. 2 for 
hy/X>100 is formed by a total of seven points: five 
vertical polarization points (two X band from the 1953, 
Golden Gate tests and three Q band from the 1955, Gulf 
of Mexico tests) and two horizontal polarization points 
(from the Q-band 10-antenna data). 

One obvious possibility for the failure of the simple 
Gaussian for AY/A>110 might be shadowing. That is, 
the ocean waves block the microwaves from illuminating 
the troughs of the waves and thus the actual scatterer 
height, say H, is less than the measured height h. Of 
these seven data points, however, four have no shadow- 
ing of the total wave height and two others have neg- 
ligible shadowing (91 per cent of the total wave height 
is seen)”°; only one point has large shadowing.” 

Other data for horizontal polarization were obtained 
by the Q-band 10-antenna system.” Plots of the median 
signal at each of nine vertically spaced antennas are 
shown in Fig. 3 for three of the days on which pro- 
nounced interference patterns were observed. Values of 
C/Dr for all of the seven days of data are plotted in 
Fig. 2. Only two of these seven values agree with the 
other data, and these two also agree with the departure 
of the experimental “curve” above the Gaussian theo- 
retical curve for hy/X> 100. 

Some possible reasons for the 10-antenna values of 
C/Dr being low are as follows: For the geometry used, 
the discrete antenna spacing is one-sixth of a lobe. Thus 
it is possible to be off a minimum by half this spacing, 
or 30 electrical degrees. There appears to be no relation 
between the low values of C/Dr and the swell and wind 
wave directions. Also, shadowing does not seem to be a 
factor since only two of the ten antenna points have any 
appreciable shadowing; 7.e., more than 11 per cent of 


26 B. L. Hicks, “Elementary Theory of Shadowing by a Rough 
Surface,” Control Sys. Lab., Univ. of Illinois, Urbana, IIl., Rept. 
R-81; March, 1956. 

27 A. W. Straiton and C. W. Tolbert, “Measurement and analysis 
of instantaneous height-gain curves at 8.6 millimeters over rough 
surfaces,” IRE TRANS. ON ANTENNAS AND PROPAGATION, vol, AP-4, 
pp. 346-351; July, 1956. 
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Measured Field Polarization Location 


(_] total field vertical Golden Gate 
O total field vertical Gulf of Mexico 
A reflected field alone vertical Gulf of Mexico 
* total field horizontal Gulf of Mexico 
+ total field (10-antenna) horizontal Gulf of Mexico 
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Fig. 3—Sample time average height-gain runs with the 10-antenna, 
Q-band, horizontal-polarization array. (a) April 29, 1955, fitted 
C/D=0.184, hv/X=90; (b) April 30, 1955, fitted C/D=0.08, 
hy /\=237; (c) May 1, 1955, fitted C/D=0.274, hy/r=148. 


the wave height is obscured.?* The fact that the X-band 
horizontal polarization point for hY/X=38, discussed 
previously, agrees with the other data shows that there 
is nothing peculiar to horizontal polarization; rather, it 
also points a finger at the extreme sensitivity of the 
Q-band measurements to water level and meteoro- 
logical changes washing out the interference pattern 
over a several minute time average, as shown previously. 


B. Incoherent Field vs Apparent Ocean Roughness 


The normalized values of the incoherent term o/Dr 
as found from the standard deviation of the fields (Sec- 
tion II) are plotted in Fig. 4 for both the 1953, Golden 
Gate data and the 1955, Gulf data. A smooth curve has 
been drawn through the aggregate of the points with a 
band showing generally the upper and lower limits of 
the data. 

The Gulf total signal data have supplied informa- 
tion on the rising part of the curve below hW/A=50, not 
available in the Golden Gate data. In addition, the 
seven reflected signal points (from the 0.3°-beam an- 


28 These two are the runs at hy/A=274 and 237, both with 40 
per cent of the wave height shadowed, 
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Fig. 4—Incoherent field vs apparent ocean roughness. Incoherent 
power = 20”. 


tenna) lie in the region 20 <M <80. About five out of 
the seven points lie within the band, although six of the 
seven points lie below the “average curve.” Thus in this 
range of hW/d a small area of water around the specular 
direction (approximately the first Fresnel zone) appears 
to be the main contributing area to the incoherent term. 
(Section IV-D also supports this deduction.) One might 
expect a greater surface area to be contributing to the 
incoherent term than to the coherent term. Possibly the 
reason such was not found here is the restricted range 
of hW/d encountered for the reflected signal alone. 

In the Gulf data, the higher values of hy/X (above 
148) were obtained at Q band. The three vertical 
polarization Q-band points, taken at maxima in the 
interference pattern (for hY/AX 2148), fall in and near 
X-band Golden Gate data. Four other Q-band vertical- 
polarization points taken at unknown positions in the 
interference pattern were treated in the following way:"" 
Since D was not known, 7.4 was divided by T instead 
of D and then multiplied by [1+(C/D) | to give the un- 
certainty range. The value of C/D was obtained by 
picking C/Dr from the experimental curve in Fig. 2 and 
multiplying it by r. The uncertainty range is shown by 
the vertical lines on these four points in Fig. 4. 

All of the points just discussed are for vertical polari- 
zation. The X-band horizontal polarization points also 
fall near the average curve. Q-band horizontal-polar- 
ization data from the 10-antenna system are available 
on seven other days. Of these seven points shown on 
Fig. 4, six fall inside the band and one falls slightly be- 
low. 

The accumulation of data at both X and Q bands, at 
different values of h and y, and with both vertical and 
horizontal polarization, supports a downtrend in ¢/Dr 
for hY/A>120. The earlier results! indicated that after 
hy/» had reached a certain value of roughness, greater 
roughness had no more effect, and o/Dr remained con- 
stant. The results of m? vs h)/X to be given later also 
favor the downtrend of o/Dr. It is interesting to note 
that o/Dr begins to decrease from its maximum at 


‘hy/ of 100 to 120, which is the same region at which the 


: 
‘ 
bs 
experimental values of C/Dr depart upwards from the 
Gaussian theoretical curve. Thus a common mecha- 
-nism may be responsible for this behavior. Again, ap- 
‘Plication of multiple scattering’ to this problem is 
needed. 

Although the spread of the experimental data pre- 
cludes the determination of the initial slope of the 
_o/Dr curve as h/A—0, this slope can be arrived at by 
substituting limiting values of the quantities into (8). 
er or the values of parameters involved, the magnitude 
Bot the smooth sea reflection coefficient r behaves as 
; (1—ay) as ~-0. For C/Dr the first two terms of the ex- 
_ pansion of the Gaussian expression are used. Assume 
that o/Dr can be expressed as a function of (hW/\)", 
where 7 is some power, as hW/A-0; then 


4 gio [1/(1 — ap)| + [1 — 2(2rhp/1000n)?] 
@ | V2 f(inp/r)” 


as hy/X— 0. (10) 


J 


eThe limit of m? as hW/A—0 is ~ (see Figs. 10 and 11). 
_ Eq. (10) at a maximum (the plus sign) satisfies this con- 
dition if »21, and is thus not productive. At a mini- 
mum (the minus sign) if M~/A—0 by h-0 or hR>~, 
again n21. However, at a minimum, if hf/A-0 by 
A w—0, then x must be equal to or greater than 2asm’— », 
- Hence, the initial slope of the o/Dr curve instead of be- 
b: ing finite, as drawn previously, should be zero by this 
argument. 
_C. Effect of Illumination on C and o 
To look into the effect of reducing the illumination 
‘of the water surface on C and o, antennas of various 
~ beam widths in the vertical plane were used (Table II). 
5 (The axes of the antenna beams are horizontal.) The 
- results for antenna No. 1 have been given in Figs. 2 and 
_ 4. The results for antennas No. 2~4 are given in Fig. 5. 
_ All data in Fig. 5 are for X band, vertical polarization, 
and an approximately constant grazing angle. The 
abscissa variation is caused primarily by the wave 
' height. 

In the plots of the coherent term, the Gaussian the- 
oretical function from Fig. 2 is plotted along with re- 
~ duced values of this function which are considered to 
approximate the data. The fractional reduction is de- 

noted on the plots. The four low points in the plot for 
antenna No. 3 are ascribed tentatively to not finding the 
~ full depth of the interference minimum. 

In the plots of the incoherent term in Fig. 5, the solid 


a ew 


line is the average line of Fig. 4. In the plots for antennas. 


No. 3 and 4, the ordinates of this line are reduced by the 
same factors found for the coherent terms for compari- 
son only and not to indicate any fit with the data points. 
(Even with all the data in Fig. 4, only a band can be 
drawn.) 
The following qualitative conclusions, however, may 

be drawn from the limited data Ghigo: 
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Measured Field Polarization Location 

(_] total field vertical Golden Gate 
© total field vertical Gulf of Mexico 
A reflected field alone vertical Gulf of Mexico 
* total field horizontal Gulf of Mexico 
+ total field (10-antenna) horizontal Gulf of Mexico 
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Fig. 5—Effect of illumination on coherent and incoherent scattered 
fields. Illumination in specular direction is down <0.2 db, 1 db, 
and 5 db, for antennas 2, 3, and 4, respectively. 


1) The coherent term decreases progressively as the 
illumination decreases. The quantitative decrease, how- 
ever, cannot be described by the quantitative illumina- 
tion decrease at the specular point. Consider Table II. 
Thus, for only one antenna (No. 4) is the coherent term 
cut experimentally by approximately the same amount 
as the beam strength in the specular direction is reduced. 
However, even this appears to be a coincidence, since 
the point on the water at which the illumination is down 
by the same amount that C(expt)/C(theo) is down is 
moving continually toward the transmitter as the re- 
ceiver beam width is narrowed. This is because the por- 
tions of the Fresnel zones between the receiver and the 
specular point are less effective as they are under-illu- 
minated in comparison to the Fresnel zone area between 
the specular point and the transmitter (¢.g., with an- 
tenna No. 4, the specular point illumination is down 5 
db, but the receiver end of the first Fresnel zone is down 
12,db, while the transmitter end is down only 3 db). 
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TABLE II 


COMPARISON OF ILLUMINATION WITH DECREASE OF COHERENT FIELD 


Beam C(expt) Poi 
- : ‘oint on Water at 
ae Vertical Seen at} C(theo) | which Illumination 
An- Beam: aes from | is down by the Ratio 
t width Point. Fig. 5. C(expt) /C (theo) 
ane db down | db down 
1 20% <a? 0 
2 lth <0.2 0.9 receiver end of 12th— 
15th Fresnel zones 
3 Day ~1 Wes) receiver end of first 
Fresnel zone 
4 12? ~S5 4.4 vicinity of specular 
region 


2) The shape of C/Dr vs hp/d approximately follows 
the Gaussian theoretical shape regardless of illumina- 
tion. 

3) The incoherent term is also reduced as the illu- 
mination is decreased. 

4) The effect of decreasing the illumination from the 
antenna at one path terminal is relatively small until 
the illumination of the first Fresnel zone is appreciably 
affected. 


D. Distribution of Coherent and Incoherent Fields Along 
the Water Surface 


The narrow 0.3°-beam width X-band, vertical polariza- 
tion receiving antenna scanned the water surface (by 
steps) to see how the scattered field was distributed 
along the water path. From the 5-minute recording at 
each pointing angle both median and standard devia- 
tion of the reflected field R were obtained. Either the 
coherent field or the directly measured median value 
shows a strong 0.6°-width peak in the specular direction 
(0.8° pointing angle below direct path) followed by a 
gradual decrease to approximately 30 db below direct 
signal at a pointing angle of —7.5°. There seem to be 
no apparent trends in either the median values or the 
coherent fields at angles other than the specular angle. 
However, when the coherent fields (or medians) are 
divided by their values at the specular angle, a trend 
does appear (Fig. 6). First, one notices that at all de- 
pression angles greater than 1.1° the relative values de- 
crease in the order of April 27, 28, 29, and 19, May 1, 
and April 30 (with April 21 not fitting). This sequence 
‘follows the trend of AW/d values (Table III) with the 
exception of April 27. Table III shows that the sequence 
corresponds to increasing values of C/D in the specular 
direction. Values of the relative medians at one angle 
(—3.5°) are given directly under the C/D values in 
Table III to show that the inverse relation holds for 
every day (except April 21). 

One type of error requires mention. (For more detail, 
see Beard.”) All the curves in Fig. 6 should have a 
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Fig. 6—Distribution along the path of the relative coherent field 
received by a 3.3°-beam width antenna. April, 1955: O=19, — 
xX =21, F=27, v=28, v=29, + =30. May, 1955: « =1. q 


TABLE III 
REFLECTED SIGNAL DATA FROM 0.3°-BEAM WIDTH ANTENNA 
(No. 5) 

Item Apr. 27 | Apr. 28 | Apr. 29} Apr. 19 |Apr. 21] May 1 | Apr. 30 
Tmea/D at —0.8°* | 0.200 | 0.458 | 0.499 | 0.506 | 0.535 | 0.570 | 0.821 
Tsa/D at —0.8° 0.101 | 0.219 | 0.162 | 0.151 | 0.161 | 0.066 | 0.068 
«/D at —0.8° ~0.13 | 0.250 | 0.167 | 0.155 | 0.165 | 0.066 | 0.068 
C/D at —0.8° ~0.15 | 0.364 | 0.467 | 0.480 | 0.506 | 0.565 | 0.816 
Tnca( —3.5°) 

— 0.665 | 0.394 | 0.222 | (0.22)+ 0.184 | 0.122 
Tmea( —0.8°) 
hy /» (mils) 54 78 65 70 59 31 21 


* The specular angle is —0.8° below the direct path. 
t+ The value in parentheses is interpolated from Fig. 6. 


maximum of unity at —0.8°. The two days of April 21 
and 27 which have maxima at —0.6° and —0.5°, re- 
spectively, lead one to suspect that the 0.3°-beam- 
width antenna was misoriented by 0.2° and 0.3°, respec- 
tively. (The data of April 26 had to be discarded be- 
cause of a 0.2° misorientation.) However, even if the 
curve of April 27 is reduced by the ratio of peaks of 
1/1.5, the points at —3.5° and —5° are still higher than — 
the respective points for the next-lower curve of April 28. 

These results mean that for large values of C/D (low 
hy/X) the coherent field is sharply peaked and comes 
from the specular direction predominantly. As hy/d 
becomes larger, more and more of the path contributes 
a larger proportion of the coherent field compared to 
that in the specular direction. Stated in other words, as 
hy /® becomes larger, the coherent “image” spreads out. 
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Fig. 7—Distribution along the path of the relative incoherent field 
received by a 0.3°-beam width antenna. April 1955: O=19, X 
=21, L)=27, A=28, V =29, +=30. May 1955: - =1. 


The standard deviations of the signals received with 


the 0.3°-beam width antenna at various depression 


angles are divided by the standard deviations (go) in 


- the specular direction and are plotted in Fig. 7. These 


4 LA) wie 88 


al 


relative standard deviations do not spread out with 
smaller C/D, or larger MW/X, as does the relative coher- 
ent field. Instead they all seem to have the same dis- 
tribution along the path between the specular point 


and the 0.3°-beam antenna (within a band of uncer- 


tainty given by the lines in Fig. 7, the order of points 
within this band shows no trend with C/D or hy/d). 

Thus the incoherent scattered energy, represented by 
2¢2, is also sharply peaked in the specular direction. If 
the ordinates of Fig. 7 are squared and replotted, the 


area under the specular peak is found to be roughly 60 


per cent of the total area. The remaining 40 per cent in 
the “tail” is distributed in the same angular fashion in- 
dependent of M/\ over the narrow range of hy /\ cov- 
ered (21-78). 


E. Azimuthal Distribution of Signal" 


The azimuthal distribution of the received signal was 
investigated by turning the 16-foot receiving antenna on 
its side to put its narrow 0.3° beam in the azimuth plane. 
If the scattered signal from the water surface is dis- 
tributed over an appreciable angle, there should be an 
apparent broadening of the antenna pattern as the 
beam is swung in azimuth. Plots of the median and the 
standard deviation of the received signal are compared 
to the radiation pattern of the antenna in Figs. 8 and 9. 
All values are normalized to unity at 0° when the 16- 
foot antenna is pointed at the transmitter. The filling in 
of the null shows a background of azimuth scattering 
some 15-to-20 db down. The standard deviations show 
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Fig. 8—Azimuthal variation of median signal received by a 0.3°- 
beam width antenna. Receiver position in interference pattern: 
x =at maximum; O=at minimum; []=unknown. 


0.8 
0.7 
ANTENNA 
0.67 PATTERN 
25 
viv 
no) 
mie as 
0.4F 
0.3f 
| fa Oo 
0.2 
R X 
a g Oo 
XN Noo. ae ae 
OL (ee ai Se ——s— | i 
0 0.2 0.4 0.6 0.8 1.0 1.2 


AZIMUTH POINTING ANGLE FROM DIRECT PATH -- DEGREES 


Fig. 9—Azimuthal variation of the standard deviation of the signal 
received by a 0.3°-beam width antenna. Receiver position in inter- 
ference pattern: X=at maximum; O=at minimum; (]=un- 
known. 


the same thing in general, although for 0.6° and greater 
off the path the standard deviation values contain only 
one significant figure. For these conditions (S22. /Cnn, 
horizontal polarization, y=13, hy/N=40) it is con- 
cluded that both the C and J signals come principally 
from a narrow band along the path (less than 0.3° off 
axis on either side of the center line). It would be in- 
teresting to have similar data for large values of hy/r 
and various ocean wave directions to see if the apparent 
antenna pattern broadens significantly. 
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Fig. 10—m? vs hy/) at a minimum in the interference pattern. Verti- 
cal polarization: O=C band, X =X band, and V =Q band. . 


F. First Probability Distribution Shapes vs hy/» 


In the previous sections, the magnitudes of C? and [2 
have been obtained by assuming that the Rice proba- 
bility distribution holds. As a check let us compare the 
shapes of the experimental probability distributions 
with those of the Rice model. A measure of the shape 
is given by the parameter m? and is found by matching 
the experimental plots to the cumulative distributions 
of the Rice model in Fig. 1. Figs. 10-12 show the results 
of this comparison of the values of the shape parameter 
m* with curves computed from the magnitudes of Figs. 
2 and 4 [the solid and dashed curves of Fig. 2 and the 
solid curve of Fig. 4 by means of (8) and (9) ]. On several 
days, an interference pattern at Q band was not de- 
tected by the procedure used."” These points are marked 

with arrows showing the direction of the correction that 
should be made in case the transmitter was not at a 
maximum or a minimum. 

The points at a minimum (Fig. 10) exhibit too much 
scatter for any conclusions to be drawn. The points 
should lie between r=0.81 and r=0.88 for vertical 
polarization at the frequencies used. The higher m? 
values in the region of hW/A<75 can be partially ac- 
counted for by assuming that the measurements were 
not made at the exact minima. (For a 20° error off 
the minimum in this region of hW/X, m® can be twice as 
great.) 


10,000 


1000 


' 


Fig. 11—m? 


| 
: 
a: : ; 
L 
| 
| 


SS os Ss | te 
50 100 150 200 250 300 


hW/X -- MILS 


vs hy/\ at a maximum in the interference pattern. 


Vertical polarization: O =C band, X =X band, and V =(Q band. 
Horizontal polarization: *=X band and 4 =Q band. 


100 


m2 
moo eR or oot 


50 100 150 200 250 300 
hw/dX == MILS 


Fig. 12—m? vs hy/2 for reflected field only. 
Vertical polarization at X band. 
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The experimental behavior at-a maximum (Fig. 11) 
is more reassuring. Scatter at low hW/\ values could be 
Partially a result of wave-height errors since the curves 
slope steeply in this region. j 
_ The fact that the Q-band points for my/A> 100 fall in 
a band between the r=0.8 and 1.0 curves lends inde- 
pendent support to the downtrend of o/Dr for hy/r 
->110 in Fig. 4. The independence arises from deter- 
mining m?’ entirely by the shapes of the probability dis- 
tributions. Curves of m vs hW/d drawn earlier® by tak- 
ing ¢/Dr as a constant! for hY/A> 100, sloped downward 
since the numerator of (8) was continually decreasing. 
Thus, the plot of experimental points, which slopes up- 
ward, did not fit. The downtrend of «/Dr for b/X> 110 
‘in Fig. 4 is sufficient to compensate for the numerator of 
(11) and to produce the agreement of Fig. 11. « 
_ X- and Q-band horizontal polarization data are also 
plotted in Fig. 11 and should lie along the r=1 curve. 
Except for two low Q-band points the agreement is 
good, showing that the presentation of Fig. 11 is gen- 
eral. 

_ The results for reflected signal alone are given in Fig. 
12. Since only X band was used, the data do not extend 
over a large range in hW/d. More extensive coverage in 
_h/) would be valuable since this curve offers one of the 
most critical tests of the model. Except for the low point 
(at hv/\=54) belonging to April 27, 1955 (which run 
‘has previously shown pecularities possibly due to mis- 
orientation), the data show correspondence with the pre- 
dicted curve. Note that m? vs hy /) appears to be insensi- 
tive to polarization from (9) since ¢/Dr has been shown 


to be independent of polarization and C/Dr also appears 


to be insensitive. 


-G. Cross Correlation of Signals from 10-Antenna Array"® 


_ The phenomenological model! assumes that the J, 
and J, components of the incoherent field are inde- 


pendent. If this is true, the cross-correlation function of 


these signals should be zero. Although the 7, and J, sig- 
nals are not available for cross-correlation computa- 

tions, cross correlation between the 10-antenna total 
signals may reveal interesting information. 


Antennas separated by only a slight distance would 


be expected to have high correlation. Antennas sep- 
arated by a quarter of a height gain lobe would be ex- 
pected to have no correlation. Data from three days, 
April 29, April 30, and May 1, were chosen for this 
analysis because these data showed pronounced height- 
gain curves (Fig. 3). For each day cross-correlation 
functions were computed for adjacent antennas, an- 
tennas at a maximum and a minimum, antennas spaced 
a lobe spacing, and antennas at a maximum and at 90°. 
For all cases the cross-correlation coefficient seems to 

be small, being a maximum of 0.4 for adjacent antennas 
on the 29th. (The highest other peaks are 0.3.) Correla- 
‘tion is positive for adjacent antennas, slightly negative 
‘for an antenna at.a maximum and at a minimum in the 
height-gain curve, and very small for antennas sep- 


- 
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arated by a lobe spacing or by a quarter of a lobe spac- 
ing (90°). This would tend to confirm the assumption 
that J, and J, are independent functions. 


H. Relation of Radio Signal Spectra to hv/r 


So far only the amplitude statistics of the radio signal 
have been described. Power spectra of the total signal 
and of the water surface variations have been made”:?93° 
for comparison with each other in the manner of a previ- 
ous investigation.” 

The spectra of vertically polarized signals are ob- 
tained from total signals recorded at the maxima of 
the interference patterns. It has been shown previ- 
ously}? for J,?/| D+ C| 2<<1, which is near a maximum, 
that, to a first approximation, the fluctuating part of T 
is that of Z,. The results are summarized by Beard.” 

Spectra of horizontally polarized signals are obtained 
from the 10-antenna, Q-band system. The nine indi- 
vidual spectra are plotted together on one graph for 
each run.’ Fig. 13 shows one example. 

The fluctuating part of the signal of each of the nine 
antennas (J,) at the various positions in the interfer- 
ence pattern is obtained as follows: 

From the vector diagram of Fig. 1 of Beard, et al.:' 


T? = (D+ Ccos¢+/, cos ¢ — I, sin ¢)” 
+ (C sin ¢ + J, sin @ + I, cos ¢)?. 
S = (D? + 2DC cos ¢ + C?).}/? 
T —S = S[i + (1/S?)(2CI, + 2DI, cos 
— 2DI, sing + 1,2 + 1,7]? — S. 


(eteie rt) KS: 
Ve = Te oS Sa 
~ (1/S,) | I,(Dcos¢+ C) —1,Dsing|+---. 


The spectrum of J, is the sum of the spectrum of J, and 
the spectrum of J, in general. For ¢=0°, only the spec- 
trum of J, is involved; for 6=90° and C<D the spec- 
trum of J, is primarily involved. Hence, if the spectra 
of I, and of I, are the same in shape, then the spectra of 
[, for all nine antennas at a particular time will be sim- 
ilar. Fig. 13 is fairly representative in that some runs 
show more variation and others show less.'* (For six- 
minute runs with 57 degrees of freedom, there is 80 per 
cent confidence that the true value lies between 0.8 and 
1.3 times the actual value.*°) An average is made of the 
nine separate spectra of one run before the 3, +, and zo 
power points are determined. (See Beard" for a sum- 
mary.) 


29 F, E. Brooks, Jr. and H. W. Smith, “Data Reduction Facilities 
at the Electrical Engineering Research Laboratory,” Elec. Energy, 
Res. Lab., University of Texas, Austin, Rept. EERL-81; April, 1956. 

30 R, B. Blackman and J. W. Tukey, “The Measurement of 
Power Spectra,” Dover Publications, Inc., New York, N. ¥. 1959. 
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Fig. 13—Example of power spectrum of total field at nine vertically- 
spaced antennas. April 21, 1955; C/D =0.04; Tsa/D =0.26; hy /X= 214, 


As has been found before, in all but a few cases, the 
peak of the radio spectrum coincides with the peak fre- 
quency of the ocean wave spectrum. The 3, 4, and ro 
power frequencies of the radio spectra are therefore di- 
vided by the peak wave frequencies and plotted as in 
Fig. 14. There are several comparisons to be made be- 
tween Fig. 14 and the similar plot of the previous work.’ 

For the first time both horizontally and vertically 
polarized total signal spectra are plotted together (Fig. 
14), and the results indicate that they are the same. 

The Gulf data have filled in the region of hY/A< 100 
which previously was essentially vacant. As hy/—0 
the 4, 1, and 75 power frequency ratios (f+/fw) seem 
to be leveling off at 1.35, 1.7, and 2.2, respectively. As 
hl /}—0, one would expect the radio spectrum to shrink 
into the ocean wave spectrum. For the Neumann ocean 
wave power spectrum,*! the ratios of $, 4+, and 75 
power frequencies to the peak frequency are 1.47, 1.77, 
and 2.18, respectively. These are to be compared with 
the limiting ratios of the radio spectra given above. 

Within the range of hW/d encountered (up to 275) the 
ratios f,/f. are qualitatively consistent with the linear 
trend found earlier? up to hy/A\=600. However, the 
slopes of the Gulf data in Fig. 14 are only 20 per cent 
of those of the Golden Gate data.’ This situation is 
thought to be a result of the equipment attenuating 
the higher frequencies. (See Beard!” for a detailed dis- 
cussion.) The numerical value of the ordinate scale, 
therefore, has been omitted from Fig. 14. 

Although the magnitudes of the slopes in Fig. 14 do 
not agree with the Golden Gate data, the relative shapes 
of the spectra do agree as shown in Table IV. This fact 
lends considerable support to the basic trends, pointing 
further to equipment limiting the scale factor of Fig. 14. 


31 G, Neumann, “On Ocean Wave Spectra and New Method of 
Forecasting Wind-Generated Sea,” Beach Erosion Board, U. S. Army 
Corps of Engrs, Washington, D, C., Tech. Memo. No. 43; 1953. 
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Fig. 14—Total field power spectral widths vs hy /d. f,(0.1) = frequency 
of 1/10 power point of total field spectrum, etc. fwaves = frequency 
of maximum of ocean wave spectrum. Vertical polarization: 
O=C band, X =X band, and VY =Q band. Horizontal polariza- 
tion: 3+ =Q band. 


TABLE IV 
RELATIVE SPECTRAL SHAPES 
hy / 250 200 150 100 
Ratios GG* | Gulf | GG | Gulf | GG | Gulf | GG | Gulf 
fr(0.1)/f-(0.25) 1:46 | 1.34 [21250 | ges4 aeasels desonlalessnl ies 
fr(0.25) /f,(0.5) 1549) | 4.55" 00048 | a52) 4) 400) £045) | Seale se 
f,(0.1)/f-(0.5) 2.18" ||) 208) | 2.22) || 2.04 | 200. et os) nae 


* Golden Gate data of Beard and Katz.? 
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APPENDIX 


TABLE V 


PARAMETERS OF A PROBABILITY DISTRIBUTION OF A CONSTANT 
SIGNAL VECTOR (S) PLUS A RAYLEIGH-DISTRIBUTED VECTOR 


4] (NoIsE)* 

mm? =—S2/202) T/o T. P K=10-| Tsa/ 
y / / ‘ea/o Tsa/T S/T Tmea/St log1/m?t| Timea 
0.00000 1.2533 | 0.6551 | 0.5227 | 0.00000 0.5564 
~ 0.10000 1.3152 | 0.6857 | 0.5214 | 0.34018 | 2.7653 10 0.5542 
~ 0.15849 1.3507 | 0.7019 | 0.5196 | 0.41682 | 2.2586 8 0.5520 
yi 0.25119 1.4060 | 0.7249 | 0.5156 | 0.50406 | 1.8724 6 0.5463 
ie 0.39811 1.4911 | 0.7568 | 0.5075 | 0.59840 | 1.5860 4 0.5348 
“4 0.63096 1.6205 | 0.7974 | 0.4921 | 0.69320 | 1.3817 2 0.5138 
oe 1. 1.8129 | 0.8446 | 0.4659 | 0.78008 | 1.2437 0 0.4802 
- .1.5849 2.0912 | 0.8925 | 0.4268 | 0.85133 | 1.1546 — 2 0.4342 
es 2.5112 2.4802 | 0.9333 | 0.3763 | 0.90360 | 1.0981 —4 0.3793 
> 3.9811 3.0062 | 0.9616 | 0.3199 | 0.93858 | 1.0622 — 6 0.3208 
o 6.3095 3.6963 | 0.9778 | 0.2645 | 0.96079 | 1.0394 — 8 0.2649 
~ 10.0000 4.5854 | 0.9869 | 0.2152 | 0.97528 | 1.0249 —10 0.2153 
15.848 5.7194 | 0.9920 | 0.1734 | 0.98437 | 1.0157 —12 0.1735 
25.113 7.1579 | 0.9955 | 0.1391 | 0.99010 | 1.0100 —14 0.1391 
~ 39.809 8.97 0.111 1.0063 —16 0.111 
uf 63.091 ils heap Ae 0.0887 1.0040 —18 0.0887 
100.00 14.19 0.0704 1.0025 —20 0.0704 


ie. 


* m2 =signal/noise power. 


+ Reference 19. 
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Electromagnetic Propagation in an Exponential 
Ionization Density* 


LEONARD S. TAYLORT, MEMBER, IRE 


Summary—The propagation of a plane electromagnetic (TE) 
wave into a plane stratified medium in which the ionization density 
varies as exp (z/z») is investigated. The solution of the wave equation. 


_ appears as a combination of Bessel functions of imaginary order and 
complex argument, but the magnitude of the reflection coefficient 


(taken at z= — ~) is given by the simple expression exp { — (4712 
cos 6;/)) tan“(v./w)}, where 0; is the angle of incidence, ) is the 
free-space wavelength and (y./w) is the ratio of electron collision fre- 


_ quency to the frequency of the field. In general, the field components 
~ must be obtained from a complex series, but at depths beyond the 


~ critical density asymptotic forms are given which display the rapid 


decay of the evanescent field in terms of elementary functions. 


INTRODUCTION 


wave in an ionized gas with an exponential varia- 


: | NTEREST in the propagation of an electromagnetic 


tion of electron density has been aroused by the 


- presence of ionization profiles resembling this variation 
_ in the vicinity of reentry and other high-speed bodies in 


the atmosphere. Although numerical solutions of the 
wave equation are sufficient to determine the field quan- 
tities to any desired accuracy, such calculations are not 


) very useful when scaling laws or other general results 


are desired. In the present instance, an extensive numer- 
ical investigation has already been made for the TM 


* Received by the PGAP, March 15, 1961. This work was sup- 
ported by AF Contract 30(602)-1968. 
+ Missile and Space Vehicle Dept., 


Forge, Pa. 


General Electric Co., Valley 


case.! We shall show that an “exact” solution is possible 
for the TE case and that the quantity of principal in- 
terest, the reflection coefficient, may be expressed quite 
generally in elementary terms. Although the evaluation 
of the field amplitudes generally requires the evaluation 
of a complex series, at reasonable depths in the plasma 
asymptotic forms may be used which reduce the series to 
simple forms, displaying the rapid reduction of the 
evanescent wave. 


THE SOLUTION OF THE WAVE EQUATION 


We consider an electrically neutral medium in which 
the electron density N(z) varies exponentially in the 
z-direction, passing through the critical value ma)”/41e” 
(Gaussian Units) at the origin (see Fig. 1). Here m and 
—e are the electronic mass and charge, and w/27 is the 
frequency of the incident electromagnetic wave. The 
index of refraction in the medium is given by 


(1) 


where gp is the “e-folding distance” which characterizes 
the electron density gradient (920), and 


y = (1 — iv-/w). 


n? = 1 — 4rye2N(z)/mw? = 1 — v exp (2/20), 


(2) 


1M. Klein, e¢ al., “Interaction of a Non-Uniform Plasma with 
Microwave Radiation,” General Electric Tech. Information Ser., 
No. R 59SD467; November 23, 1959. 
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47 Ne*/mw- 


Fig. 1. 


y, is the quantity usually described as the collision fre- 
quency. The negative sign appears in (2) by virtue of 
- the choice of the positive exponential for the time de- 
pendence. 

Assuming a TE wave incident at — © in the yz- 
plane, it may be shown from Maxwell’s equations? that 


Eye 2,t). = ul2) exp [i(Ro sindy + wd) ], (3) 
where ky)=w/c and 
d?u(z) ne 
ee + ko?[n?(z) — sin? 0;|u(z) = 0. (4) 
ge 


In these expressions 6; is the angle of incidence. Eq. 
(3) is the result of a separation of variables; the form 
of the y-dependence is an expression of the generalized 
Snell’s Law, k sin 6=const. Combining (1) and (4) we 
obtain 


d?u(z) 
dz? 


+ ko?[cos? 6; — y exp (2/20) |u(z) = 0. (5) 


We now observe that (5) resembles the standard form* 
au 
de Ge [Ro? expi(26) aC? |u a0 (6) 
with solution 
u(é) = Fr(Roe*), (7) 


where & is the cylinder function of imaginary argument. 
Eq. (7) suggests the transformations 


P= 2z0Ro cos 0; (8) 
2zokov/ —Y exp (2/220), (9) 


Qa 


2M. Born and E. Wolf, “The Principles of Optics,” Pergamon 
Press, New York, N. Y. , pp. 50-54; 1959 

2G. N. Watson, “Theory of poser Functions, 2nd ed.,” Cam- 
bridge University Press, Cambridge, England; 1952. 
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which reduce (5) to the Bessel equation 


, ou du ; 2) 
a PRRAGG USES. + v?)p 


I 
= 


(10) 


- 


We write, therefore, 
ula) = AT i(e) + BJ-w(a). (1198 
| 


Eq. (11) represents more than just a formal solution 
of the problem. Although it is true that the Bessel 
functions of imaginary order are not tabulated, the 
theory of Bessel functions (or at least that part of the > 
theory which is independent of the nature of the order) . 
may be applied to the problem. As shall be demon-— 
strated shortly, this procedure leads to elementary ex- 
pressions for some important quantities of interest. 


THE REFLECTION COEFFICIENT 


We adopt the normalization associated with the Bes- | 
sel functions and write the series solution of-(11)# 


(—1)"(@/2)ir#2 
nish ad m!T(m + 1 + i) 
(=1)"(@/2)-#42™ | 
+P ie hISw a2) 


It will be simplest, henceforth, to use a coordinate scale | 
in which z)=1. Then with this convention, from (8) 
and (9), we have 


(—1)™(ko/— 7. 
m!T(m + 1 + iv) 
+ B exp (—izko cos 0;) 


7) iv+2m 
u(z) = A exp (izko cos 4;) > 


m=0 


exp (mz) 


0 (—1)™(ko/ —¥ Y)~ iv+2m 
2» pall Cried Ao See Sash 


(13) 


For large negative z, we have only the first term of each 
series, so that 


u(z) = A exp (izko cos 0,)[(Ro/—y)®/TA + iv) | 


+ B exp (—izky cos 0;)[(ko/ —y)~*/T(1 — iv) ]. (14) 


Thus, at the low density limit, the wave resolves into 
the forward and backward components as required. 
Since we have taken the positive exponential for the | 
time dependence, the second term in (14) is the incident 
wave. Normalizing to unity, we take 


Bi (hon? 9 Pe ta) (15) 


Thus, writing 


A = Rko/—7) Tl + ), (16) 


we have for <0 


u(z) ~ exp (—izko cos 6;) + R exp (izko cos 6,). (17) 
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Therefore, R is the (complex) voltage reflection coeffi- 
cient, and (11) becomes 


u(z) = (ko —7)*T (1 — iv)J_w(a) 


+ RkoV —y) 7°11 + wv) Jo (a). (18) 


We now assume that the plasma is bounded to the 
right by a perfect reflector at z=2z,. Then 


R= (ko —7)?°T(— iv) Jina) /T (iv) F(a), (19) 


where use has been made of the relation T'(é+1) =éT'(é) 
and 
Oy = als = Sy). (20) 


We now move the boundary out to large positive 2, 
‘so that |a,|>>1. Then noting 

q | arg Ow | = 1| arg 7 | = tan i(pe/@) <n, (21) 
we may employ the asymptotic forms? 

| J +iv(aw) = (2/may)!!? cos (aw — 1/4 + imv/2). (22) 


‘Thus, (19) may be written as 


R = (ko —y)2"T(—iv) cos (ay — 1/4 + inv/2)/T (ir) 


© -cos (ay — 7/4 — inv/2). (23) 
‘But, 

ae a /4 ny?) 

cos (a, — 7/4 — inv/2) 

_ exp (—2v/2+ tay —in/4)+exp (rv /2—iow+im/4) (24). 


exp (1v/2+-iay—in/4)+exp (—1v/2—iaw+in/4)- 
Thus, since for _Im(a,,) >0, we have for 2, sufficiently 
large 
cos (a, — 7/4 + imv/2) 


25 
- cos (ay — 7/4 — imv/2) 2°) 


~ exp (7), 


therefore, for distant positions of the reflecting wall 


; R® [ho —7)2*T (—iv)/T(i»)] exp (wr), (26) 


and , does not appear in the expression for Rk. Thus the 
wall has “disappeared” back into the plasma, as ex- 
pected. 
We are interested in the main only in the absolute 
value |R. From (26) 


| R| = | Gov =7)?*| exp (x). (27) . 
But | (pet#) 2” =e 20" whence 
| R| = exp (av — 2y arg /—7). (28) 


54. (2) and (8) now yield the principal result of this 
analysis 


| R| =Sexp [— (4 cos 0;/) tan? (v./w) J. (29) 
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Eq. (29) is a general expression for the reflection co- 
efficient. The only approximations used in its deriva- 
tion involved the position of the reflecting wall in the 
plasma; (29) will be accurate whenever the wall is 
sufficiently deep in the plasma. 

In Fig. 2, the power reflection coefficient | R|? is 
plotted vs cos 0;/X for various values of v./w. (It should 
be recalled that \ will be measured in multiples of 2.) 
Since the TE and TM waves are identical for normal 
incidence, (29) provides a check on numerical calcula- 
tions previously made! for the TM wave. In these calcu- 
lations the wave equation was integrated numerically, 
assuming that the field could be taken as zero slightly 
forward of z=0. The results are in agreement to within 2 
or 3 per cent, confirming the adequacy of the assump- 
tion. We also note in passing that the reflection coeff- 
cient is unity at w=0 and decreases monotonically with 
increasing w, assuming vy,=constant. The limiting forms 


are 
| R| ~ exp (—4r cos 6;/d), y,/'o >> 1 


(v./w)? <K 1 


(30) 
| R| =~ exp (— 2p, cos 6;/z0c), (31) 


where c=3 X10" cm sec7}. 


RP 


0.0l 0.1 10 
(cos 8;)/X 
Fig. 2. 


0.001 


Finally, we notice that since 
a = 2'!ky exp (2/2) X [1 + (-/w)?}-1 
{-[ V1 4 (Jo)? — 1]? + [V1 + G./o)? + 1]*} (32) 


the condition which assured the disappearance of the re- 
flecting wall, Im(a,,)>>1, is given by 


Zy >> 2 1n (A/z), (v./w)? K1 (33) 


and 


Zp o> 2 In [2° @./w) A/a], (fw)? > 1. (34) 


Eqs. (33) and (34) might be regarded as “probes” 
which determine the limits of penetration of the 
radiation field, viz., we regard the reflecting wall as 
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a probe which we insert from + ©. When (33) or (34) 
is fulfilled, an approximate measurement of the reflected 
field amplitude at — will not detect the presence of 
the reflecting wall, so that we might assert that the 
field does not penetrate to the depth of the wall. For 
(v./w)?<1, the measurement would not be very sensi- 
tive to the presence of the perfect reflector in the nearly 
perfectly reflecting medium, but in the lossy case 
(v./w)*>>1, however, the assertion does appear reason- 
able. We notice then that for (v./w)?>>1, the field does 
not appreciably penetrate even up to z=0 for wave- 
lengths which are very short compared to 2p. 

The phase of the reflection coefficient may be calcu- 
lated from (26) if the phase of the gamma function of 
imaginary argument is known. This is a tabulated func- 
tion. Moreover, for short wavelengths, the approxima- 
tion® 


arg {T'(iv)} ~ 3/4 + o[Invy — 1] — 1/127, »>>1 (35) 


may be employed. For long wavelengths, the expres- 
sion® 
: 1-2---(m— 1) 
n 
no 2(2-+ 1) eA ie (gtn— 1) 


(36) 


leads to 


arg { T'(iv) \ 


= 7/2 +yLim {in n—v >> tan om) See Gye) 
n> 0 m=1 

Thus, we have upon expansion of the arc tangent and 

regrouping 

; ; ed (— 1) y2nrt 

arg {T'(iv)} = 2/2—y'v+ ms ag 1), 


m=1 2m 


pel 38) 


where y’=Euler’s number ~0.5772, and ¢ is the Rie- 
mann Zeta function. Thus, ignoring terms in v* or 
higher, we have from (26) 


arg R ~ w + 2r(y’ + In ky), (39) 


It is apparent that in the limit »=0, arg R=7, cor- 
rectly. We point out that (37) may readily be applied to 
reasonably large values of v. Thus the foregoing, plus 
the relation® 


p<omle 


| T'(iv) | = (n/v sinh rp)!/2, (40) 


may be combined with (26), (16), and (15) to evaluate 
the series (12). We forego this general calculation in the 
present work and restrict our attention in the next sec- 
tion to the evanescent fields beyond the critical density 
level. 


*E. J. Gumbel, “Tables of the Gamma Functions for Complex 
Argument,” Natl. Bur. of Standards Publication 9-34; 1954. 

5 E. Jahnke and F. Emde, “Tables of Functions, 4th ed.,” Dover 
Publications, Inc., New York, N. Y., pp. 10-11; 1945. 

§ KE. T. Whittaker and G. N. Watson, “A Course of Modern Anal- 
ysis, 4th ed.,” Cambridge University Press, Cambridge, England, 
pp. 237-239; 1958, 
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ENERGY ABSORPTION 


Further knowledge of the local details of the field may 
be obtained by considering the rate of energy absorp- 
tion in the plasma as a function of z. Since from Max-| 
well’s equations (with w= 1) . 

OB, ey OE 
k 


Oz oy 


ikH =VX E=f , (41) 


the rate of energy absorption is given by ~. . | 


V-S = (c/8r)V- Re { E* X H} 


Of heals | 
= — (c/8rko)V-Im {ies Pe + RE,* rs . (4290) 


From (3) 


V-S = — (c/8rko)V 


Peay s 
-Im {fit sin 6,u*u + ku* =I : (43) 
Z 


The first term on the right now drops out as expected, © 
and we have 


d du 


Employing (1), (2), and (4), we may now show that for 
any general electron number density variation | 
(v-/w) 


ee pel 
327?mw L1 + (»,/w)? 


|r@lual. 4s), 

We observe that for | | >1, it is possible to employ 
the asymptotic forms for the Bessel functions and write _ 
the solution of the wave equation in terms of elementary 
functions. Thus in this case, which is of practical in- 
terest, the field a few e-folding distances beyond the 
critical density may be derived in a simple manner 
from our previous results. 

Combining (18), (22), and (26), we have after some 
recombinations, 


u(z) = iv(2/m)!!?(RoW/ —y)”T(—iv) sinh ay 
{art exp [i(a — 1/4) + rv/ 2}. 
Thus, making use of (32) we find 
C (v-/w) 
<|- “i ao Hoey) 
‘exp (2/2 — 2Ima). 


(46) 


Vi. Geer 


(rv) e™’ sinh rv | 


(47) 
Therefore, the z-dependence of V-S is given by the 


factor 


exp {2/2 — 2%/ko[1-+ V1 + (@/w)?]!2e/2}. (48) 


The very rapid decrease with increasing z is in accord- 
ance with our expectations. For (v./w)?<«1, (47) becomes 


1961 
z 
V:-S = — — (v,/w)|(rv)e™ sinh m| 
821? 


-exp (2/2 — 4koe?!?). (49) 


The bracketed expression yields the 6,-dependence of 
: - S. At grazing angles, the variation is proportional to 
cos? 0;. 

_ The formula (49) is not valid when » is large. The rea- 
son for this discrepancy is that the asymptotic forms for 
the Bessel functions which were employed are not valid 
unless the argument of the functions is large compared 
to the order. In this instance, we must employ the more 
recondite forms, 


Tis(iv sec 8) 
i ~ (2/miv tan B)'/? cos [iv(tan 8 — 8) — 7/4], (50) 
J_i,(iv sec B) 

; ~ (2/miv tan B)'!? cos [iv(tan 8 — B + 7) — 1/4]. (51) 
We justify the use of these expressions on the basis of 
the Meissel expansions? 

: J,(v sec B) = (2/mv tan B)!/2e-?» cos (Q, — 4/4), (52) 
Y,(v sec B) = (2/zmv tan B)'/e-?» sin (Q, — #/4), (53) 
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where 
P= ee (4 sec? 8B + sect B) 
16v? 
cots : 
ere (32 sec? B+---)-: (54) 
cot? 6 
Q, = v(tan 6 — B) — (2-3 Sec: 8) 
24v 
aes ieee ee (55) 
5760v? 


The variables in these expansions are not restricted to 
real quantities. In the present instance we take in (48), 


sec B = a/iv = v/y sec 0; exp (2/2). (56) 


We note that in order to employ (52) and (53) success- 
fully, 8 must be large. Then, as with (46), we now find 


u(z) = v(2/m)!2(ko/ —y)?T (—w)e™ sinh ay 


_ sexp f(s — tan )| 
\ (vy tan )1/? \ 


(57) 


p Some Propagation Characteristics of High UHF Signals 
7 in the Immediate Vicinity of Trees* 


4 A. H. LaGRONE}, FELLOW, IRE, AND C. W. CHAPMAN{ 


Summary—Results are reported of measurements made at very 
low angles of 2880-Mc vertically-polarized signals over wooded 
areas, with the elevation angle to the transmitter the principal vari- 
‘able. The effects of one tree and of many trees on the apparent 
ocation of a signal source, as determined with a narrow-beam an- 

tenna, are reported. A hypothetical direction-finding system is as- 
sumed and its pointing characteristics determined. 


I. INTRODUCTION 


4 HE operation of radio and radar systems at fre- 
quencies in the 3000-Mc range is often limited or 
greatly influenced by the presence of vegetation 

on the path. This is particularly true of a field-tracking 
radar system. This system must, at times, be expected 
to operate effectively under conditions considerably 
less than ideal. Thus, any advance knowledge of the 


* Received by the PGAP, March 24, 1961; revised manuscript 
received, June 7, 1961. This work was done under the sponsorship 
of the AF Research Div., Air Research and Dev. Command, Bedford, 
Mass., Contract AF 19(604)-5504, Task No. S0015-. 

+ Dept. of Elec. Engrg., University of Texas, Austin, Tex. 


effect of vegetation on the path on the performance of 
the tracking system is highly desirable. 

This paper presents measurements of signal strength 
and antenna patterns for various receiver locations be- 
hind trees and for various transmitter elevation angles. 
The effects caused by local vegetation are deduced and 
reported. A hypothetical direction-finding system is as- 
sumed and its pointing characteristics determined. All 
measurements were made at 2880 Mc using vertical 
polarization. 


II]. SIGNAL STRENGTH MEASUREMENTS 


Vertically-polarized signal-strength measurements 
were made with trees on the path very close to the re- 
ceiving antenna. The geometry of the two paths used 
in making measurements is shown in Fig. 1. The Bal- 
cones path employed a 280-foot tower to lower a trans- 
mitter behind a small grove of live oak and cedar trees. 
The transmitting antenna was a wide-angle horn. The 
receiving antenna was a 4-foot parabolic dish having a 
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Fig. 1—Field-measurement paths. (a) Two-mile path at Balcones 
Research Center, Austin, Tex. (b) Path at Bastrop State Park, 
Bastrop, Tex. 


beam angle of 33°. The receiving antenna was located 8 
feet above ground, approximately 500 feet from the edge 
of the grove, and was fixed in elevation pointing at the 
center of the tower. The trees were about 20 feet tall. 
The Bastrop path employed an airplane carrying the 
transmitter with a wide-angle horn mounted on the 
landing gear. The receiving equipment was the same 
as used at the upper path except for a tracking mech- 
anism, added to keep the receiving antenna pointed 
toward the airplane. Tracking was accomplished opti- 
cally to the line-of-sight. For one degree below line-of- 
sight, tracking was accomplished by decreasing the an- 
tenna elevation angle at the same angular rate as ob- 
tained at line-of-sight. The airplane flew at an eleva- 
tion of 1025 feet above local terrain. The pine trees on 
this path were 50-to-60 feet tall. 

It should be noted that the tall trees at Bastrop re- 
sult in a much larger line-of-sight angle a, than at Bal- 
cones. This larger line-of-sight angle, hereafter called 
screening angle, apparently makes a significant differ- 
ence in the signal-strength level when the transmitter is 
below line-of-sight. This is discussed further in a later 
paragraph. 

The results of the signal-strength measurements at 
Balcones and at Bastrop are shown in Fig. 2. Theoretical 
knife-edge diffraction? and smooth-sphere diffraction 
curves’ are shown for comparison. All curves have been 
normalized to their line-of-sight value, and the dis- 
tance factor‘ in the Bastrop data has been removed. 

It will be noted that the signal at Balcones has about 
the same slope as the theoretical smooth-sphere diffrac- 
tion curve, Fig. 2. This is in general agreement with 


1A. Sommerfeld, “Optics,” Academic Press, Inc., New York, 
N. Y., pp. 249-253; 1954. 

? “Summary Technical Report of the National Defense Research 
Committee,” U. S. Gov’t. Printing Office, Washington, D. C., vol. ily 
pp. 125-131; 1946. 

*K. Bullington, “Radio propagation fundamentals,” Bell Sys, 
Tech. J., vol. 36, pp. 593-626; May, 1957. 
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Fig. 2—Signal-strength comparison. 


signal-strength measurements made by other re- 
searchers at different frequencies but under similar cir- 
cumstances.*® Trees are found not to present an opaque 
diffracting edge at the optical line-of-sight, but rather | 
an effective diffracting surface somewhere below the 
optical line-of-sight, the depth depending on the density 
of the foliage and limbs, and the screening angle. : 

The signal-strength curve at Bastrop is found to be 
considerably stronger at the same relative elevation — 
angle than that at Balcones and significantly stronger 
than the theoretical knife-edge diffraction curve. The | 
difference in the screening angle at the two sites is be- 
lieved to be the principal difference in the two cases. 
At Balcones, the propagation path is almost parallel to 
the tops of the trees, while at Bastrop this is not so. At_ 
Balcones, for example, a path from the receiver to the 
transmitter, 1° below the line-of-sight path, passes 
through approximately 1000 feet of oak and cedar 
trees while a similar path at Bastrop passes through 
only 87 feet of pine trees. By this comparison one does 
not mean to imply that the received signal actually 
travels the above paths between the transmitter and 
the receiver, but rather that the diffraction problem is 
different in the two cases. It appears that the effective 
diffracting edge at Bastrop may be well below the 
visual line-of-sight path, while at Balcones it may be 
very near the visual line-of-sight path. More measure- 
ments of the type reported here would have to be made 
before the interrelationship of the variables could be 
put on an empirical basis. 

The time characteristics of the vertically-polarized 
signal at Balcones for different elevation angles of the 
transmitter are shown in Fig. 3. The average wind speed 
was recorded at 10 mph at tree-top level during the 


°H. Head, “The influence of trees on television field strengths 


at ultra-high frequencies,” Proc. IRE, vol. 48, pp. 1016-1020; 
June, 1960. 
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Fig. 3—Time characteristics of signal at Balcones. 
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“measurements. The curves show the signal-strength 
_ variations to increase as the transmitter drops further 
_ below line-of-sight. The variations are observed to be 
_ relatively slow and in the order of 1-to-6 cycles per 
b: minute. 


z: III. RECEIVING ANTENNA RADIATION PATTERNS ON THE 
BALCONES PATH 


Receiving antenna radiation patterns in the horizon- 
tal plane were made at Balcones for the same trans- 
mitter elevation angles shown on the time-characteris- 
tics study in Fig. 3. It was felt that any deterioration in 
the approximate point source, caused by the trees, 
_ would be detected in the patterns by the beamwidth 
_ broadening, the signal maximum shifting, and the nulls 
- filling in. . 

The beamwidth was found to be somewhat greater 
as the transmitter dropped deeper behind the trees, but 
not a significant amount. The signal maximum ap- 
peared, as nearly as could be determined, to come from 
the same point for all cases. No significant fill-in in the 
~ nulls was observed. 


7 
; 


IV. ABSORPTION AND DIFFRACTION STUDY 
OF A SINGLE TREE 


It was decided that a detailed study of a single tree 
be made to see what could be learned about the diffrac- 
tion and absorption of radio signals for such an ob- 
stacle. A preliminary search of the literature revealed 
little information on the subject; however, one signif- 

icant paper by McPetrie and Ford® was found. 

Fig. 4 is a diffraction study of a narrow opaque ob- 
stacle made by McPetrie and Ford. As indicated in the 
sketch of the tree, shown on the figure, the signal- 
strength measurements were made laterally from a 
path passing through the center of the tree, 59 inches 
above ground and approximately 100 feet behind the 


6 J. S. McPetrie and L. H. Ford, “Experiments on propagation of 
9.2 cm wavelengths, especially on the effects of obstacles,” J. Inst. 
Elec. Engrs. (London), vol. 93, pt. 3-A, pp. 531-543; March-May, 

- 1946. 
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Fig. 4—Diffraction study of a narrow opaque obstacle. After J. 3 
McPetrie and L. H. Ford, “Experiments on propagation of 9.2 
cm-wavelengths, especially on the effects of obstacles,” J. Inst. 
ie Engrs. (London), vol. 93, pt. 3-A, pp. 531-543; March—May, 
1946. 


tree. The principal path did not at any time pass 
through the branches of the tree. The trunk was 5 feet 
7 inches in diameter. The Fresnel diffraction curve,® 
shown as a solid line, was computed on a point basis for 
a very long, opaque strip of the same width as the tree. 
The dotted curve is an integrated curve and represents 
the signal received in a 4-foot strip normal to the axis 
of the tree. This curve should approximate the signal 
received by a 120-cm parabolic antenna. The measured 
signals, using the 120-cm parabolic antenna, are shown 
as circles. A reasonable agreement with the theoretical 
curve is obtained. McPetrie and Ford made similar 
measurements at other distances behind the tree with 
comparable results. 

No lateral signal-strength measurements were found 
in the literature where the path passed directly through 
the foliage of the tree. Since this seemed to be im- 
portant to the present study, such measurements were 
made. 

The diffraction and absorption characteristics of a 
signal passing through a large live oak tree, dubbed 
Montopolis oak, is shown in Fig. 5. The oak tree was 
approximately 30 feet tall and had a foliage spread of 
about 50 feet. The trunk was 5.4 feet in diameter, 3 
feet above ground. Three lateral runs were made be- 
hind the tree at distances of 300, 600, and 900 feet. The 
approximate transmission paths through the branches 
of the tree are shown on the tree insert. The approxi- 
mate shape of the foliage is also shown. 

It will be noted that the signal drops a surprising 
amount when directly behind the Montopolis oak and 
considerably more than would be predicted by absorp- 
tion alone. At a frequency of 2880 Mc, the average ab- 
sorption should be approximately 0.18 db/foot, giving 
9 db for the thickest part of the foliage.’ 


7A. H. LaGrone, “Forecasting television service fields,” Proc. 
IRE, vol. 48, pp. 1009-1016; June, 1960. 
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The signal strength is observed to rise sharply about 
18-to-20 feet to the left of the tree trunk. This is due to 
the thinness of the foliage in this part of the tree. The 
transmitter was actually visible through the leaves. 

The terrain behind the Montopolis oak was not table 
top level, hence the height at which the path passes 
through the tree on the various lateral runs was dif- 
ferent. A calibrating run, not shown, was made in front 
of the tree and showed no serious variations in signal 
strength. 

Antenna patterns in the horizontal plane taken di- 
rectly behind the Montopolis oak at distances of 300, 
600, and 900 feet are shown in Fig. 6. The free-space 
pattern of the antenna is included to show the extent of 
pattern distortion. The approximate point at which the 
transmitter would be found if the tree were removed is 
shown on the tree insert. 
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Fig. 5—Diffraction and absorption study at various distances 
behind a live oak tree. 
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Fig. 6—Antenna patterns directly behind Montopolis oak. 
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The deterioration of the approximate point source is" 
clear. The antenna beam patterns are no longer recog- 
nizable, and pointing for a signal maximum would be 
meaningless insofar as indicating the true direction of 
the transmitter. 

The Montopolis oak patterns were all taken with the 
transmitter at a fixed height. No facilities were available | 
at this site for raising or lowering the transmitter more | 
than a few feet at most. In order to vary the height of © 
the transmitter, a tree near the 280-foot tower at 
Balcones was selected. { 

The Balcones tree, a hackberry, was much smaller — 
than the Montopolis oak, but perhaps more nearly — 
represented average field conditions than the Montop-_ 
olis oak. The geometry of the path is shown in Fig. 7. 

Antenna patterns in the horizontal plane taken be- _ 
hind the hackberry tree at Balcones are shown in Fig. 8. 
The patterns were separated in preparing the drawing, 
for the sake of clarity, and their relative positions do — 
not indicate their relative amplitudes. 
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Fig. 7—Path geometry for hackberry tree measurements. 
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Fig. 8—Antenna patterns behind leafless hackberry tree 
at various elevation angles, 
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In these patterns, the deterioration of the approxi- 
mate point source is again evident though not to the 
extent found behind the Montopolis oak (Fig. 6). At 
the highest-elevation angle, ~=9.28°, the shape of the 
beam is reasonably close to the free-space shape, al- 
though some deterioration is observable here. At 
a=9.28°, the center of the main beam is approximately 

'2.85° above the top of the tree. At the lowest-elevation 
angle, a=2.14°, the beamwidth is 94° between half- 
power points. This, compared to the 33°-beam angle for 
the free-space pattern, indicates the distortion taking 
place in the approximate point source. The remaining 
“curves show the pattern constantly changing as the 
transmitter moves up the tower. 


V. PorInTING ERRoR SIGNAL ANALYSIS_ 


6 and 8, on direction-finding radar equipment is severe. 
- Add to this the absorption due to transmission through 
the branches and foliage of the trees, if any, and the 
-problem of low-angle tracking is even more difficult. 


‘4 The effect of pattern distortion, as observed in Figs. 


- Ordinary search radar tracking an airplane through a . 


“tree would see a distorted picture of the target echo. 
The echo would change with elevation angle of the 
target and would be very broad at times. It would be 
_ difficult to get a true bearing on such a target. 

Radars that employ a lobe-comparison technique for 
! improving their pointing accuracy must have lobe pat- 
_ terns that are of a known shape and patterns that do not 
change under operating conditions. Such radars employ 
narrow beams and frequently do their lobe switching 
by employing a nutating feed. A discussion of the lobe- 


- switching techniques employed, of course, is beyond the 


~ scope of this paper. 
_ While it is not desirable to discuss a particular radar, 
a an analysis can be made of a hypothetical radar em- 
_ ploying the 32° pattern of the receiving antenna used 
_ in these measurements. Let it be assumed, for example, 
4 that the antenna beam can be switched continuously 
and rapidly in azimuth through an angle of +23°. 
“On Target” conditions would be indicated by equal 
_ echo signal amplitudes with the beam in the left- and 
right-azimuth positions. With this hypothetical antenna 
- system, the direction of the transmitter was determined 
- for each of the antenna patterns taken at Balcones. 
Fig. 9 shows the theoretical pointing error for four 
patterns measured behind the grove of scrub oak and 
cedar trees on the Balcones path shown in Fig. 1. The 
antenna patterns used in this study are discussed in 


Section III.8 The pointing error for a=+0.26° was. 


assumed to be zero, and the other pointing errors were 
determined relative to this reference. 

Fig. 10 shows the theoretical pointing errors deter- 
mined for the hackberry tree at Balcones. The direc- 


8 A. H. LaGrone and C. W. Chapman, “Some Effects of Trees on 
Low Angle Propagation at a Wavelength of 10.4 Centimeters,” 
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Fig. 9—Theoretical pointing error behind grove of live oak 
and cedar trees at Balcones. 
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Fig. 10—Theoretical pointing error behind leafless 
hackberry tree at Balcones. 


tion obtained from the pattern taken at the highest 
elevation angle (a=9.28°) was assumed to be correct, 
and the other pointing errors were determined relative 
to it. The error curve is observed to be S-shaped, in 
that it shifts from one side to the other as the trans- 
mitter moves down the tower. This pointing error 
would, of course, be different for another tree, and in 
all probability, quite unpredictable. 


VI. CONCLUSIONS 


1) The signal strength for vertical polarization be- 
hind an obstruction of trees is found to vary with the 
screening angle, and for small screening angles, to ap- 
proach the theoretical curve for diffraction over a 
spherical surface. 

2) Time variations in the vertically-polarized signal 
are observed to be relatively slow and of the order of 
1-to-6 cycles per minute. 

3) An approximate point source is distorted when a 
tree or a portion of a tree is intercepted by the main 
beam of the receiving antenna. 

4) Radars employing beam-comparison techniques 
are subject to significant pointing errors when a tree or 
a portion of a tree is intercepted by the main beam of 
the tracking antenna, 
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Tropospheric Scatter Propagation and Meteorological Conditions 
in the Caribbean* 


R. E. GRAY}, MEMBER, IRE 


Summary—Results are given of radio propagation measurements 
which have recently been made in the Caribbean. The transmission 
loss on beyond-the-horizon paths in this region was found to be con- 
siderably less than that measured over paths of similar length in 
temperate climates. 

The fact that the path loss is, in general, relatively low in the 
Caribbean is believed to be chiefly due to the large gradient of the 
refractive index of the atmosphere which normally prevails in the 
first few hundred feet above the earth’s surface in that area. 


INTRODUCTION 


of the radio propagation and meteorological 

conditions in the Caribbean area. More than 
6000 hours of transmission loss measurements have 
been made in the UHF band over eighteen beyond-the- 
horizon paths in this region. This study showed that, 
in general, the median transmission loss was about 14 
db less over typical paths in the Caribbean than on 
overland paths in temperate climates. 

Beyond-the-horizon radio propagation in the Carib- 
bean is influenced by a slow vertical sinking of air from 
anticyclonic areas, and particularly from that area 
known as the Bermuda High, by warm dry air from 
land masses blowing over a cooler sea, and by the rela- 
tively large decrease of the refractive index of the atmos- 
phere with height, a condition which normally pre- 
vails in the first few hundred feet above the earth’s 
surface in this region. 

It is shown that the relatively low path loss measured 
in the Caribbean may be explained by the large refrac- 
tive index gradient normally existing near the earth’s 
surface in that area. The geometry of typical tropo- 
spheric scatter paths, and the nonlinear nature of the 
normal refractive index gradient, indicates that it is the 
gradient in the first few hundred feet, rather than in the 
more usual first kilometer, which determines the modi- 
fied radius of the earth for tropospheric scatter trans- 
mission. 

Measurements show that a single empirical curve may 
be drawn representing the transmission loss as a func- 
tion of the effective distance for both the Caribbean 
and temperate climate areas. 


D evtte the last few years a study has been made 


I. Pato Loss MEASUREMENTS 


Transmission Paths Measured 


This paper describes the results of UHF transmission 
loss measurements made on thirteen beyond-the-hori- 
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zon paths between islands in the West Indies, and on | 
five such paths between Florida and neighboring is- 
lands.'2 These transmission paths varied in length _ 
from 80 to 540 statute miles. 


Conditions of Test 


All transmission loss tests were made on frequencies 
between 800 and 1000 Mc, except those between Flor- 
ida and Nassau,” where a frequency of 1970 Mc was em- 
ployed. Ten-foot in diameter paraboloid antennas were 
used for transmitting and receiving for all tests except 
those on a path longer than 500 miles where a 28-foot 
in diameter paraboloid was used for transmitting; the 
usual 10-foot antenna was employed for receiving. 

Tests on each path were conducted continuously, 24 
hours per day, for periods ranging from two weeks to 
several months; however, the median path loss figures 
referred to in this paper were determined from tests 
lasting a total of between 300 and 800 hours on each 
path. 

The measurements were carried out at various times 
of the year; however, the received carrier level during 
two years of operation on the Florida-Cuba link showed 
that the change in the monthly median path loss in this 
area was only a few db. Summer and winter tests made 
on a longer path in the eastern Caribbean showed a sea- | 
sonal change of 8 db. It is believed, therefore, that the 
measured median loss for all paths may be considered to 
be within about +4 db of the annual median value. 

The recorded level of the received carrier on each 
path was divided into 15-minute intervals, and the 
median carrier level was determined from the distribu- 
tion of the 15-minute medians. It is believed that the 
over-all accuracy of all path loss measurements was of 
the order of +2 db. 


Results of Path Loss Measurements 


The value of the median scatter loss measured on 
each path has been plotted as a function of the effective 
distance. 

By scatter loss is meant the measured median path 
loss in excess of the computed free space loss. The effec- 
tive distance in statute miles is defined as the angular 
distance in radians (angle of intersection of the antenna 
beams) multiplied by the radius of the earth modified 


‘K. P. Stiles, “Report on over-the-horizon radio transmission 
tests between Florida and Cuba,” 1956 IRE ConvENTION RECORD, 
pt. 8, pp. 212-215. 

* K. P. Stiles, “Tropospheric scatter path loss tests Florida- 
Bahamas,” IRE Trans. ON COMMUNICATION Systems, vol. CS-7 
(no. 3.), pp. 205-208; September, 1959, 
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Fig. 1—Scatter loss in the Caribbean area as a function of the 
effective distance computed for a value of K of 1.33. 


3 for normal refraction, which is usually taken as 5280 


- miles. The effective distance of a transmission path is 


thus the distance on the surface of a smooth earth of 
radius 5280 miles which subtends an angle at the center 
of the earth equal to the angle of intersection of the an- 
tenna beams as computed for the actual transmission 
path. It should be noted, however, that the computed 
angular distance for a particular path is a function. of 
the value assumed for the earth’s radius modified for 
atmospheric refraction. 

The value of the median scatter loss measured be- 
tween Florida and Nassau has been reduced by 3 db to 
correct for the fact that this path was measured on ap- 
proximately twice the frequency employed for tests on 
the other paths. 

Scatter Loss and Effective Distance with K = 1.33: The 
values of scatter loss have been plotted in Fig. 1 as a 
function of the effective distance. This distance has been 
computed from the profile of each path, and assumes a 
modified earth radius of 4/3 the true radius, or a value 
of K of 1.33. In Fig. 1 there is also shown an empirical 
curve of the annual median loss at 900 Mc derived from 
the average of a large number of path loss measurements 
made on overland paths in temperate climates. 

Scatter Loss and Effective Distance with K=2.4: In 
Fig. 2 the scatter loss values have been plotted as be- 
fore, as a function of effective distance, but in this case 
the distance has been computed assuming the ratio of 
the modified to real radius of the earth (K) to be 2.4 
instead of 1.33. In Fig. 2, the empirical path loss curve 
for temperate climates is again shown, the effective dis- 
tance in this case corresponding to a value of K of 1.33. 
It will be seen that the scatter loss values, with the ex- 
ception of the four points numbered 1-4, now lie close 
to the temperate climate curve. 

In the case of path No. 1, both horizons were on land 
and it is seen that the measured loss corresponds to an 
effective distance of approximately 220 miles, which is 
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Fig. 2—Scatter loss in the Caribbean area as a function of the 
effective distance computed for a value of K of 2.4. 


‘equivalent to a value of K for this path of about 1.7 in- 


stead of 2.4. 

On path No. 2, a 4000-foot mountain formed a com- 
mon horizon from each end of the path, and it is be- 
lieved that some obstacle gain was obtained on this 
path. 

In the case of path No. 3, a 4000-foot peak on an is- 
land near the middle of the path formed a common hort- 
zon from the two terminals, and considerable obstacle 
gain was no doubt obtained. 

Path No. 4 showed an exceptionally high transmis- 
sion loss, about equivalent to an effective distance cor- 
responding to a value of K of 1.4. The reason for the 
high loss measured on this path is not known, but it is 
suggested that it may have resulted from the fact that 
one end of the path was only about 60 feet above sea 
level, whereas the other was at an elevation of 2000 feet. 


Il. METEOROLOGICAL CONDITIONS 
Refraction of Radio Waves* 


If there were no atmosphere, radio waves would 
travel in straight lines, and the curvature of the rays, 
relative to the earth’s surface, would be 1/R where R is 
the radius of the earth. Refraction in the atmosphere 
causes a bending of the rays, and with a constant re- 
fractive index gradient the curvature is equal to the 
variation in the refractive index with height, that is 
to dn/dh. With a constant refractive index gradient, 
we may consider the radio rays to be straight lines with 
the earth having an effective radius Ri, where 


il i dn 


a a ee 
Riv KR dk 


3D. E. Kerr, “Propagation of Short Wave Radio Waves,” M.I.T. 
Rad. Lab. Ser., McGraw-Hill Book Co., Inc., New York, N. Y., vol. 
152 1951- 


494 
If the ratio of R; to R is denoted by K, then 


1 
K= 


14k dn 
dh 


The Gradient of the Refractive Index 


If one assumes the usual value of K of 1.33 for land 
paths in temperate climates, then transmission loss 
measurements show that the scatter loss on beyond-the- 
horizon paths in the Caribbean area corresponds to an 
effective value of K of 2.4 for oversea paths, and about 
1.7 for overland paths. Actual measurements of the re- 
fractive index gradient in this region, however, indicate 
that the value of K is considerably lower, being approxi- 
mately between 1.45 and 1.65. Fig. 3 shows the average 
value of K at San Juan, Puerto Rico; very similar values 
are also obtained at Swan Island, some 400 miles west 
of Jamaica. At Miami, Fla., the average monthly value 
of K varies from approximately 1.41 to 1.65. The figures 
quoted above were derived from the difference in the re- 
fractive index of the atmosphere measured at the surface 
of the earth and at 1 km above the surface.’ 

The values of K derived in this way are thus in fair 
agreement with the value of 1.7 estimated from path 
loss tests of an overland path, but are much lower than 
the value of 2.4 estimated for oversea paths in the Carib- 
bean. This lack of agreement is believed to be due to 
the fact that the values of K quoted above for San Juan 
and Miami were derived from the gradient of the refrac- 
tive index measured over land for the first km of height, 
and do not, therefore, correspond to the much larger 
gradients which are likely to exist in the first few hun- 
dred feet above the surface of the sea. 


The Geometry of Transmission Paths 


From profiles of typical beyond-the-horizon oversea 
paths, it is seen that a large percentage of the trans- 
mission path is relatively close to the surface of the sea. 
If one assumes a smooth earth with a linear refractive 
index gradient and with antenna elevations of only 30 
feet, then Fig. 4 shows the percentage of the transmis- 
sion path which is below 500 feet for path lengths of 
from 100 to 400 miles and with values of K of 2.4 and 
1.33. If the heights of the antennas are increased to 500 
feet, then the percentage of the path which is below 500 
feet is also increased. For example, if one assumes a 
value of K of 2.4, a path length of 200 miles and the 
terminals 500 feet above sea level, then 84 per cent of 
the transmission path will be below 500 feet. Thus, 
from the geometry of typical beyond-the-horizon trans- 


4 Private Communication from B. R. Bean, Natl. Bur. of Stand- 
ards, Boulder, Colo, 
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Fig. 3—Modified Earth’s radius factor K for San Juan, Puerto Rico. 


PERCENTAGE OF PATH 
BELOW 500 FEET 


200 250 
PATH LENGTH IN MILES 


100 150 300 350 400 


Fig. 4—Percentage of a transmission path below 500 feet with 
antennas 30 feet above a smooth earth. 


mission paths, it is seen that if the refractive index 
gradient is not linear, then it is the average gradient, 
not in the first few thousand feet but in the first few 
hundred feet, that is likely to be of chief importance. 

It can be shown that if the variation in refractive 
index with altitude is parabolic, then the effective gradi- 
ent of the refractive index on a beyond-the-horizon 
transmission path is given by® 


ee 2G, + G2 
3 


where 


G=effective gradient, 
Gi =gradient at the lowest point on the path, 
G.=gradient at the base of the common volume. 


Refractive Index Gradients Near the Surface of the Sea 


The refractive index gradient near the sea surface is 
likely to be high on oversea paths particularly in warm 


» Bs Misme, “Influence Radioclimatiques Sur Les Liasons Trans 
Horizon,” L’Onde Electrique, vol. 394, pp. 116-123; January, 1960. 
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climates. This is so because of the large gradient of the 
relative humidity near the surface of a warm sea, and 
because of the temperature inversion which is often 
caused by warm dry air from neighboring land masses 
blowing over a cooler sea.® 
Measurements made of the refractive index gradient 
over the Irish Sea in the summer of 19447 showed that 
this gradient was greater than normal below 1000 feet 
for wind speeds less than about 18 miles per hour, and 
-under conditions of low wind speed ducting conditions 
often occurred below about 500 feet. 


_ Measurements reported by Katzin and by Bean of the 


refractive index gradient near the sea surface in the 
Caribbean have shown a high percentage of surface 
ducting conditions.®:? Further measurements made in 
-that area by Fitzsimmons and Stephenson’? have 
shown that on the Caribbean side of the Panama Canal 
Zone the large lapse rates of temperature and specific 
humidity, observed in the lowest 50 feet above the 
~ Caribbean, result from the high water temperature of 
80° to 83°F, and from the steady trade winds blowing 
spray into the lower air. Soundings showed that the re- 
_ fractive index gradient in the first 300 feet above the 
_ surface of the sea often showed the existence of a surface 
duct within the first 100 feet with an effective value of 
_K, for the first 300 feet, of about 2.4. 

_ Misme reports measurements made on the south 
- coast of France by observing radar echos from the island 
4 of Corsica. These showed the following seasonal values 
Sof K:4 


ON 


_ Season K 

e . Spring Dick, 
~ Summer 20 
e Fall tae 
' ~ Winter te 
a Annual Dal 


- Thus, measurements of the refractive index near the 
_ sea surface in warm climates indicate that the average 
_ value of K is likely to lie between 2 and 3. 


, 


6 R. B. Montgomery and R. H. Borgoyne, “Modified Index Dis- 
tribution Close to the Ocean Surface,” M.I.T. Rad. Lab., Cambridge. 
_ Mass., Rept.; February 16, 1945. ; 
; 7P. A. Sheppard, “The Structure and Refractive Index of the 
_ Lower Atmosphere,” presented at the Conf. on Meteorological Fac- 
- tors in Radio Wave Propagation, London, Eng.; 1946. 
| 8 M. Katzin, R. W. Bauchman, and W. Binnian, “3 and 9 cen- 
timeter propagation in low ocean ducts,” Proc. IRE, vol. 35, pp. 
891-905; September, 1947. : 
9B. R. Bean, “Climatology of ground-based radio ducts,” J. R. 
Natl. Bur. of Standards, pp. 29-34; July-August, 1959. 
q 10 K. E, Fitzsimmons and S. T. Stephenson, “Low Level Me- 
 teorological Soundings and Radar Correlations for the Panama 
- Canal Zone,” N.D.R.C., Ottawa, Can., project no. PDRC 647; 
June, 1944. ; . 
Pp. Misme, “Le Gradient Equivalent, Measure Directs et 
Calcul Theorique,” Annales des Telecommunication, vol. 15, pp. 92- 
~ 99; March-April, 1960. 
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Refractive Index Gradient Near the Surface of the Land 


It has been shown by Bean and Thayer! that the 
average value of the surface refractivity is well corre- 
lated with the refractive index gradient as measured at 
45 U. S. weather stations. Fig. 5, taken from N.B.S. 
Rept. No. 6025, shows the relation found between the 
surface value of refractive index N,, and AN, the dif- 
ference between the value of the refractive index at the 
earth’s surface and at 1 km above the surface. Since 
the surface value of the refractive index in the Carib- 
bean is usually between 350 and 380, the corresponding 
value of AN would be expected, from Fig. 5, to lie be- 
tween 50 and 60, corresponding to a value of K of be- 
tween 1.5 and 1.65. 


100 
“if liece aoe 


pp lee 


Fig. 5—Refractive index gradient in the first km of height (ANY) 
as a function of the surface value of the refractive index (Ns). 


It has, however, been shown by Misme* that the gra- 
dient of the refractive index within the first 1000 feet of 
the atmosphere may in general be higher than the value 
obtained from Fig. 5. For example, measurements made 
at Trappes, near Paris, showed that the value of AN at 
the 0,35-km level varied from 44 in February to 66 in 
August, the corresponding values of N, being about 310 
and 340. It is probable, therefore, that for overland 
paths in the Caribbean—that is, for paths where the 
two radio horizons are land areas more than about five 
miles from the coast—the average effective value of K 
is somewhat greater than 1.6. This agrees with the figure 
of 1.7 estimated from transmission loss measurements 
made over path No. 1 which had land horizons. 


2 B. R. Bean and G. D. Thayer, “On Models of the Atmospheric 
Radio Refractive Index,” Natl. Bur. of Standards, Rept. No. 6025; 
November, 1959. 

13 P, Misme, “Essai De Radioclimatologie D’Altitude Dans le 
Nord de la France,” Annals des Telecommunications, vol. 13, pp. 
303-310; November—December, 1958. 
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CONCLUSIONS 


Transmission loss measurements made on beyond- 
the-horizon paths in the Caribbean have shown that in 
general the path loss in that area is considerably less 
than that found in temperate climates. These measure- 
ments indicate that if the appropriate values of the 
earth’s radius factor K are assumed, then a single em- 
pirical curve can be drawn, for both Caribbean and 
temperate climate areas, representing the annual me- 
dian scatter loss for effective distances of between 50 
and 300 miles. 

The effective values of K for overland and oversea 
paths are estimated in Table 1. The values of K given 
above for oversea paths are greater than have been gen- 
erally assumed. It is believed, however, that such high 
effective values may be explained by the geometry of 
typical beyond-the-horizon paths, and by the relatively 
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TABLE I 


Value of K 


Caribbean 
tev 
Vee 


Overland paths 


Oversea paths 


large refractive index gradient normally found near the: 
sea surface. 

Transmission loss and meteorological measurementss 
indicate that it is the effective length of the path, cor-: 
responding to the refractive index gradient in the first’ 
few hundred feet of the earth’s atmosphere, and not in! 
the first km, which chiefly determines the median scat-. 
ter loss on typical beyond-the-horizon transmission 
paths. 
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Comments on “Scatter 


Communications with Radar 
Chaff”* 


The recent paper by Hessemer' proposes 


the novel use of (radar) chaff for the purpose 


of scattering electromagnetic energy for non- 
line-of-sight communication links. While in- 
ventive in its approach, this paper makes no 
reference to several well-known papers of 
particular importance to the problem con- 
sidered. 

The omitted references have particular 
bearing on 1) the considerations of scattering 
from an ensemble of chaff dipoles randomly 
oriented in space, and 2) the evaluation of 
the energy scattered from a short chaff di- 
pole. These references should be made in 
order to verify (or refute) the formulas and 
data presented by Hessemer. 

Van Vleck, Bloch, and Hamermesh? have 
written the first definitive paper concerned 
with the scattering of electromagnetic en- 
ergy from thin wires (chaff). The Van Vleck 
work gives the first evaluation of scattering 
by an individual dipole and also evaluates 
the back-scattering cross-section averaged 
over all angles of orientation of the dipole. 
The peak (resonant) value of ¢ (back-scat- 
tering cross section for normal incidence 
averaged over polarization angles) obtained 
by Van Vleck, et al., is about 0.32, meas- 
urably more than the 0.272 obtained by 
Hessemer when his angles are specified to 
correspond to the back-scattering case. 

Further work reported in the open lit- 
erature relating to the cross sections of indi- 


* Received by the PGAP, April 25,1961. 

1R, A, Hessemer, Jr., “Scatter communications 
with radar chaff,” IRE TRANS. ON ANTENNAS AND 
PROPAGATION, vol. AP-9, pp. 211-217; March, 1961. 

2J. H. Van Vleck, F. Bloch, and M. Hamer- 
mesh, “Theory of radar reflection from wires or thin 
metallic strips,” J. Appl. Phys., vol. 18, pp. 274-294; 


March, 1947. 


vidual thin wires,’-® for lengths near the 
first resonance, bear an even greater impor- 
tance on the Hessemer paper. Hessemer uses 
a heuristic approach to this problem, involv- 
ing antenna-dipole input impedances for 
which no terminals are defined. It has been 
shown that reliable results are obtained for 
scattering from thin wires only when a thor- 
ough treatment of the boundary-value prob- 
lem is made. 

Tai? published the first paper using the 
variational technique for treating the 
boundary-value problem posed by the thin- 
conducting wire less than 4.0) in length. Tai 
also discusses the current distributions on the 
wires’ and the validity of the thin wire ap- 
proximation.‘ Kouyoumjian® uses a similar 
technique and compares the results favor- 
ably with experiment. Sevick® obtains ex- 
cellent agreement with experiment, using 
the variational technique with the Morita’ 
current distributions and making corrections 
for end conditions. The writers*:® have ex- 


2C. T. Tai, “Electromagnetic back scattering 
from cylindrical wires,” J. Appl. Phys., vol. 23, pp. 
909-916; August, 1952. 

4C, T. Tai, “A new interpretation of the integral 
equation formulation of cylindrical antennas,” IRE 
TRANS. ON ANTENNAS AND PROPAGATION, vol. AP-3, 
pp. 125-127; July, 1955. 

5R. C. Kouyoumjian, “The Calculation of the 
Echo Areas of Perfectly Conducting Objects by the 
Variational Method,” Antenna Lab., The Ohio State 
University, Columbus, Tech. Rept. No. 444-13; No- 
vember, 1953. 

6 J. Sevick, “Experimental and Theoretical Re- 
sults on the Back Scattering Cross Section of Coupled 
Antennas,” Cruft Lab., Harvard University, Cam- 
bridge, Mass., Tech. Rept. No. 150; May, 1952. 

77T. Morita, “The measurement of current and 
charge distributions on cylindrical antennas,” Proc. 
IRE, vol. 38, pp. 898-904; August, 1950. 

8 EB. S. Cassedy and J. Fainberg, “Back scattering 
cross sections of cylindrical wires of finite conduc- 
tivity,” IRE TRANS. ON ANTENNAS AND PROPAGA- 
TION, vol. AP-8, pp. 1-7; January, 1960. 

9B, S. Cassedy and J. Fainberg, “Electromag- 
netic Cross Sections of Cylinders of Finite Conduc- 
tivity,” Rad. Lab., The Johns Hopkins University, 
Baltimore, Md., Tech. Rept. No. AF-81; August, 


1960. 
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tended the above results to include effects 
of loss, obtaining good comparison between 
theory and experimental results. 

An example indicating the error of the 
approximate method of Hessemer may be 
drawn from the authors’ work.® A dipole of 
length 0.4 has a back-scattering cross sec- 
tion (per unit length of dipole) which is re- 
duced by a factor of 1.74 from that of a 
resonant dipole. This dipole is a circular 
cylinder about 0.02 in diameter, which is 
comparable to the flat strip of a length-to- 
width ratio of 20 to 35 given by Hessemer. 
Hessemer obtains a reduction of 6.5 for his 
0.4 dipole. 

E. S. CASSEDY 

Microwave Research Inst. 
Polytechnic Inst. of Brooklyn 
Brooklyn, N. Y. 

J. FAINBERG 

Radiation Lab. 

The Johns Hopkins University 
Baltimore, Md. 


Author's Reply” 


The comments of Cassedy and Fainberg 
are well taken but need to be put into proper 
perspective with reference to the work under 


discussion. 
The author was well aware of the work 


of Van Vleck, Bloch, Hamermesh, and 
Phillips" and made reference to it in an 


10 Received by the PGAP, May 24, 1961. 

11, Bloch, M. Hamermesh, and M., Phillips, 
“Return Cross-Sections From Random Oriented Reso- 
nant Half-Wave Length Chaff,” Radio Research 
Lab., Harvard University, Cambridge, Mass., 411- 
TM-127; June 19, 1944. 

2 J, H, Van Vleck, F. Bloch, and M, Hamermesh, 
“Theory of radio reflection from wires or thin metallic 
poe J. Appl. Phys., vol. 18, pp. 274-294; March, 
1947, 
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earlier report. However, the sophisticated 
method that they used does not lend itself 
to a rotation of coordinates, which is the 
central point of the paper in question. In 
fact, very little work is reported in the lit- 
erature for the case where the transmitter 
and receiver are at different locations. As to 
the values for the singular back-scatter case, 
Mentzer shows values ranging from a low 
of 0.242 to a high of 0.37d2, depending on 
the method of analysis. It is rather significant 
that a direct “antenna” approach yields an 
answer in this same range as well as the very 
important “bistatic” data. 

The second point seems to be more a 
matter of interpretation. A half-wave reso- 
nant dipole actually resonates at a length 
slightly less than 0.5\ depending on the an- 
tenna geometry. The data presented by the 
author assumed the resonant point at 0.5); 
thus the 0.4 data represent a point 0.1A 
down from resonance. With this interpreta- 
tion, the work of Cassedy and Fainberg”® 
gives a value of about 5.5 rather than 1.74. 
Considering the differences in length-to- 
width ratio plus the fact that the numerical 
values given by this author were merely 
approximations used to demonstrate a re- 
sult, it is difficult to compare the values of 
5.5 and 6.5. 

R. A. HESSEMER, JR. 
Elec. Engrg. Dept. 
University of Arizona 
Tucson, Ariz. 


1% R. A. Hessemer, Jr., “Scatter Communications 
with Radar Chaff,” Appl. Research Lab., University 
of Arizona, Tucson, S-5-80146; March, 1959. 

uJ. R. Mentzer, “Scattering and Diffraction of 
Radio Waves,” Pergamon Press, Inc., New York, 
N. Y., p. 98; 1955. 

ib EK, S. Cassedy and J. Fainberg, “Back scattering 
cross sections of cylindrical wires of finite conductiv- 
ity,” IRE Trans. ON ANTENNAS AND PROPAGATION, 
vol. Ap-8, pp. 1-7; January, 1960. 


Calculated Equatorial Plane 
Radiation Patterns Produced by a 
Circumferential Slot on a Cylinder* 


Summary—The equatorial plane radi- 
ation fields produced by a sinusoidal distri- 
bution of axial electric field along a half 
wavelength circumferential slot on a metallic 
cylinder of radius a have been computed by 
means of an IBM-704 electronic computer. 
Radiation patterns in this plane are given for 
sizes of cylinder from 1<ka<30. It is 
believed that these computations are the 
most accurate made to date for these pat- 
terns. From these patterns it is seen that as 
the cylinder increases in size, for a given 
wavelength, the patterns sharpen, and in the 
limit of infinite cylinder radius the patterns 
approach that from a slot on an infinite plane. 


Recent work involving the determination 
of the radiation fields produced by a circum- 


* Received by the PGAP, May 1, 1961. 
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ferential slot on a dielectric coated cylinder! 
considers, as a special case, the noncoated 
cylinder. It is believed that the calculations 
made for this special case are the most ac- 
curate made to date since they were per- 
formed by means of an IBM-704 computer. 
Although some of these calculations were 
manually performed in the past by others 
and have been reported,? it seems worth- 
while to report these more accurate results 
here in a concise form over a large range of 
cylinder radii. The field calculated Eg is given 
by (131) p. 44 of Wait? for the equatorial 
plane 6=7/2 with the slot centered at ¢o=0. 


Power radiation patterns in decibels, 
plots of 


1C, M, Knop, “The Radiation Fields from a Cir- 
cumferential Slot on a Metal Cylinder Coated with a 
Lossy Dielectric,” Res. Lab., Bendix Systems Divi- 
sion, Bendix Corp., Ann Arbor, Mich. Research Note 
No. 21; October 5, 1960. 

2J. R. Wait, “Electromagnetic Radiation from 
Cylindrical Structures,” Pergamon Press, New York, 
N, Y., p. 3; 1959, 
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On Axis FIELD STRENGTH VALUES 


F(o) = IF (6) |e 


c=ka 


COIAMPANE 


are shown in Fig. 1, where f(¢) is defined by 


|F(b) |p—o OF |p=o Terms 
0.31677 296.90 5 
0.32480 6.7027 10 
0.32349 69.231 15 
0.32214 129.39 12 
0.32118 188.50 15 
0.32052 247.05 18 
0.32005 305.25 21 
0.31972 3.2459 18 
0.31947 61.090 27 
0.31927 118.82 30 
0.31901 234.09 30 
0.31884 349.16 38 
0.31872 104.12 42 
0.31864 219.00 46 
0.31858 333.82 50 
0.31854 88.605 54 
0.31850 203.35 58 
0.31847 318.08 62 
0.31845 72.78 66 
0.31843 187.47 70 
0.31831 —90.000 Plane 
Solution 
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and where Vo=voltage across the center of 
the slot,  =27/do, and \y)=free space wave- 
length of wave exciting slot. In the equa- 
torial plane 
f 
% H. .— Eee 5) 
a Vv bo/' €0 
and all other radiation fields are zero. 
Tabulations of the magnitude and phase 
of f(¢) for ¢=0, .e., of the on axis fields are 
given in Table I which also gives the number 
of terms taken in the summation of (131). 
‘It is believed that with this number of 
‘terms, five significant figure accuracy is 
‘obtained in both the magnitude and phase 
‘of f(¢). Computations* were made for 
1<c<30 in increments of 1 up to c=10 and 
in increments of 2 for 10<c<30, where 
¢=ka. Not all patterns for these cases are 
‘shown in Fig. 1 for the sake of clarity. The 
angular increments in ¢ were in 5 degrees 
for 0<¢<90 and 2 degrees for 90 <$<180. 
All the patterns are symmetrical with re- 
“spect to ¢. Tabulation of all these calcula- 
‘tions should appear in a forthcoming report. 
“From these calculated power patterns, it 
is seen that as the cylinder becomes larger 
_the patterns sharpen. 
A The equatorial plane radiation pattern 
for the case of a half wavelength slot in an 
_ infinite metallic plane is also shown in Fig. 1. 
For this slot 


7 
=. 


| A Pre; cos (5 sin ) . 


4 ar COS 
_ The radiation fields from this slot have been 
' obtained in the past by using the magnetic 
- current concept,®-° but can also be obtained 
"as a special case of the radiation fields pro- 
_ duced by an aperture in an infinite plane 
coated with a dielectric-magnetic slab.’ It is 
seen that the pattern for a large cylinder 
“approaches that from the plane. In par- 
ticular, for cylinders of size c2 30, the pat- 
- terns are identical within approximately 3 db 
for angles off the axis of +65 degrees. In 
addition, from Table I it is seen that for 
cylinders of c>30, the absolute magnitudes 
of the on axis field strength from the plane 
‘and from these large cylinders are identical 
_ within three significant figures. 


; C. M. KNop 
4 A. R. BATTISTA 


a Hallicrafters Co. 
; Chicago, Ill. 


3 Basically, the program for computation utilized 
~ the following methods: SHARE NUBES-3 and the 
technique of Goldstein and Thaler.‘ The technique 
employed in the electronic computation was to output 
the results after each sum, where a sum was formed 
- for each value of ¢ used. Each sum was then exam- 
ined, and when it was evident (as ascertained by in- 
spection of future sums) that any more summing 
would not affect the seventh significant figure of |f (?) | 
the process was stopped. The time required for the 
IBM-704 computer to process each case was approxi- 
mately 20 seconds to one minute depending on the 
size of the cylinder. This program is currently avail- 
aye Goldstein and R. M. Thaler, “Recurrence 
techniques for_the calculation of Bessel functions,” 
Mathematical Tables and Other Aids to Computation, 
vol. 13, pp. 102-108; tees 
5 Wait, op. cit., DD. : f 
6 y Dt Ck “Antennas,” ees Book 
., Inc., New York, N. Y., D- ; 1950. 
a 7C. M. Knop, “The Radiation Fields Produced 
by a Specified Tangential Electric Field Distribution 
Over an Aperture in a Perfectly Conducting Infinite 
Plane Coated with a Dielectric-Magnetic Slab,” Res. 
Lab., Bendix Systems Div., Bendix Corp., Ann Arbor, 
Mich., Research Note No. 30; February 1, 1961, 
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Dielectric Lens for Second-Mode 
Spiral* 


The operation of a double Archimedes 
spiral antenna in the second mode was de- 
scribed recently.! For a particular applica- 
tion its omnidirectional pattern and _hori- 
zontal polarization were suitable, but the 
vertical beamwidth was too large. An at- 
tempt was made to reduce the vertical beam- 
width while increasing the antenna gain. 

It had been suggested? that one type of 
horizontally polarized omnidirectional an- 
tenna of small vertical beamwidth could be 
obtained by embedding a small horizontal- 
loop antenna in the center of a flat disc of a 
dielectric material. This same idea has been 
applied to the second-mode spiral with the 
results that are given here. 

Our double spiral is similar to that of 
Donellan! except that the central feed sec- 
tion is a 4-inch-thick sandwich transmission 
line between two ground planes instead of 
being microstrip. The outer diameter of the 
spiral is 12 inches, and the ground plane is 
2 inch in diameter. The whole spiral is em- 
bedded in a dielectric material with e=2.56. 
The dielectric is one-half inch thick at the 
center and tapers to a line at 8 inches in 
diameter. The assembled view is given in 
Fig. 1 and the parts are shown in Fig. 2. 


Fig. 1. 


Fig. 2. 


The vertical radiation pattern of the 
antenna has a beamwidth of about 50° over 
the frequency range 4000-5600 Mc (essen- 
tially the bandwidth of the basic spiral). 
A typical pattern is shown in Fig. 3. By the 
application of the lens, the vertical beam- 
width has been reduced by a factor of about 
3, with a consequent increase in gain of be- 
tween 4 and 5 db. 


* Received by the PGAP, May 3, 1961. 

1J. R. Donellan, “Second-mode operation of the 
spiral antenna,” IRE TRANS. ON ANTENNAS AND 
PROPAGATION, vol. AP-8, p. 637; November, 1960. 

2C. Ducot, “Omnidirectional dielectric aerials for 
ultra short waves,” Le Recherche Aeronaut., no. 9, 
pp. 33-36; May-June, 1949, Ministry of Supply, 
Royal Aircraft Establ., Farnborough, Eng., RAE 
Translation No, 546, 
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Fig. 3—Antenna power pattern—vertical plane, 
no reflector. 


Fig. 4—Antenna power pattern—vertical plane, 
reflector spacing 1.25 inches. 


Incidental to the development of this an- 
tenna, patterns were also taken of the double 
spiral with dielectric lens mounted at dif- 
ferent spacings from a large ground plane. 
A spacing of 11 inches from the center of the 
spiral to the ground plane gave the pattern 
shown in Fig. 4. 

J. H. CRAVEN 

Radio and Elec. Engrg. Div. 
Natl. Research Council 
Ottawa, Can. 


Comments on “Diffraction of Scalar 
Waves by a Circular Aperture”* 


The authors would like to call attention 
to an error in sign and an omission of a cru- 
cial sentence in the above paper.1 Both the 
error and the omission occur in the unique- 
ness proof which appears on page S17, col- 
umn 1. As the proof stands, the reader may 
be left with the impression that the unique- 
ness of the solution follows from the equation 


J, (vole +] olpar = 0. 


* Received by the PGAP, May 4, 1961. 

1J, Bazer and A. Brown, Inst. Math. Sci., New 
York University, N. Y., IRE TRANS. ON ANTENNAS 
AND PROPAGATION, vol. AP-7, pp. S1 2-S20; December, 
1959, 
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However, this equation is incorrect and 
should, in fact, be 


if (| vol? — #| »|>av = 0. 
D 


To’conclude uniqueness from this corrected 
equation, it is sufficient to require that k 
have a small positive imaginary part and a 
sentence to this effect should have been in- 
cluded at the beginning of the proof. For a 
uniqueness proof which does not make use 
of this added assumption, we refer the 
reader to the forthcoming report by Levine.” 
J. BAZER 

A. BROWN 

Inst. Math. Sci. 

New York University, N. Y. 


2L. Levine, “Some Uniqueness Theorems for the 
Reduced Wave Equation,” Inst. Math. Sci., Div. of 
Electromagnetic Res., New York University, N. Y., 
Res. Rept. No. BR-33. 


On the Determination of the Disk 
Temperature and the Flux Density 
of a Radio Source Using High-Gain 
Antennas* 

INTRODUCTION 


One of the major observational tasks in 
radio astronomy is the determination of the 
disk temperature and flux density of radio 
sources. A highly directive, pencil beam an- 
tenna is often used in the observations. The 
power output available from the antenna is 
used to calculate the disk temperature and 
flux density of radio sources. When the 
beamwidth of the antenna is very large com- 
pared with the angular extent of a radio 
source under observation, the disk tempera- 
ture 7% and the flux density F of the radio 
source are related to the observed antenna 
temperature 74 by! 


IL AO 
Dies oe 1 
: Nr Qa ( ) 
and 
2kT 4 
F =a , 
ie (2) 


where Qu is the solid angle subtended by the 
source, Q4 is the antenna solid angle, k is 
Boltzmann’s constant and A,y is the maxi- 
mum effective aperture of the antenna. 
The antenna solid angle Q,4 is given by 


oh = if (6, ®)d2 = =, (3) 


M 


where f(0, ®) is the radiation power pattern 


* Received by the PGAP, May 18, 1961. This 
work was carried out under Contract AF 19(604)-4079 
of the AF Cambridge Res. Labs., Electronics Res. 
Directorate through The Ohio State University Res. 
Foundation. 

1jJ, L. Pawsey and R. N. Bracewell, “Radio As- 
es The Clarendon Press, Oxford, England; 
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of the antenna. normalized with respect to 
the maximum radiation intensity. Thus, 
f(@, 6)<1. Dy is the maximum directivity 
of the antenna. It is well known in antenna 
theory? that A-wQ4=7,\, where d is the 
wavelength and 7; is the radiation efficiency 
of the antenna. 

In order to detect weaker radio sources, 
an antenna of larger effective aperture must 
be used which also provides a narrower 
beam. As the beamwidth of the antenna 
becomes comparable with the diameter of 
the radio source, (1) and (2) are no longer 
valid. The purpose of this communication 
is to show the accurate relations which exist 
between the observed values and the physi- 
cal quantities of the radio source. Conversion 
factors are derived for three typical antenna 
patterns. These are useful in calculating 
the disk temperature and flux density of ex- 
tended radio sources from their observed 
antenna temperature. 


Disk TEMPERATURE AND FLUX DENSITY 


Let us consider a highly directive, pencil 
beam antenna pointing at the center of a 
radio source. It is assumed that the radio 
source is a circular disk with an angular di- 
ameter 62, having a uniform brightness tem- 
perature 7); across the disk. The radiation is 
randomly polarized and the source is located 
against a background of zero brightness. 
The antenna temperature due to the radio 
source is then given by 


ih J eee &)dQ 


tN —— Tn 
f f(, ®)da 
4r 
ff se, a0 
= mT Qa EO ae oe 
Q@Qa 
or 
IAD 
T,= A M4 Qa (4) 


Sen of 
Ge Pa f if (6, ®)d2 


The flux density of the radio source is 
given by 


2RT, 2kT. 
Fai| | = da = ——* ny. 
source 2 x? 


Substituting (4), we have 


ne 2RT 4Qa Qa 
nr 
{f feaa0 
source 
i 2kT 4 Q 


d 


Comparisons of (1) and (2) with (4) and 
(5) show that a correction factor inside the 
brackets in (4) and (5) is required to obtain 
the accurate disk temperature and flux 
density. This correction factor becomes 
unity when the source is so small that 


2J. D. Kraus, “Antennas,” McG -Hi 
Co,, Inc., New York, N, Y,; 1950 oC eee 


, 


ee 
e 

f(0, ®)=1 across the source. Then (4) and 
(5) reduce to (1) and (2). { 
The correction factor may be readily 
computed if the antenna power pattern and | 
the angular diameter of the radio source 
are known. In order to compute the correc- 
tion factor, it is sufficient to know the an- 
tenna power pattern for the regions of the H 
major lobe and its immediate vicinity. { 


\ 
CoRRECTION FACTORS 
Three mathematical functions are se- | 
lected for computing their correction fac- 
tors. These functions represent typical | 
shapes of the antenna major lobe encoun- © 
tered in practice. The correction factors are — 
derived in terms of the ratio of the angular 
diameter 6, of the source and the half-power _ 
beamwidth 6y of the antenna. The antenna 
beams are assumed to have a circular sym- _ 
metry with respect to its axis of maximum | 
directivity. f 
Case 1: Let f(0, 6)=e—®® , where p is a 
constant. If 6% is the angle between the 
half-power points of the antenna main lobe, 
f(0, ®)=4 at 0=64/2, from which 


=2,/In 2/éx. Thus we have 


ib i ___ Le, aaa 
afr 


T 
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6a /2 a 
[e-@0) | sin odode 


[1 ee (v8 /2)"| 


assuming sin @=6 for small angles. 
The correction factor K,; becomes* 


i Qa ___ (in 2)@a/On)? 
f f f(@, ®)dQ 


K = 
Case 2: Let 


1 [— en 20(64/4,)2] 


sin (p sin 6) 7? 


IGA [ psin@ : 


where p isa constant. At 6= 6/2, f(@, &) =4, 
from which we obtain p=1.392/sin (04/2) 
=2.784/0% assuming sin 6=@ for small 
angles. Thus we have 


J pee: (0, &)dQ 


2r 6/2 j 7 0 2 
= 1) it i EY amt) sin 6d6d®. 
0 0 psin@ 


Substituting w=p sin @ and sin 6=6, the 
above integral becomes 


Qr ¢ ?9q!? sin? uv aw (?94 1—cosu 
= du = — — du 
peo u pJo u 


= i Cin (p64), 


CE ILAe, Ko, “On the Analysis of Radio Astronomi- 
cal Observations Made With High Resolution Radio 
Telescope Antennas,” The Ohio State University, 
Columbus, Radio Observatory Rept. No, 21; Feb- 
Tuary, 1961, 


_ 


? 1 
1961 


where Cin (p62) is the modified cosine- 
integral whose values are tabulated in the 
literature.?4> 
‘ The correction factor Kz becomes 
K,= Ee ee 
a[Cin (p64)]/p? 
__(2.78404/62)? 
4 Cin (2.78462/6z) 


Case 3: Let 
a 70,®) = [ 


where J; is the Bessel function of the first 
ind. of order one and # is a constant. At 
6=6n/2, f(0@, &)=4, from which we obtain 
p=1.616/sin (01/2)=2X1.616/6%. Substi- 
tuting w=? sin @ and sin 6=8, we have 


2Ji(p sin 6) 7? 
p sin @ 


ff se, a0 


4 for pve" 47) 
ee i 

oa 8a  Pq'? J1?(u) 

: = — SS SS du 

ec . peo ud 

Be 4a 

= ro [1 — J:2(p0a/2) — Jo*(p04/2)], 


« 


~ where J; and Jy are the Bessel functions of 
- the first kind of order one and zero, respec- 
tively. 


“a Thus the correction factor Ks becomes 
4 Q 

K3 = hl da 

© ff sone 

4 source 

2 Pe Qap? 

© Salt = F.2(p04/2) — Jo2(004/2)] 

or (1.61602/6x)* 
; 4[1 — J:2(1.616 02/ox) — Jo2(1.616 62/0x)| 


It can be shown that when 64/6H be- 
_ comes very small, all these factors reduce to 
~ unity as expected. The order of corrections 
' required may be seen by the following nu- 
merical examples. Thus, for 04=6z, Ky 
P =1.39, Ke=1.38, and Ks=1.38; for 64 
_ =26n, Ki=2.96, K,=3.19, and K;=3.15. 

j H. C. Ko 
Radio Observatory 
3 Dept. of Elec. Engrg. 
‘ The Ohio State University 
i Columbus, Ohio 


4C, T. Tai, “Table of Modified Cosine-Integral,” 
Stanford Research Institute, Menlo Park, Calif., 
Tech. Rept. No. 15, SRI Project No. 188; January, 


1951. : 
5 “Tables of Generalized Sine- 


= and Cosine-In- 
 tegrals,” Harvard University Press, 


Cambridge, 
Mass.; 1949. 


Operators for Wave Equations in 
_ General Linear Media* 


This communication presents some linear 
vector and scalar operators for electromag- 
netic waves in stationary, 1sotropic, and 


* Received by the PGAP, May 19, 1961. 
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linear media. Commutability of two linear 
operators are discussed in deriving wave 
equations. These waves are directly obtained 
from a set of Maxwell’s equations: 


curl E(q, t) = LiH(q, 0), (1) 
curl A(g, 2) = J, t) + LE, 4), (2) 
E-grad e(q, t) + e(q, t) div E = p(g,é) (3) 
A-gradu(g,t) tug) divH =0, (4) 


where linear operators Z; and Lz stand for 
—p:i—u(d/dt) and e+ ¢(0/dt); the subscript 
t stands for the partial derivative with re- 
spect to time; the symbol qg designates the 
magnitude of a general position vector @; 
the current density vector J can be either 
o(q, t)E or pi(q, t), where oa is a linear iso- 
tropic conductivity and @ is the electron 
velocity. The four equations above are com- 
bined to give 


H1H-O. 6 
bn Ln ) 3): . 
where L. is V?2—(em):—[em+(eu):](0/0t) 
—pe(d?/at?)+V(Ve-/e) and the last_ term 
produces (Ve-E/e). J, is used for —(VuX )s, 
and source terms S, and Sm, respectively, 
represent grad(p/e)+(uJ): and —curl J. 
Indicated interchanges between ¢ and yp» 
provide Lm and Jp, such as Tom = Le ep, 
pe) and Im=1.(u—e), respectively. 


If may be hard, if not impossible, to ob- 
tain a matrix of inverse operators so that 


Ne ine 1 
eee Ga). 5° 
Nm Nin In ia (Oi 

If this N matrix is determined, the pre- 
multiplication by N matrix gives E and H 
explicitly. This process has been shown 
later for restricted cases. In the past, the 
identical forms for linear scalar operators, 
L, and Im, have been widely treated for 
homogeneous media. However, other media 
can assume the identical form for ZL, and 
Lm; this is realized if (¢/u); is zero. 
» The restricted class of source-free (0 =0, 
7 =0) physical systems assumes the propa- 
gation of waves in the z direction, as along 
a cylinder parallel to the z axis having ar- 
bitrary cross-sectional geometry. Let the 
subscripts s and z stand for the transverse 
and z-direction parts. Eqs, (1) and (2) are 
rewritten as 


Li Vix ie) =v,X ek (7) 
—VzX I» Es Hz 

The problem here is how to derive linear vec- 
tor wave equations for H.(q, t) and E.(q, t) 
when the z-direction field components are 
given in general forms as Eq, t) and 
Hq, t). This problem will be discussed for 
four medium types separately. 


A homogeneous medium, 1.€., € and pw are 
independent of q and t 

The inverse matrix is readily obtained, 
giving! 

1W. K. H. Panofsky and M. Phillips, “Classical 


Electricity and Magnetism,” Addison-Wesley Pub- 
lishing Co., Inc., Reading, Mass., pp. 192-193; 1955. 
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E, 
n(z) 
A. 


to) . 
—VzXx Vex BVsX% ‘Ot Ez 
= i he) 


(6) = 
—aoVs pe SAE 8 z 
o TR VzeXVsxX / \H, 


where Ls is for {V-XV.X +eu(d?/a*)}. 
A Medium of &t) and p(t) 


In this case the NV matrix for (7) is ob- 
tained, resulting in 


Be Vix Vex LiVsX E, 
Tea ee NEES 

Hs DPVsX VzeVsX/\Hz 
where Ls is for {V-XV2X —LiL}. 


A medium of &q) and w(q), and a medium of 
e(q, t) and u(q, t) 


Equations similar to (8) and (9) are not 
easily derived. Because operators V, and 
L; do not commute each other, 


(LiveX — VeXL1)He = (usoXHe)4, (10) 


where pu, designates the partial derivative of 
p with respect to z and is a unit position 
vector. For this reason E;(q, ¢) is not deriv- 
able from general functions, Eq, #) and 
Hq, t). Similarly, the noncommutability 
between V, and Zz is shown by 


(L2V2% sa VX Lo) Es =. (eGoX E,) t- (11) 


An Exceptional Case of <q, t) and y(q, t) 
for Harmonic Fields 


If Z, and H, vary with a specified time- 
dependence, ¢.g., a harmonic variation of 
exp(jwt), the separability of field quantities 
between the time-dependent and the posi- 
tion-dependent parts requires that the time 
dependence of ¢(q, ¢) and u(q, ¢) be the same 
as the time-dependence of E, and H.. Sup- 
pose £, and H; vary like exp (jwt); the in- 
verted matrix of this exceptional case pre- 
sents 


EQ 
Le H.(q) 


1 
—V:X—VeX —2jwoVsX\ / E:(g) 
3 i : , 2) 
2jwVeX VeX —VeX/ \A-(Q) 
€ 


where Ls is for an operator (VzX1/pV.X 
—4w%e) and Le=Ls(e—u, ue). 

The association between the commutator 
and the simultaneous measurements of two 
physical quantities corresponding to two 
linear operators has a sound basis for the 
quantum mechanics, whereas macroscopic 
experimental facts (Maxwell’s equations) 
yield physical quantities which are always 
measured simultaneously. Eqs. (10) and (11) 
only say that if wu, and ¢, are nonvanishing, 
linear wave equations for E, and H, are not 
simply related to E.(q, t) and Hq, t), as in 
(8) and (9). 

A few things taken for granted in 
homogeneous media are not directly ap- 
plicable to other medium types. For in- 
stance, Ohm’s law and Maxwell’s equations 
are valid in isotropic linear media, yet the 
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relaxation time of the medium, which is al 
combined result of the law and equations, 
is subject to examination, as in: 


(dp/dt) + (op/e) + D-grad (c/e) = 0. (13) 


Hoffman? studied the case of o(q), €(q), and 
a harmonic field for £. 

Phenomena of a nonlinear system in- 
clude those of its linearized system as a 
special case. Physical facts encountered in 
electromagnetic wave propagation through 
media of ¢(g, t) and: p(q, ¢) embrace many 
specialized cases, such as equations* of 
Mathieu for periodic variations of medium 
parameters and nonuniform transmission 
line problems.*:5 Although Maxwell’s equa- 
tions may not be separable in general, those 
interested in problems of satellite com- 
munications and controlled inhomogeneous 
devices (like plasma-filled waveguides) 
would find useful aproximate numerical so- 
lutions. Vector linear operators presented 
here may also be of some use in these prob- 
lems. 

I. SuGAI 
ITT Federal Labs. 
Nutley, N. J. 


2 W. C. Hoffman, “A possible mechanism for ra- 
diation and reflection from ionized gas clouds,” PRoc. 
IRE, (Correspondence), vol. 47, pp. 1274-1275; July, 
1959, 

3A, M. Markovkin, “Some properties of media 
with periodically time-varying parameters « and uy,” 
Radiotekh. i Elektron, vol. 4, pp. 785-791; May, 1959. 

4H. Osterberg, “Propagation of plane electro- 
magnetic waves in inhomogeneous media,” J. Opt. 
Soc. Am., vol. 48, pp. 513-521; August, 1958. 

5L. M. Brekhovskikh, “Waves in Layered Me- 
dia,” Academic Press, Inc., New York, N. Y., pp. 
168-189; 1960. 


On the Transposed Radiating 
Systems in an Anisotropic Medium* 


It has been known for some time!” that 
the Lorentz form of reciprocal theorem, or 
the Rayleigh-Carson form of reciprocal 
theorem, can be generalized to two sets of 
fields in an anisotropic medium provided 
that one of the fields is defined in a medium 
which is transposed to the other. Let us con- 
sider, for example, the case of a magneto- 
ionic medium. with a dielectric tensor 
medium represented by [e]. If we denote the 
two sets of fields by (A,, Hi) and (Ey, H2) and 
assume that they are harmonically varying 


* Received by the PGAP, May 29, 1961. This re- 
search was supported in part by Aeronautical Sys- 
tems Div., Air Force Systems Command, Wright-Pat- 
terson Air Force Base, Ohio, under Contract No. 
AF 33 (616)6782, and The Ohio State University 
Research Foundation. 

1M. H. Cohen, “Reaction Concept for Electro- 
magnetic Field in an Anisotropic Medium,” The 
Ohio State University, Columbus, Antenna Lab. In- 
formal Rept., 1954; see also V. H. Rumsey, “Reaction 
concept in electromagnetic theory,” Phys, Rev., vol. 
94, pp. 1483-1491; June, 1954. 

2R. F. Harrington and A. T. Villeneuve, “Re- 
ciprocal Relationships for Gyrotropic Media,” IRE 
TRANS. ON MICROWAVE THEORY AND TECHNIQUES, 
vol. 6, pp. 308-310; July, 1958. 
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with respect to time, then they satisfy 


VX A = — jou (1) 

VX HM = it jolelE (2) 
and 

VX Ep = — jou (3) 

VX Mp = I +jole|Ey, (4) 


where [é] denotes the transposed tensor of 
[e]. The currents J; and J2 are the ones not 
of magneto-ionic origin. Applying the tech- 
nique used in deriving the ordinary Lorentz 
theorem one finds that 


ff [vex me — Bx Hilde 
Vv 


-[ffie- Jy E;)dv (5) 
f[@xm- E, X Hh) -dS 


=fffam- Jy-E;)dv. (6) 


If in the region of integration J; and J2 are 
zero then one obtains the Lorentz form of 
the reciprocal theorem, that is, 


or 


Jf [ext Bx m)-as =0. (7) 


If we let S go to infinity in (6), then the 
Rayleigh-Carson form of the reciprocal the- 
orem is obtained, that is, 


fJfGe- Iy-Ejdo=0. ©) 


Eq. (8), when interpreted in a more sophis- 
ticated manner, forms the basis of Rumsey’s 
reaction concept. 

We apply now (8) to the problem of re- 
lating the transmitting and the receiving 
properties of an antenna. Assume that the 
axis of the medium, corresponding to the 
direction of the stationary magnetic field 
Bo, is in the direction of the z-axis as shown 
in Fig. 1(a). Let (E,, H,) be the field created 
by a transmitting antenna 74 with the re- 
ceiving antenna Rg placed at the central 
region of the coordinate system. The receiv- 
ing antenna Rg is represented schematically 
by three heavy lines. For field (Ey, He), we 
interchange the role of transmitting and 
receiving as shown in Fig. 1(b). The sta- 
tionary magnetic field is reversed so that 
the medium would correspond to the trans- 
posed tensor [é]. Applying (8), one obtains 


Vale = Vala. (9) 


Eq. (9) is the reciprocal relation as applied 
to two antennas in an anisotropic medium. 
So far we have merely repeated the exposi- 
tions which are already known. 

Now suppose that the Bo field in Fig. 
1(a) is a fixed field which could not be re- 
versed, the system in Fig. 1(b) would then 
be impractical to reproduce. In order to 
remedy this difficulty a transposed radiating 
system can be formed that is shown in Fig. 
1(c). The system is created by rotating the 
entire structure in Fig. 1(b), including the 
Bo field, around the x-axis by 180°. Electro- 
magnetically, the transposed radiating sys- 


(a) (b) 


Fig. 1—(a) A transmitting antenna and a receiving © 
antenna in an anisotropic medium. (b) Therole 
of transmitting and receiving being changed as 
well as the direction of the Bo field. (c) The trans- 
posed radiating system with the same Bp field as 
in (a). 
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(c) 
Fig. 2—A more complicated system where (a), (b), 


and (c) have the same significance as those in 
aT ae 


tem has the same property as the one before 
rotation. The medium in Fig. 1(c), however, 
is not transposed. In fact it has the same 
characteristic as the one in Fig. 1(a) if [e] is 
symmetric with respect to the x-y plane. 
The reciprocal relation as stated by (9) now 
can be interpreted as the one between Fig. 
1(a) and 1(c). The above demonstration can 
be illustrated by a more complicated system 
shown in Fig. 2, where the circle represents 
a spherical body upon which an antenna B 
is mounted. It is obvious that the receiving 
property of B is now directly related to the 
transmitting property of the same antenna 
mounted at the opposite side of the sphere. 
The application of the transposed radi- 
ating systems to some specific problems will 
be reported later in a subsequent paper. The 
author acknowledges the valuable discus- 
sions which he had with Dr. R. G. Kou- 
youmjian and Dr. H. C. Ko on this problem. 
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